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INTRODUCTION 


I AM DEEPLY HONORED to have been chosen to deliver 
this Seventeenth Wright Brothers Lecture before the 
Institute of the Aeronautical Sciences on this Fiftieth 
Anniversary of Powered Flight. It was my privilege to 
present the Wilbur Wright Memorial Lecture before 
The Royal Aeronautical Society in London in 1931, and 
I am indeed proud to be numbered among those few 
men who have been twice honored in this series. 

It is a humbling experience to face the task of review 
ing, however briefly, the incredible story of the first 
half century of flight in America. 

I am happy to assay this task, however, because of the 
deep sense of obligation I feel toward aviation in 
America. It has been the central theme of my life for 
a half century, and it has been the most demanding yet 
the most rewarding of any field of work I could possibly 
have chosen. It has enriched my life immeasurably in 
friendships, in honors and awards, and in that priceless 
feeling of successive jobs well done. Certainly no man 
can ask more of his chosen profession. 

I hope that I have made some contribution to it, no 
matter how small, out of the depths of that obligation. 
If I have, however, it has been only as a member of a 
continuing team of enthusiastic and loyal coworkers 
who have shared with me the triumphs and the failures 
of these interesting years. It is my hope, too, to make 
contributions to it in the future, for I have no intention 
of ceasing my labors in aviation in America. 


THE WRIGHT BROTHERS 


It has been said that nothing in this first half of the 
Twentieth Century the atomic bomb 
transcends in importance the magnificent accomplish- 


not even 
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MARTIN 


half 


Having had some experience with air 


ment of the Wright brothers at Kitty Hawk a 
century ago. 
planes and none with the atomic bomb, I am not pre 
pared to dispute this high tribute to these famous 
brothers. 

The death of Wilbur Wright in 1912 robbed our in 
dustry of its greatest pioneer, and I have no question 
that we have the Fortunately, 
Orville Wright survived his brother by 35 years, and it 
was he who discharged the very great responsibility of 
holding the Wright name aloft through periods of 
severe trial] and who preserved its dignity and integrity 


been poorer for it. 


in the world’s esteem. 

I was closely acquainted with Orville Wright, for 
many years, first as friend and later as business associate 
and I can say that he was the most honest man I ever 
convinced that it 
that lay at the root of the Wright 
It has been 


met. I am was this devotion to 


truth—in all things 
brothers’ success in creating the airplane. 
honest men who have built our great industry and who 
will guide its sound growth in the future. 

The quality of honesty in daily dealings is an infallible 
guide to a man’s capabilities in the engineering and 
scientific fields. Structures and machines are unfor 
giving of the cheater and inevitably indict those who 
toy with the facts. It was their indomitable faith in 
truth that guided the Wright brothers through all their 
investigations and_ that perhaps, 
brought them to success at Kitty Hawk 50 


ultimately——and, 
inevitably 
vears ago. 
It would be a vain man, indeed, who would attempt 
to add to the stature of the Wright brothers and the 
technical magnificence of their 1903 Flyer. You are 
all familiar with the fact that this first airplane weighed 
750 Ibs. fully loaded yet was lifted into the air by an 
engine developing only 12 hp. This phenomenal power 
loading meant that each horsepower was lifting 62? 


Ibs. of weight, 
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Their achievement was more than merely lifting 
themselves from the ground, for as we study their 
what truly 


have learned 


You are familiar, I am 


writings and records we 
scientific workers they were. 
sure, with their search for suitable testing methods 
during which they developed a number of ingenious 
devices for getting at the scientific truths of flight. 
Their wind tunnel, for exampie, used a honeycomb grid 
for straightening the airflow, and in its 27-m.p.h. throat 
they tested model wings that varied over a wide range 
of aspect ratio, camber, and taper. They made more 
than 1,000 glider flights at Kitty Hawk prior to their 
first powered machine. Unquestionably they read all 
of the written records of aeronautical thought available 
tothem. Thus, they gave us not only the first airplane 
but an entire scientific method for work in the field of 


aeronautics. 
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Production history of the U.S. aircraft manufacturing industry during its first 50 years. 


But the Wrights did not work in a scientific wilder- 
They stood on the foundation erected by their 
Yet all of these pearls did 


ness. 
predecessors in science. 
not make a necklace until the inspired jewelers came 
along. 

The Wrights 
approach to invention in which intuition and logic are 


brought to aviation the scientific 
used as a leavening for the diverse technical facts at 
hand. 
century prior to the Wrights, but it was their course ol 


This method, of course, had prevailed for ¢ 


action following a conclusion that made them im 
mortals and their predecessors only dreamers and vain 
experimenters. Once having fashioned an idea, the 
Wrights proceeded to test it in scale form, either in their 
wind tunnel at Dayton or in a kite or glider at Kitty 
Hawk. no matter how sound or 


complete the idea, they were prepared from the start t 
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FIRST HALF-CENTURY 


modify it, to improve it, or to abandon it. For them 
there was no perfection. 

rhe real legacy of the Wright brothers to their 
‘ontemporaries was the single truth that flight had 
been achieved. With this fact in hand, the rest of us 
were able to sweep all difficulties before us and rush to 
our individual successes. It is the haunting fear that 
what he is trying to do is impossible that corrodes the 
will of the inventor. Once this fear is removed, the 
spirit overrides all problems and moves quickly to the 
fruition of the dream. 

All of us who followed the Wrights lived beyond such 
fears. Our only problem was how best to do the job. 
Freed of this worry, Glenn Curtiss on the East Coast 
sped on to success with an advanced design and I 
on the West Coast, built, taught 
myself to fly my first airplane without the benefit of 
ever having seen one. When I started, I did not know 
how to do the job; I only knew that the job could be 
But it was this single bit of information that 


designed, and 


done. 
dwarfed in importance all the other knowledge at my 
disposal in its contribution to my success. 

I have had many years in which to recall my first 
flight and to attempt to express the experience in 
words, but I have always failed. Leaving the ground 
for the first time in a machine of your own design, which 
you have built with your own hands, is simply an in- 
describable experience. One of my genuine regrets, in 
the face of the miraculous progress we have made in 
aviation in the past half century, is that the exhiliration, 
pride, humility, fear, and triumph of such an experience 
now seems forever barred from human experience. 

The Wrights unleashed a torrent of thought which 
even today shows no signs of abating. Once we had 
solved the initial problem of simply getting ourselves 
into the air, we turned to the great debate on the form 
of the airplane. I have often felt that aviation’s most 
inventive period was its first. It was the period of the 
grand idea as the airplane quickly evolved from its 
original form in the 19083 Wright Flyer to the basic form 


in use today. 


BEFORE WORLD War I 


The 1911 
books as the age of the death-defying exhibition pilot, 


1912 period is described in the history 


and in those two brief years such names as Lincoln 
Arch Hoxsey, Ralph Johnstone, Phil 
Parmalee, Frank Coffin, Howard Gill, J. C. “Bud” 
Mars, ““Becky’’ Havens, Charles Willard, and a host 
of others became immortal. But while their exhibitions 
provided thrills for the crowds and a fabulous revenue 
for the Wright and Curtiss Exhibition companies, these 


Beachey, 


pioneer pilots were daily creating flight-test data on 
their machines and amassing important performance 
stability, and control experience for the designers. 
During every exhibition flight I ever made, I was alert 
to the performance and handling characteristics of my 
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airplane and was making mental notes on its short 
comings and the possibilities of their improvement. 
Aeronautical science was being rapidly evolved, even 
while ‘ladies fainted and strong men quelled.’ 

Until World War I aircraft manufacture consisted of 
building one airplane and, after that one was com 
pleted, transferring the personnel to start construction 
on the second, and soon. When I received a large order 
for the Model TT airplanes from the Signal Corps, it 
was perfectly clear that our way of doing things would 
have to change. We therefore adopted the new prin 
ciple of making all of the parts first and assembling the 
This 


led, of course, to the assignment of workmen to specific 


airplanes as they moved along an assembly line. 


tasks, such as rib-making or engine mount fabrication, 
with no other duties. Our factory floor had to be laid 
out into separate departments and the work performed 
on a predetermined schedule. We prepared and used 
rudimentary flow charts for the first time and were able 
to pinpoint delays by determining the department in 
difficulty and ascertaining the preceding department 
from the chart, which spotted the trouble quickly. 


MATERIALS 


Wood, of course, was our staple material in those 
days and I will certainly not defend it as a material of 
aircraft construction. Its principal difficulty was its 
variation in quality, almost from one board to the next. 
This naturally made it impossible to establish firm 
allowables for design or manufacture. Its strength 
characteristics literally changed with the weather, but 
we always believed it to be in excess of the required 
strength. 

However, we soon learned many effective ways of 
treating, storing, and handling woods and veneers, and 
a visit to any major museum today will reveal wooden 
airplanes 30 or 40 years old still in good states of 
preservation. 

The introduction of sheet metal into aircraft con- 
struction in the late twenties turned a great adventure 
into a science and revolutionized not only aircraft 
manufacture but design as well. We had been building 
metal airplanes for some years, of course, in the form of 
welded steel tubing fuselages introduced during World 
War I. The first high-strength aluminum alloy, 175, 
was made and tested early in 1917, and it was first 
for construction of 
Naval Aircraft 
Our first use of it in air- 


produced commercially in the U.S. 
the airship “Shenandoah” at the 
Factory during 1920-1921. 
craft construction was in the fabrication of formed ribs 
and spars, which, when combined with the welded steel 
tubing fuselage, produced an all-metal airplane struc 
ture fabric-covered. 

The Goodyear Tire & Rubber Company conducted 
an extensive research program on the forming and heat- 
treating of 17S in preparation for construction of the 
RS-1 dirigible, and I engaged their services for the 
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fabrication of parts for the Martin MS-1, the first all- 


metal airplane in America. 


METAL AIRCRAFT 


The cantilever wing marked another great revolution 
in aircraft design and manufacture. This problem had 
long been one simply of choice rather than possibility. 
A wing had to be thick in order to permit sufficient 
structure for internal bracing, and thick wings, of 
course, were high-drag wings. However, through wind- 
tunnel research we finally developed a profile thick 
enough to permit internal bracing, yet with acceptable 
drag characteristics. With these profiles available in 
the late twenties, we faced the problem of fabricating 
stress-skin structures in production. 

It required development of a whole new science in 
manufacture. Heat-treating equipment, metal forming 
and cutting machines, and the entire new science of 
tooling had to be developed before mass production of 
metal aircraft was possible. So revolutionary was this 
new development that I solved the problem only by the 
layout and construction of an entirely new factory 
especially designed and equipped for metal aircraft 
fabrication. 

With the all-metal cantilever monoplane and the 
semimonocoque fuselage, designers were able to make 
their calculations with great precision and aircraft per- 
formance moved ahead rapidly. I think it was at this 
period, in the early thirties, that science assumed the 
leadership over experiment and permitted aircraft 
design, for the first time in its history, to move forward 
to its full potentialities. 


AERODYNAMICS 


The man that founded our whole modern body of 
aerodynamic theory was Dr. Ludwig Prandtl, whose 
classic mind was lost to us only a few months ago. 
Prandtl presented his concept of the boundary layer in 
1904, which clarified the puzzling relation between the 
classical hydrodynamics of perfect fluids and real air- 
flows. The application of perfect fluid theory to many 
important practical aerodynamic problems resulted 
from this clarification of the role of viscosity. Con- 
tinued research on the boundary layer by Prandtl and 
his associates at Géttingen resulted by 1920 in a 
qualitative understanding of skin friction, form drag, 
and stalling. 

However, Prandtl’s works did not become available 
in America until after World War I, and our authority 
in the early days was Frederick W. Lanchester, one of 
the great geniuses of aviation. His classic textbook, 
Aeronautics, appeared in 1907 and presented his concept 
of ‘‘circulation’”’ and the existence of vortices at the tips 
of wings. I must admit that all of us found Lanchester 
very hard going, since he used very little mathematics 
in his work and relied largely on page after page of 
laborious discussion. 


AERONAUTICAL 
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Succeeding papers by Prandtl and Lanchester between 
1907 and 1919 developed the theory of the finite wing, 
the thing we were all interested in, and explained 
scientifically what we all had known from experience 
that the best lift-drag characteristics are obtained from 
a high aspect ratio wing. The Lanchester-Prandtl 
theory, as some call it, provided the first rational basis 
for the aerodynamic design of wings and also gave load- 
It also 
clarified the aerodynamics of biplanes and explained the 


ing information needed for structural design. 


existence of interference between the two wings, which 
we had begun to suspect and which led, eventually, to 
the monoplane. One of the important results of this 
theory was the fact that it broke wing design down into 
two separate parts: the effect of air foil section and the 
separate effect of wing geometry. 


AIRCRAFT CONSTRUCTION PROBLEMS 


The choice of materials for aircraft construction has 
never been a simple one because it lies at the root of 
nearly everything that is done in the manufacture of 
aircraft. The stress analysts, usually enthusiastically 
supported by the aerodynamics people, have always 
rushed to embrace a new and always promising aircraft 
material, and the fervor of their arguments has often 
severely tested the prudence of management. I have 
always admired this enthusiasm and this constant 
eagerness for progress, but the introduction of a new 
material in aircraft construction is a long and extremely 
costly business that must be considered from every 
possible angle. Usually, such materials require entirely 
new machinery, complex handling and processing, and 
invariably a substantial increase in costs in the early 
stages. Titanium, of course, is a current example of 
these promises and problems. 

But progress must never be denied, and experience 
has shown the long-range economy of developing and 
using the very highest quality materials and fabrication 
techniques available. This is because of the unique 
structural requirements of the airplane in which its 
materials must be loaded just as near their limits as 
seasoned judgment will permit in order that structural 
efficiency be raised as high as possible. 

It has always been wise to bear in mind that the 
actual costs of the materials in an aircraft structure is 
only 10-20 per cent of its total cost, the balance being 
in labor. Therefore, it has always been essential to 
choose a material that could be worked easily and 
rapidly, and it is good economy to pay twice as much for 
a new material if it reduces fabrication time by one- 
half or more. Unfortunately, this has rarely been the 
case, and, instead, we have attacked the problem of 
fabrication time by the more difficult approach of 
developing special tools and procedures. 

Again, however, this approach requires assessment of 
the probable number of aircraft to be built of a new 
model, the extent of tooling complexity and cost vary- 
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Fic. 4. Variation of wing loading with time for five basic types of U.S. aircraft. 
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FIRST HALF-CENTURY 

ing directly with the size of the anticipated order. 
Here is the crux of the aircraft manufacturing economic 
problem and the source of the greatest demands on 
management acumen. It seems reasonable to assume 
that about SO per cent of the structure of a new airplane 
model will not change throughout its life and that tool- 
ing can safely be ordered for such parts. Despite the 
growth potential and rapid development of a new model 
in its early years, these things actually effect only about 
20 per cent of the major assemblies, and these units may 


properly be lightly tooled. 


AIRCRAFT ENGINES 


The airplane, from its inception, has grown in 
capacity and performance only at the pace permitted 
by the aircraft engine, and we have historically selected 
the engine first and designed our airplane around it. 
The search for high power has never relaxed for an in- 
stant. After their flights in December, 1903, at Kitty 
Hawk, the Wrights returned to Dayton and im- 
mediately began the design of a more powerful engine. 
Aiter Curtiss made his first flights in 1908 with a 40-hp. 
Curtiss engine, he immediately began construction of a 
50-hp. engine for his next machine. In 1909 I made my 
first flight with a modified Ford engine of 12 hp. and 
wasted little time securing an Eldridge marine engine of 
30 hp. Thus, from the very first, we have energetically 
pursued more engine power and our magnificent air- 
craft engine, propeller, and accessory industries are 
keeping pace with our demands just as well as science 
and time permit. 

The subject of the turbojet and turboprop forms of 
aircraft propulsion are well known to you all, and I will 
touch on them only to express my conviction that we 
can now have all the power we demand, since the air- 
craft gas-turbine form has no inherent limitations to its 
size or Output. The rapid maturity of the turbojet in 
its few brief years of development has been one of the 
truly amazing stories of this half century of flight. We 
now have turbojet engines with reliability and dur- 
ability comparable to our piston-engined types and, 
when flown at sufficiently high speeds, economy closely 
approaching that of the piston engine. The rocket and 
ram-jet forms of propulsion know utterly no limit to 
their output, and the aircraft designer has had the 
barrier of power removed from his list of major 
problems. 

Today, on this 50th Anniversary of Powered Flight 
in America, we have flown at 1,655 m.p.h. and to an 
altitude of 89,000 ft. in man-carrying aircraft. We have 
lifted a total disposable load of 200,000 Ibs., and we have 
carried a load of 10,000 Ibs. for a distance of more than 
10,000 miles. We have flown around the world non- 
stop, we have carried 308 people in an airplane, and 
we have flown a distance of 11,236 miles nonstop with- 
out refueling. Only a few days ago we witnessed the 
inauguration of scheduled nonstop transcontinental air 
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service in less than 8 hours. Today we have a combat 
airplane in quantity production that flies regularly and 
smoothly at supersonic speed in level flight. The list 
of such miracles is endless. 


CHARTS SHOWING TECHNICAL PROGRESS 


Fig. 1 shows the production history of our industry 
during its first 50 years. I specifically had this plotted 
on simple graph paper, rather than the more popular 
semilogarithmic paper usually used, in order to illustrate 
the really furious upheavals in our annual rates of 
production. We see here not only the familiar peaks 
and valleys in our military output, which is geared to 
wars, but the identical slopes and proportional magni- 
I think that here, in 


this simple graph, is portrayed the fundamental problem 


tudes in commercial production. 


of the aircraft manufacturing industry and the root of 
nearly all of it difficulties. 

Fig. 2 is the most popular index to our technical 
progress: the top speed of our aircraft. Although it is 
customary to pinpoint all of the speeds and then lay 
down a “‘trend’’ curve through a sea of such dots, I 
have used here the precise maximum speed of each of 
the vears for each of the various aircraft types, since, in 
my opinion, “trends” are often misleading and the 
state of the art must always be judged by the best 
that we have—not the average. 

Fig. 3 illustrates our engine progress over 50 years. 
It shows how the various engine types have progressed 
and tends to indicate that a radically different power- 
plant type has been born just as the previous type was 
beginning to show its age. It might follow that nuclear 
energy power plants will be available when the turbojet 
curve starts to flatten. I think it significant, too, that 
our rate of development of a new engine type has always 
been much more rapid than the last. At the upper right 
corner, for example, we see that the flick of an after- 
burner switch in our currently most powerful turbojet 
engine instantaneously produces the same amount of 
power it took us 40 years to obtain with the piston 
engine. 

So far we have seen measurements of our absolute 
progress, but now let us see some qualitative data on 
that progress. Fig. 4 shows what has happened to 
wing loading in 50 years. This chart shows its in- 
evitable and rapid increase through the years, regardless 
of airplane type. Only the horizontal or time displace- 
ment is different for the different types, each following 
the other within a period of 5 years or less. One point 
of interest revealed here is that it is the bomber 
that has led the way in higher wing loading 
Today's jet bomber wing load- 


not 
the fighter 
for the past 20 years. 
ing, aS a matter of fact, is substantially higher than 
even our supersonic rocket-powered research aircraft. 
Our postwar transports have far exceeded even our 
fastest racing planes of prewar years in wing load- 


ing. 
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Fig. 5 shows the most consistency in power loading of 
any of the parameters examined here. It is seen that 
the transport and the fighter have decreased only 
slightly in power loading during the past 25 vears. 
This is probably explained by the equal consistency with 
which we match a power increase in a new aircraft with 
an almost proportional weight increase. Only the 
bomber continues to show improvement in power load 
ing, reflecting its multiengined character. With the 
modern pod arrangement we have reached a point where 
we can literally “throw in another engine or two”’ when- 
ever the gross weight shows signs of getting out of hand. 
For those purists who might quarrel with the apparent 
interchangeability of horsepower and pounds of thrust 
in these charts, I might point out that, right or wrong, 
it certainly simplifies chart-making. 

Fig. 6 illustrates a subject that has always been of 


paramount importance to me. Early in my career I 
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Increase in aircraft useful load and gross weight during the first half century of flight. 





came to the conclusion that only the airplane that could | 


lift a good load was of any long-range importance. | 
This figure shows just how much weight we have been 
able to lift and is, in my opinion, the most important | 
index of our technical progress. Only quantity produc- 
tion airplanes have been used in this study, and none 
of the experimental ‘‘giant’’ airplanes of the past has 
been included. The dark hatched area in the lower 
right portion of the figure shows that we have managed 
to lift well over 100 tons in useful load. 

This research work was performed by Robert Mc- 
Larren, the Editor of Aero Digest magazine, so we cat 
be assured of its accuracy. Like all engineers, Bob| 
wanted to plot these things on fancy logarithmic paper, 
but I am old-fashioned enough to prefer simple graph 
paper. When something is twice as big as something! 
else, I want to see it that way and not only one and one- 


third as big. 
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FIRST HALF-CENTURY 
These are the dramatic accomplishments of the first 
half century of flight in America. More important, | 
think is our realization that these are only the mani- 
festations of our progress and not progress itself. I 
have traced briefly the progress of man’s thoughts on 
aviation, rather than his deeds, for it is our thinking 
that builds the lasting foundations for continued future 


progress. 


DISCUSSION 


In a free independent society there can be no grand 
guiding hand controlling the destiny of an industry or a 
people. We in the aircraft industry have guided our 
own progress through our faith in, and our loyalty to, 
each other. The modern airplane constitutes an 
integration of the thoughts and lifetime work of more 
individuals than any other industry on earth. This 
integration was accomplished not by fiat but by under- 
standing cooperation between its many and diverse 
elements. 

This interplay of really diverse avenues of human 
effort is a monument to man’s dependence upon other 
men, and it has been only the unfettered freedom of 
America which has created this giant team of co- 
operating individuals which is our modern aircraft in- 
Today there appears to be no limit to the 
We have the aerodynamic knowl- 


dustry. 
power of this team. 
edge, the structural materials, the power plants, and 
the manufacturing capacities to perform any con- 
ceivable miracle in aviation. 

But miracles must be planned, nurtured, and ex 
ecuted with intelligence and hard work. The miracles 
of the second half century will prove progressively more 
difficult, more costly, and more time-consuming. They 
will demand great skills from all of us. 

First we must have a clear definition of the miracles 
to be expected. Our magnificent aircraft manu- 
facturing industry exists only to serve its customers 
not itself. It is the customer's requirements that have 
shaped all that we have done and which have carefully 
controlled our rate of progress. My files have been 
bulging for years with proposals for new aircraft, the 
fruit of engineering imagination and ambition. These 
are the fundamental ingredients of our capacity for 
progress, but such beautiful and arresting sketches are 


not progress until a customer says: ‘“‘I want that.” 


Our product must meet someone’s need as fully and 
exactly as we know how to make it. The more widely 
useful our product, the higher will be our level of 
activity and the lower the cost to the customer. To- 
day and throughout all our 50 years of aviation history, 
the Government is the customer for more than 90 per 
cent of our effort. The grave responsibility of the Air 
Force and Naval aviation for our national defense is 
Their weapons must be 


fashion. The 


awesome in this atomic age. 


the human effort can 


most effective 


OF 
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narrow and restrained requirements of the services 
inevitably create extremely costly aircraft. 

The current criticism of the high cost of combat air- 
craft is nothing new. I heard it in 1918, in 1926, and in 
1940. The aircraft is a vehicle, and its cost increases 
directly with its usefulness. This is true in all forms 
of transportation, being lowest for the barge and 
highest for the airplane, but their usefulness varies in 
direct ratio with their cost. It seems fundamental to 
me that an airplane capable of operation at 600 m.p.h. 
at 50,000 ft. is going to cost considerably more than one 
limited to only 400 m.p.h. at 25,000 ft. An aircraft 
capable of seeking out and firing upon an unseen target 
at night in the stratosphere is certainly a more useful 
and, therefore, more costly weapon than one operable 
only in daytime at moderate altitude and using the 
pilot’s limited vision and perception for its effective- 
ness. 

I have thought carefully for many years about the 
pronounced apathy of the American people for the 
flying machine they invented. While living in Santa 
Ana, our family Doctor, then living in Los Angeles, 
mailed the following postcard to my mother: 

“September 30, 1910 
“For Heaven's Sake, if you have any influence with that weld 
eyed hallucinated young man, call him off before he is killed 
“Have him devote his energies to substantial, feasible, and 
profitable pursuits, leaving Dreaming to the Professional 
Dreamers 


H. H. Sutherland” 


And _ this one of Throughout 


my participation in aviation, which spans its life 


was only many. 
history, I have witnessed a public criticism of the high 
cost of aircraft. There seems to be a pervading public 
mustrust of its aircraft industry and unrelenting pressure 
on the Government to keep vigil on aviation profits 
which resulted, eventually, in severe profit limitations 
by law. Perhaps this concern is a mistaken heritage 
from World War I when the term “‘war profiteer’’ be- 
came almost synonymous with weapons production. 
I do not know why our citizens heartily approve sub- 
stantial profits in the manufacture of automobiles and 
household appliances but become incensed when the 
word ‘“‘profit’’ is even associated with the design and 
production of vital weapons of national defense. 
Without profits there can be no applied research, no 
advanced machinery and equipment, no progress of any 
sort. Without profits for expansion of 
potential, that expansion becomes the responsibility 
of the Government. And that responsibility, un- 
checked, leads inevitably to the political and economic 
system we are attempting to defeat thoughout the 


available 


world. 

There is a proper role for the Government in aviation 
progress. Asa customer, it must establish clearly and 
precisely the goals that we will be required to meet. 
It must assume the financial cost of the research and 
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development work required in the perfection of that 
goal. It must make certain that this nation has, at all 
times, a strong aircraft manufacturing industry, a 
healthy and progressive accessory and equipment in- 
dustry, an experienced and capable air-line industry, 
and a reservoir of educated and skillful young men 
coming forward to man these industries. The burden 
of this role, costly as it may seem, will be borne cheer- 
fully by an enlightened American citizen once he under- 
stands its full significance to him and his progeny. 

When I experienced the indescribable thrill of my 
first flight, I was convinced that here was a new art 
that men of vision everywhere must rush to embrace. 
I was unbounded in my gratitude for having been born 
into a new race of man privileged to witness the birth 
of the air age and to share its triumphs during their life- 
times. For me, this age has fulfilled all of its promises, 
but I am only one of a few hundred thousand of us for 
whom the dream has come true. For more than a 
hundred million Americans it seems the dream has 
proved to be an elusive and a frightening thing. 

The salvation of this great nation lies solely in the 
revival of that dream and its fruits for the American 
people. That nation will be strongest on earth whose 
peoples embrace this new age to their hearts. Its 
destiny must become their destiny, else they perish. 
The achievements and the problems, the victories and 
the costs, the dangers and the promises of aviation—all 
of these things must become a part of the soul of a 
nation if it is to find security with its prosperity. We 
must seek out the roots of public apathy and mistrust 
and, in destroying them, create the omnipotent Air 
Power we have the brains and the tools to do. 

It is not enough to have a core of aviation experts 
with the sympathy and understanding of Air Power 
required for its progress. We must have a whole people 
rising in a great tide of belief in this new art, convinced 
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of its demonstrated ability to usher them into an 


entirely new age of mankind. 


THE FUTURE 


With such support assuredly no miracle is impossible 
for aeronautical industry. With such 
support I have no hesitancy in predicting, for the 


science and 


second half century of flight in America 

(1) Interstellar space ships with speeds of 25,000 
m.p.h. 

(2) Giant 200-passenger jet air liners crossing the 
continents and the oceans in literally zero time in the 
westerly direction. 

(3) Nuclear-energy-powered 
bombers capable of multiple nonstop circumnavigation 
of the globe and second as passenger air liners. 

(4) Helicopters carrying all air-line traffic over dis- 


aircraft, first as 


tances of 150 miles or less. 

(5) Flying boats with speeds matching that of the 
equivalent landplane but with greater capacity for 
cargo, passengers, or weapons. 

(6) Fully air-line flight 
through electronic guidance and control equipment. 

(7) Reductions in the cost of air travel well below 


automatic operations 


that of any other form of transportation. 

These predictions are testimony to my faith in the 
conviction of all of us in aviation today that ‘“‘the 
difficult we can do immediately the impossible takes a 
little longer." The only real problem we face in the 
second half century of flight in America is the same one 
faced by the Wright brothers 50 yearsago: todetermine 
the scientific truths and to follow them with bold 
honesty through all difficulties to inevitable success 
This is our great heritage from the two men we honor 
here, and may God grant us the courage to honor it as 


well. 
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Skin Friction in Slip Flow’ 


S. A. SCHAAF? anpb F. S. SHERMAN? 
TM niversily of California al Berkeley 


ABSTRACT 


Skin-friction data in the slip flow régime are presented, cover 
ing the range: 0.008 < M/V/R. 0.38; 34 < Re 2,020 for 
25 M<38: 3 < Re 500 for Ve~ 0.2 and 0.6 
are discussed in terms of the various relevant theories, and it is 
inferred that, in the neighborhood of Re = 1,000, the dominant 
effect is the interaction layer and external 
flow which increases the skin friction, but that, in the neighbor 
hood of Re = 50, the dominant effect is that of slip that decreases 


the skin friction 


The results 


between boundary 


NOMENCLATURE 


= plate breadth 


drag coefficient = total drag force (both sides of 


plate)/dynamic pressure X area (one side of 
plate 


drag force 


7 


2 
erfe )= e~ 72 dx 
Va Jz 


= fraction of diffusely reflected molecules 
= plate length 
Mach Number (free 


static pressure 


stream 


impact pressure 
stagnation pressur¢ 
cone surface pressure 
) = Prandtl Number 
Reynolds Number = ual vy 
= temperature (absolute 
= free-stream temperature 

= plate thickness 
r = tangential velocity 

= free-stream velocity 

= tangential coordinate 
y = normal coordinate 
15M 


= thermal accommodation coefficient 


= VR 


= specific heat ratio 
= Euler Number, 0.577 
= boundary-layer thickness 
bs = boundary-layer displacement thickness 
\ = molecular mean free path 
Ml = viscosity coefficient (free stream 
= po/p 
= density (free stream 
a = kinetic theory constant, 0.996 
¥ = stream function 


= viscosity-temperature exponent, up» ~ 7“ 
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(1) INTRODUCTION 


GENERAL FIELD OF RAREFIED GAS DYNAMICS 


— 
i.e., the mechanics of a gas so rarefied that mean 
free path effects become important 
Theoretical and experimental re 


has been dis 


cussed by Tsien.' 


sults for several cases of free molecule flow—4i.e., a 
flow in which the mean free path is large compared to a 
typical body dimension——have been presented subse 
quently by Stalder, et al.,° 
Kane 


Experimental results are, in general, in excellent agree 


Heineman,°* Ashley,' 


and Estermann,’ and Sauer, among others 
ment with theory. 

In the slip-flow range, where the mean free path is 
small but not negligible compared with a typical body 
dimension or the boundary-layer thickness, experi 
mental results for sphere drag have been obtained by 
Kane,’ while heat transfer to spheres and cylinders 
has been determined by Drake and Backer" and Stal 
der, et al." 
able body of literature'’ 
slip-flow solutions, mostly based on simplifications of 
the Navier-Stokes equations, and first-order slip and 
Comparisons 


There has also appeared a consider 
dealing with theoretical 


temperature-jump boundary conditions. 
between theory and experiment have been difficult be 
cause the complex flow fields for which the experimental 
results are available are not amenable to easy theo 
the 
which has 


retical analysis. (An exception is the case of 
Couette flow between rotating cylinders, '* 
been investigated experimentally by Kuhlthau." 

It is the purpose of this paper to report the results of 
a series of experiments to determine the drag for a flat 
plate at zero angle of attack in slip flow. It was felt 
that skin friction phenomena, which could thus be iso 
lated, were of considerable intrinsic importance and 
also that there was some possibility of eventual com 
parison between theory and experiment. The experi 
ments cover the range of Mach Number / and Reyn 
olds Number Re (based on plate length and free-stream 
conditions): 3 < Re < 500 for J ~ 0.2 and 0.6; and 
34 < Re < 2,020 for 2.5 < M < 3.8. 
range of the rarefaction parameter .//+/ Re (roughly 
the ratio of mean free path to boundary-layer thickness 
of 0.008 < M/4/Re < 0.38. 


These cover a 


(2) EXPERIMENTAL RESULTS 


The experiments were conducted in the low-density 
wind tunnel at Berkeley.’ The technique of Kane® 


was closely followed and will only be summarized here. 





86 FTOURNAL OF THE AERONAUTICAL SCIENCES FEBRUARY, 1954 





TABLE | 
Summary of Flat Plate Drag Data 
Pressures Static Net 
Drag Plate Dimensions Measured Corrected Mach Pressure, Drag, Drag Reynolds 
i, b t p " F pe No., p D Coefficient, No., 
(in. ) (in.) (in. (micron of mercury M (microns) (mg Cp Re Ps 
1.005 0.5005 0.009 L127 118 1,133 108 3.05 91 350 0.134 904 
933 104 935 96 2.92 82 318 0.149 708 
653 81 650 74 ye 64 267 0.182 162 
455 64 446 59 2.51 52 234 0.230 294 
272 37 259 33 2.55 29 162 0.276 173 
0.812 0.4995 0.009 1,127 118 1,133 108 3.05 9] 297 0.141 734 
933 104 933 96 2.92 82 269 0.155 575 
653 81 650 74 2.73 64 226 0.190 375 
155 64 446 59 2.51 52 187 0.228 239 
272 3¢ 259 33 2.55 29 140 0.294 140 Q) 
0.597 0.5005 0.009 L327 118 1,133 108 3.05 91 266 0.171 540 
933 104 933 96 2.92 82 238 0.187 23 
653 81 650 74 3.73 64 197 0.226 276 
155 64 446 59 2.51 52 158 0.263 176 
272 37 259 33 2.56 29 112 0.321 103 +) 
0.398 0.5010 0.009 L127 118 1,133 108 3.05 91 198 0.191 360 
933 104 933 96 2.92 82 176 0.207 282 
653 81 650 74 2.73 64 144 0.248 184 
55 64 446 59 2.51 52 119 0.297 117 
272 37 259 33 2.55 29 87 0.373 69 
*0.197 0.4999 0.009 1,127 118 1,133 108 3.05 9] 148 0.289 178 
936 104 933 96 2.92 82 132 0.314 139 I 
653 81 650 74 2.73 64 106 0.369 9] 4 
155 64 446 59 2.51 52 81 0.409 58 set 
272 37 259 33 2.55 29 55 0.478 34 2.00 
*(). 4008 0.5010 0.018 1,127 118 1,133 108 3.05 91 229 0.219 362 
933 104 933 96 2.92 82 203 0.237 284 
653 81 650 74 2.473 64 163 0.278 185 
455 64 446 59 2.51 52 131 0.324 118 
272 37 259 33 2.55 29 94 0.401 70 
0.3995 0.5001 0.0045 1,127 118 1,133 108 3.05 91 193 0.186 361 
933 104 933 96 2.92 82 177 0.207 283 1.00 
653 81 650 74 2.73 64 145 0.249 185 : 
455 64 446 59 2.51 52 119 0.296 117 
272 37 259 33 2.55 29 SS 0.377 69 
0.3994 1.1976 0.009 1,127 118 1,133 108 3.05 91 196 0.199 361 
933 104 933 96 2.92 82 448 0.219 283 
653 81 650 74 2.73 64 369 0.264 185 
455 64 446 59 2.51 §2 300 0.311 117 
*0.4000 0.1998 0.009 1,127 118 1,133 108 3.05 9] 89 0.214 362 
933 104 933 96 2.92 82 79 0.231 284 
653 81 650 74 2.73 64 65 0.279 185 0.50 
455 64 446 59 2.51 52 52 0.323 118 
272 37 259 33 2.55 29 4] 0.439 70 
1.005 0.5005 0.009 1,936 143 1,936 130 S.44 103 417 0.092 2,020 
1,339 114 1 ,346 103 3.47 84 357 0.114 ] ,246 
994 87 999 78 3.43 64 305 0.131 911 
541 57 531 51 3.03 43 224 0.182 422 
0.398 0.5010 0.009 1,936 143 1,936 130 3:77 103 266 0.148 800 
1 ,339 114 1,346 1038 3.47 S4 220 0.177 494 
994 87 999 78 3.43 64 184 0.199 361 0.20 
541 57 531 51 3.03 13 128 0.263 167 
396 44 388 38 3.00 a2 105 0.295 120 
(Table 1 continued on page 87) 
Gross and tare drag forces were determined at each of sures by means of the Kopal tables.*? In the sub- 
several tunnel flow rate settings for a series of sting- sonic tests the Mach Number and free-stream static 
mounted flat plates of various length, breadth, and pressure were determined from the measured stagnation 
thickness as indicated in Table 1. The setup for a chamber pressure, the measured nozzle wall tap static T 


aia termination is shown in Fig. 3. The forces ae . ‘ ; 
typical determination 1s shown in Fig. 3. The forces pressure, and the assumption of isentropic flow. Free- expe 

were detected by means of the balance and quartz- : : MI: 
: stream temperatures were determined from the meas- Mac 


spring arrangement used previously by Kane’ and indi- ‘ fr +1¢ 
ted in Fie. 3 sli ured stagnation chamber temperature (standard ther- | ~ 
cated in Fig. 3. , — ' . 
The Mach Number and free-stre: on an mocouple) and the assumption of adiabatic flow. For } the 
1e Mach Number and free-stream static pressure supersonic flow all pressures were measured with the | °'"! 
were calculated for the supersonic case from measured SUPETSON ee Se a ee am 
impact and 5° half-angle cone-static pressures. These Precision U-tube manometer designed by Maslach oun 
readings were corrected for viscous effects™ *! to give OT with the earlier oil manometer used by Kane.° For |} nds 


the ideal values indicated in Table 1 and then were re- subsonic flow the pressures were measured with a pre- | duce 
duced to Mach Numbers and free-stream static pres- cision McLeod gage. cent 
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SEIN FRECTION IN SLIP FLOW 
TABLE 1—Continued 
Pressures Static Net 
Drag Plate Dimensions Measured Corrected Mach Pressure, Drag, Drag Reynolds 
L b t p p pi p No., D Coefficient, No., 
in (in (in (micron of mercury M (microns mg Cp Re 
0.197 0.4999 0.009 1,936 143 1,936 130 Riva 103 205 0.231 396 
1,339 114 1 346 103 3.47 S84 163 0.265 244 
994 87 999 78 3.43 6 136 0.298 78 
541 57 531 51 ,.03 13 RO 0.370 83 
396 $4 388 38 3.00 32 70 0.398 59 
0.3995 0.5001 0.0045 1.936 143 1,936 130 Sa 103 250 0.139 S03 
1,339 114 1 , 346 105 >. 4; S84 208 0.167 195 
994 87 999 78 3.48 64 175 0.189 32 
541 57 531 51 >. 08 43 123 0.252 168 
396 44 388 38 3.00 32 49 0.278 121 
0.3994 1.1976 0.009 1,936 143 1,936 130 3.404 103 648 0.150 S03 
1 339 114 1 346 108 3.47 S4 531 0.178 195 
994 87 999 78 3.43 64 135 0.196 32 
541 57 531 | 3.038 43 303 0.259 168 
396 4 388 38 ; 00 32 244 0.286 121 
0 400 0.1998 0.009 1,936 143 1,936 130 Re wg 103 111 0.155 SO 
1 339 114 | 346 103 3.47 S4 SY 0.179 196 
994 87 999 78 3.43 64 79 0.201 2 
d41 57 531 51 3.03 13 16} 0. 236 168 
396 44 388 ite) , OO 32 ‘bel 0. 267 121 
Static Stagnation Net 
Drag Plate Dimensions Mach Pressure, Pressure, Drag, Drag Reynolds 
L b t No., p p D Coefficient, No., 
in (in in M microns of mercury (mg Cr Re 
2.009 1.0011 0.003 0.566 559.3 695.1 308.5 924 
0.194 1,708 1,753 111.8 0 502 
0.683 313.7 292.0 268.7 0 50) 
0.205 799.1 822.9 82.7 0 251 
0.577 67.15 84.15 131.9 0 64.2 
0.186 930 .7 936 .3 16.6 0 65.6 
0.460 $1] 46 17.93 RO 29 0 10.6 
0.163 121.8 124.1] 34.5 0 30.4 
1.0022 1.0022 0.008 0.546 1,121 1,373 287 .0 0 106} 
0.199 3,188 3, 277 108.5 0 {833 
0.590 541.3 685.1 220.0 0 263 
0.194 1,706 1,751 79.8 0 252 
0.72 202.0 285.6 187.1 125 
0.205 806 .2 $30 .2 59.6 0 126 
0.597 65.33 83.14 89.6 0 32.2 
0.186 232.5 238.2 33.5 0 32.9 
0.466 41.49 18.14 53.9 0 15.3 
0.164 122.8 25.1 25.3 l 15.2 
0.5029 1.0010 0.003 0.196 3,282 3,37 82.0 0 242 
0.586 548.5 692.1 163.9 0 132 
0.196 1,708 1,754 61.1 0 126 
0.707 209.5 292 .4 139.7 0 62.9 
0.207 804.2 828.6 15.8 0 63.0 
0.626 62.39 81.25 66.0 0 16.1 
0.188 231.4 37.2 25.4 l 16.5 
0.516 37.32 44.75 39.7 l 7.69 
0.166 121.6 124.0 18.4 l 7.64 
0.2070 1.0005 0,003 0.586 543.8 686.2 117.4 0 4.2 
0.195 1,714 1,760 $4.4 0.537 52.0 
0.710 205.5 287 .6 97.9 0.745 25.8 
0.207 805.1 829.6 34.3 0.785 25.9 
0.622 62.03 80.52 412.2 1.387 6.64 
0.188 232.8 238.6 18.5 1.776 6.79 
0.521 36.43 43.83 25.9 2.064 3.17 
0.167 122.5 124.9 14.1 3.257 3.15 


The No. 2 and No. 3 supersonic nozzles'* used in the 
experiments produced test-section flows for which the 
Mach Number and static pressures varied by the order of 
+10 per cent in general. For a few of the worst cases, 
the values given in Table | are the result of averaging 
over the space occupied by the model. The departure 
from uniformity in these test-section flows probably 
represents one of the major sources of random experi- 
The subsonic nozzle pro- 
+ (),2 


mental error in these tests. 


duced a flow that was uniform to within per 


cent. 


In the supersonic tests, a systematic attempt was 
made to eliminate aspect ratio and thickness effects. 
Analysis of the data obtained with plates of varying 
thickness and span indicated that such effects become 
less than the general overall order of uncertainty for 
all the unstarred data of Table 1 for supersonic condi- 
tions. Only this material has been used in the figures. 
No systematic attempt to eliminate thickness and as- 
pect ratio effects was made in the subsonic tests. It 
is probable that there are no serious thickness effects. 
The effect of the varying aspect ratio probably con- 
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tributes considerably to the relatively large scatter of 
the subsonic data, but no systematic trends could be 
determined. For the worst case, the approximate 
theory of Howarth*! indicates an aspect ratio effect of 
32 per cent. This theory predicts a decrease in the 
effect with increasing Reynolds Number, a conclusion 
that apparently agrees roughly with the reduction in 
the scatter of the present data for higher Reynolds 
Numbers. 

Summarizing the sources of error, it is estimated that 
the indicated experimental skin-friction coefficients are, 
in general, good to within +5 per cent. It will be ob- 
served that the skin-friction coefficient itself varies by 
approximately 25 to 1 over the range of these experi- 
ments. 

The results are tabulated in Table | and are also pre- 
sented graphically in Figs. 1 and 2, along with various 
theoretical curves and a curve summarizing the experi- 
mental work of Janour® for incompressible flow. 
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(3) THEORETICAL CONSIDERATIONS 


There is as yet no satisfactory theory for predicting 
skin friction in the range of Mach and Reynolds Num- 
bers covered by these experiments. The range of Reyn- 
olds Numbers is too low for the Prandtl boundary- 
layer theory to be completely valid and is at the same 
time too high to warrant use of the Stokes or Oseen ap- 
proximations based on partial or complete neglect of 
inertia forces. On the other hand, the rarefaction 
parameter J//Re is not high enough to warrant use of 
free-molecule flow theory. These three theories con- 
stitute known end points for any applicable theory, 
however, and in a sense bracket the possibilities. The 
corresponding total drag coefficients are: 


"6 (Re > 1), M/Re ~ OJ. 


Boundary-layer theory 


+ 0.365177 X 


w-—1)/2 
= 1)Pr } | 


Cy = (2.6356 V Re) }1 


[(Re x 0), Re ~ O}. 


Oseen theory :*’ 


Cp = 8/Re[1 — y’ — In (Re/16)] 


|(Re ~ 0), 17/Re 


Free-molecule flow theory 


Cy = 1.35/M diffuse reflection) 3 


None of these theories considers the effects of slip, 
which is known to be an important phenomenon in the 
range covered by these experiments. From kinetic 
theory considerations, one has the boundary condi- 


tion at a wall, vy = 0, 


ee 
tain = €2 — 


f\A(Ou/Oy) +... } 


where f is the reflection coeflicient (approximately 
unity for most engineering surfaces) and o is a numerical 
constant very close to 1. There is also a corresponding 


temperature jump boundary condition, 


AT = Tan — Trani = 


») = 
~ a Y 


A OF 
a y + 1 Pr oy 
where a is the thermal accommodation coefficient and 


Pr is the Prandtl Number. 
been calculated by Schamberg 


Higher order terms have 
> and Grad.” 

The effects of slip by itself on incompressible flat- 
plate boundary-layer flow were determined to first 
order by Lin and Schaaf.'’ The first two terms in a 
perturbation series for the stream function were ob 
tained in the form 


YV=WtAMWtr... 6 


There is a decrease in the boundary-layer displacement 
thickness, 


bx /x = (1.73/V Re) [1 — 0.866 (W/V Re)...] (7 
but no corresponding change to relative order \ 6 
(or \/\/Re) in the skin-friction coefficient. This is 
apparently associated with the fact that the tangential 


velocity distribution function for the flat plate boundary | 
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SKIN FRICTION 


layer with no free-stream pressure gradient has an in- 
flection point at the surface. In the general case, a 
correction term of relative order \ 6 to the skin fric 
tion could be expected. Maslen' has extended these 
results to the compressible boundary-layer flow case. 

It should be emphasized at this point that rarefac- 
tion effects are not confined to the phenomena of slip 
and temperature jump, even in the so-called slip-flow 
range. These latter phenomena are associated with 
the boundary conditions only, and in general it is to 
be expected that modifications in the differential equa 
tions used to describe the flow will also have to be 
made. These modifications in the differential equa- 
tions consist 1n the addition of neglected terms of two 
general kinds: (1) terms contained in the complete 
Navier-Stokes equations but neglected in the first ap 
proximation u(O*u Ox") for 
boundary-layer theory or the inertia term p(Du Dt 


e.g., the stress term 


for Stokes theorv—and (2) correction terms to the Na 
vier-Stokes stress and heat flux formulas themselves 
terms! obtained from the kinetic 


e.g., the Burnett 


theory of monatomic gases. (It might be noted at this 
point that it is quite possible that corrections of this 
second kind should be made by considering more equa- 
tions in more dependent variables, rather than adding 
more terms to the Navier-Stokes stress and heat flux 
formulas, following the treatment of Grad.*) 

Writing Eqs. (4) and (5) in dimensionless form, 


Usiix WT O(u u Me O(u u 
~ (S 
uu V Re Oy 4) Re O(y/L) 
Al M ooT/T.) | aes 
aoe = ~ Oo(7 7 (9 
1 V Re Oy 4) e 
one observes that terms of relative order lJ, y/ Re 


neglected in the boundary-layer flow equations must 
be considered if slip at high Re is to be considered, while 
terms of relative order 17, Re neglected by the Oseen 
equations-must be considered for slip flows at low Re. 
Che problem is further complicated by the fact that 
the flow equations depend upon .V/ and Re independ- 
ently and not merely on a rarefaction parameter. 
Corrections to boundary-layer theory in the form of 
asymptotic series expansions of the flow equations for 
the no-slip case and based on the Navier-Stokes ex- 
pressions have been investigated for 1J = 0 by Alden’® 
M > Alden 


finds correction terms to the Blasius expression for the 


and for | by Lees and Probstein.*® 
ocal skin-friction coefficient in the form of an additive 
term proportional to 1 Re Lees and Probstein find 
additive terms for the /oca/ skin-friction coefficient pro- 
portional to 1/* Re for 1J* Re < 1 and proportional 
to M’*/Re’* for M*/~’Re > 1.. (In latter 
range the term dominates the ‘“‘first- 


this 
“corrective” 
order’ term.) The physical basis for the increased 
skin friction is the interaction between the boundary 
laver and the free-stream pressure disturbance induced 
by the The corresponding local 


boundary _ laver. 


ix Shiv? FLO se so 


drag coefficients, it will be noted, are, except for the 
last case, nonintegrable. The analytic solution for the 
flow field must thus apparently be obtained by ‘‘match 
valid near the leading edge to the 


ing’ a solution 


asymptotic solution valid downstream. Some pre 


liminary results of this kind for the case ./ 0 have 
Leen obtained by Lewis and Carrier 

In addition to these terms, which must be retained 
in the flow equations, it is to be expected, as indicated 
above, that higher order approximations to the stress 
tensor and heat flux vector than given by the Navier 
Kinetic 


theory permits calculation of such correction terms for 


Stokes expressions must also be considered. 
the case of a monatomic gas; these are the so-called 
Burnett terms. However, the equivalent expressions 
for air are not available, and the work of Truesdell 

indicates that they may be of considerably greater 
complexity than the monatomic gas relations of Bur 
nett. 

To summarize, it appears that an asymptotic series 
type of solution valid for the range covered by these 
experiments may be possible, correcting boundary 
layer theory with respect to neglected terms in the 
flow equations, including the still-undetermined stress 
and heat flux expressions for air, as well as slip and 
temperature jump in the boundary conditions. It is 
further to be expected, from the partial results of Alden 
and Lees and Probstein, indicated above, that only 
local nonintegrable skin-friction values will result. 
The complete analytic determination for the total drag 
coefficient would probably involve a matching of the 
asymptotic solution to one that is valid near the lead 
ing edge. A direct numerical analysis is probably not 
warranted for air until the higher order stress and 


heat flux terms have been determined. 


(4) APPROXIMATE OVERALL SOLUTIONS 


In view of the great complexity attached to the solu 
tion of the problem of skin friction in the range covered 
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Photograph of experimental setup 
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by these experiments, it has seemed desirable to de- 
velop a rough approximate theory that spans in closed 
form the entire range from continuum to free-molecule 
flow. Such a theory gives an idea of the overall trend 
and assists in the analysis of the data. 

The Rayleigh equation** 


Pott ~(OU/OX) = wao(O°u/Oy") (10) 


is used in conjunction with the first-order slip boundary 
condition, Eq. (4) with f = 1, to obtain the solution 


u 14 ( Vx . *) f ( Ux 
= exp erfe 
1h « 7u x NU. 
Te ; hes" 
ve) = erie ( *) (11) 
4dvx : fyx 


The corresponding total drag coefficient is 


7 2.67 § » , 26 } 
CoM = —  se* erfc (gz) —1+ 
2? | Vr (12) 
3? = Re/2.25M? | 
This has the limiting forms of 
Cp > 8/V rRe, M?/Re—> 0 | 12 
(io) 


M? Re of 


Cp > 2.67/M, 


which are correct in form although not in numerical 
value. Mirels'® has extended this calculation to in- 
clude the effect of compressibility and has introduced 
an empirical relation so as to ensure a correct matching 
for Re > o. In view of the discussion above such 
extensions seem of Both 
the foregoing analysis and that of Mirels imply no 
first-order effect of slip on skin friction, in agreement 
with the more detailed analyses of Lin and Schaaf and 


somewhat limited value. 


Maslen. 

A casual inspection of Fig. | would seem to indicate 
that, within reasonable scatter, the data seem to lie 
on the curve defined by Eq. (12) for high rarefaction, 
the boundary-layer curve of Blasius for low rarefaction, 
and an interpolating curve in between. However, 
closer inspection reveals additional interrelations and 
strongly indicates that no such single correlation of 
CyM vs. M/+/ Re in fact exists. The departure of the 
skin-friction coefficient from its continuum boundary- 
layer theory value seems to be associated with at least 
two phenomena—namely, slip effects and interaction 
effects. (The word interaction is used to include all 
effects stemming from modifications in the flow differ- 
ential equations.) It can be guessed from the work 
of Alden and Lees, et al., that these interaction effects 
will increase the skin friction, and do so increasingly 
with higher Mach Numbers. It appears in Fig. 2 
that in the vicinity of Re ~ 1,000 this phenomenon 
dominates the slip effect, in agreement with the theories 
of Lin and Schaaf and of Maslen which indicate that 
slip effects are initially negligible. However, at lower 
Reynolds Number, Re ~ 50 say, slip apparently begins 
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to be the dominant effect, and the skin friction de- 
creases with increasing Mach Number. Slip can cer- 
tainly be expected to be a dominant effect for 1//+/ Re 
The fact that the Rayleigh limiting value of Eq. 
(13) happens to be approximately as much higher than 


ot, 


the Blasius value as the actual increase at the higher 
Reynolds Numbers in these experiments, due pre- 
sumably to the interaction effects, is probably only a 


fortuitous coincidence. 
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Turbulent Boundary Layers in Adverse 
Pressure Gradients’ 


FRANCIS H. CLAUSER? 
The Johns Hopkins U'nwersity 


SUMMARY 


An experimental study has been made of turbulent boundary 
layers in adverse pressure gradients which were established in 
such a way that for each layer the profiles, when plotted in a set 
of universal coordinates, are of similar shape. A comparison is 
made with the corresponding set of laminar profiles. Results 
ire presented for the effects of Reynolds Number, pressure gra 


dient, and roughness on skin friction. It is shown that boundary 


layers in moderately high adverse gradients will not exhibit what 


is termed downstream stability. A comparison is made of the 


present experimental results with conventional methods of pre- 
dicting the effect of pressure gradients on turbulent boundary 


layers, and it is shown that the agreement is poor 


INTRODUCTION 


O- OF THE MOST CHALLENGING PROBLEMS in fluid 
mechanics today is that of understanding the 
behavior of turbulent boundary layers, particularly 
with regard to their separation or detachment. A 
great deal of our technology is concerned with devices 
that are fluid mechanical in nature such as fans, com- 
pressors, turbines, pumps, torque converters, recipro- 
cating engines, jet engines, aircraft, missiles, ships, and 
the like. In almost all cases, the limiting performance 
of these devices is determined by flow separation. The 
ultimate lift of a wing or thrust of a propeller is deter- 
mined by stalling-i.e., separation. Similarly, the 
maximum compression of an axial-flow fan and the 
maximum torque of a converter are limited by separa- 
tion. Even the simplest water faucet depends on 
separation for its throttling action. 

Near every solid surface immersed in a moving fluid 
there is a layer of fluid retarded by friction. For the 
fluids that are commonly used, together with the ve- 
locities and sizes that most frequently occur, the Rey- 
nolds Numbers are sufficiently large that the layers are 
not only thin but also turbulent. Thus we find that 
the greater number of our fluid mechanical machines 
are directly limited in their capabilities by separation 
of turbulent boundary layers. This separation is 
caused by the action of adverse or rising pressure gra- 
dients, retarding the fluid in the boundary layer, 
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bringing it to rest, and forcing the outer stream to leave 
the surface and ride out over this stagnant fluid. 

In spite of the great technical importance of turbu 
lent separation, our knowledge of the action of pressure 
gradients on turbulent layers is poor. First, no theory 
worthy of the name exists for any turbulent shear 
flow. Second, several empirical methods exist for pre 
dicting the effect of pressure gradients on turbulent 
boundary layers, but experience indicates that none 
of them is reliable. They do little more than correlate 
the set of data on which they were based originally. 

Early, it was found that the state of affairs in the 
turbulent boundary layer depends not only upon the 
local environment such as the local values of the pres 
sure gradient, the wall shear, the surface roughness, the 
boundary layer thickness, etc., but also upon a long por 
The forces acting on 
Skin- 


stresses 


tion of the history of the layer. 
surprisingly 
friction coefficients that 
are only a few thousandths of the kinetic energy per 


a boundary layer are small. 


show shearing 
unit volume. Furthermore, because the layers are 
thin, the pressure forces that act upon an effective 
area across the face of the layer are also very small in 
size. With such diminutive forces present, the bound- 
ary layer reacts slowly to a changing environment 
It is for these reasons that the previous history of the 
layer is so important. 

When we began our turbulent boundary-layer re 
search at Johns Hopkins, we wanted to study in detail 
a few cases that would have simple and well-defined 
pressure histories and which could be measured with a 
In specifying the history, non- 
The 
pressure gradient acts across some effective face area 
6’ of the layer with a force per unit length and width 
of 6’(dp dx). The only other gross force on the layer 
is the wall shear, 7». Consequently, the boundary 
layer will have a constant history, nondimensionally 
speaking, if (6’/7)) (dp/dx) is held constant throughout 
its length. Since a constant history is a particularly 
simple and well-defined case, we decided to attempt to 
obtain such cases, even though we were not sure that it 
would be easy in practice to set up such ‘“‘simple’”’ 
histories. Throughout the rest of this paper we shall, 
for convenience, call this class of profiles ‘‘equilibrium”’ 
profiles. This paper the results of the 
work accomplished so far on layers having such pro- 


fair degree of precision. 
dimensional parameters must, of course, be used. 


describes 


files. 
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NOTATION 

c. = aconstant 

cf = local skin-friction coefficient 

F = power spectrum 

F,, Fe = Reynolds Number factors for effect of roughness 
and pressure gradient on skin friction 

fig = functions used in describing velocity profiles 

G = family parameter of equilibrium profiles 

H = §*/6 

k = height of roughness elements 

m = constant in Falkner-Skan solution 

n = spectral wave number 

p = local static pressure 

dp/dx = static pressure gradient 

q = local free-stream dynamic pressure 

R = Reynolds Number 

u,v,’ = components of velocity 

U = free-stream velocity 

Ux = friction velocity = V ro/p 

4, 2 = position coordinates 

Au, Ay ‘ . s * = 

7 = shifts of semilogarithmic profiles caused by rough- 

ae ness and pressure gradient 

3 = constant in Hartree solution 

6 = distance to outer edge of boundary layer 

6* = displacement thickness 

6 = effective boundary-layer thickness used in 
(6'/ 7) )(dp/dx) 

A = universal turbulent boundary-layer thickness 

0 = displacement thickness 

p = fluid density 

v = kinematic viscosity 

qT = shearing stress at wall 


EXPERIMENTAL PROGRAM 


Description of Boundary-Layer Tunnel 


Our program of turbulent boundary-layer research 
is a long-range one which has as its goal both the deter- 
mination of the behavior of turbulent layers and an ex- 
ploration of their internal structure. A relatively in- 
expensive wind tunnel was built especially for this re- 
search. In establishing conditions for the tunnel de- 
sign, it was considered necessary that the layers them- 
selves and their associated sublayers should be suffi- 
ciently thick that microinstrumentation would not be 
required. This led to a test section 37 ft. long, 4 ft. 
high, and 3 ft. wide with velocities ranging from 6 to 
10 ft. per sec. Measurements are made on the bound- 
ary layer that develops on the floor. The layers range 
from 2 to 20 in. thick. Because separating boundary 
layers were to be one of the subjects of investigation, 
there was great likelihood that adverse gradients would 
induce separation on the walls and ceiling before it oc- 
This indicated the need for bound- 
To ac- 


curred on the floor. 
ary-layer control on the upper three surfaces. 
complish this without suction ducting and at the same 
time to simplify the pressure distribution control, the 
tunnel was arranged as shown in Fig. 1. A large 
blower is used to force air into the tunnel under pres- 
sure. A screened diffusor and a large settling chamber 
with four single-piece precision screens makes the air 
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stream uniform and of low turbulence. The tunnel 
is open-return, and a heavily screened diffusor at the 
exit maintains a controllable overpressure in the work- 
ing section. 

On the 
rounded edges are provided every S in. for bound- 
ary-layer removal, the air simply being blown into the 


Flaps 


sides, simple variable-width slots with 


atmosphere by the overpressure in the tunnel. 
are provided on the top for exhausting relatively large 
amounts of air to slow down the central stream and es- 
tablish a variable adverse pressure gradient. <A_por- 
tion of the floor is a flexible masonite sheet, controlled 
by automobile jacks for more detailed pressure con- 
trol. <A traversing gear for both pitot tubes and hot- 
wire probes can be used at any location in the test sec- 


tion. 


Preliminary Analysis 

At first sight it would appear simple to establish a 
pressure distribution giving a sequence of equilibrium 
turbulent boundary-layer profiles. Falkner and Skan! 
obtained the theoretical solution for the corresponding 
set of laminar profiles. The Falkner and Skan solutions 
show that, if (6/79) (dp/dx) is held constant throughout 
the length of the layer, then the profiles of u/ UL’ versus 
y/6 are all identical, even though the Reynolds Number 
(and, hence, c,;) is changing with x. Further, if x is 
measured from the origin of the boundary layer, then 
(x/q) (dp/dx) throughout the 


length of the layer. 


is likewise constant 
This shows that L’ is proportional 
to x”, where m = —(1,2) (x/q) (dp/dx). It is of in- 
terest that neither (6/g) (dp dx) nor (x/7) (dp dx 
is constant. Since for these laminar layers the profiles 
are all identical throughout the length of the layer, 
the choice of thickness in the expression (6/79) (dp dx 
is arbitrary since all thicknesses will be in constant 
ratio for a given layer. A simple special case of the 
above analysis is the Blasius solution, obtained when 
the pressure is constant. 

When we turn to the case of turbulent layers, the 
situation is not so simple. In contrast to the Blasius 
case, turbulent profiles are not similar, even when the 
pressure is constant. At lower Reynolds Numbers, 
the profiles are roughly one-seventh power profiles, 
and as the Reynolds Numbers increase, the profiles 
are more nearly one-eighth or one-ninth power pro- 
files. The typical difference can be seen in Fig. 2 for 
the profiles of Hama and Klebanoff and Diehl. Fur- 
ther, it is seen that a profile on a rough wall shows even 
greater lack of similarity. 

Normally, one might overlook the deviation for the 
the smooth plate profiles. However, with our long 
test section, such discrepancies were likely to show up, 
and, since we had a chance to make careful measure- 
ments under well-defined conditions, we did not want 
to have effects due to pressure gradients confused with 
such Reynolds Number effects. 
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DETAIL OF FLEXIBLE SECTION 


A- ADJUSTABLE FLAP 
6 - GLASS PANEL 

C- FLOOR SUPPORT 
D- TURNBUCKLE 

E- JACK 











SKETCH OF WIND TUNNEL FOR 


BOUNDARY LAYER TESTS 





Now it has been known for many years that, if we 
measure velocities relative to the free stream and 
abandon Ll’ as a nondimensionalizing velocity, using 
in its place ue = V 7/p, the friction velocity, then all 
constant pressure turbulent profiles are similar. In 
Fig. 3, a collection of profiles has been plotted, and it 
will be seen they all fall together, including the dis- 
similar profiles of Fig. 2. This way of plotting retains 
the typical shape of turbulent profiles. The customary 
zero velocity point is located at the now variable 


height, — ‘V 2 ‘cy. 


It should be emphasized that this plot correlates the 
profiles independently of both Reynolds Number and 
roughness. When 
parable to the thickness of the layer, then the origin 
Moore? has shown that, even for 


the roughness size becomes com- 


of y is undefined. 
very large roughness elements, an origin exists which 
will give agreement with this universal law, as can be 
seen from his data plotted on Fig. 3. It is believed 
that all turbulent profiles for which reliable wall shears 


are known are in agreement with the method of plotting 


used in Fig. 3. In fact, it is seldom, either in turbulence 
research or boundary layer research, that a collection of 
data from varied sources shows such consistency. 

The method of correlating turbulent profiles given 
in the preceding paragraphs offers a possible foundation 
on which to build an analysis of equilibrium profiles 
It appears reasonable that a family of equilibrium pro 
files should be characterizable by a dimensionless 
parameter containing the pressure gradient as one of 
the factors in the parameter. If this is true, then all 
of the constant pressure profiles must correspond to a 
single member of the equilibrium family, since for all 
these profiles the parameter will be zero. This in turn 
necessitates that a correlation scheme be adopted for 
which all the constant pressure profiles will be a single 
member of the family—i.e., they are all identical when 
plotted nondimensionally. This has been achieved 
by the method just described. 

When we consider equilibrium turbulent layers for 
which the pressure is not constant, it is necessary to at 
tempt to define more precisely what we mean by equilib- 


rium layers. What we should like, of course, is to have 
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a layer for which at every cross section the flow is com- 
pletely dynamically similar. By this is meant that at 


each section both the mean velocity patterns and the 


statistical patterns of the fluctuations will be similar. * 


It is almost certain this cannot be achieved, because 


as the layer moves down the plate its Reynolds Num- 
ber increases and this changes the high wave number 
end of the spectrum compared to the low wave number 
end. This is clearly seen in the recent measurements 
of Laufer,'® which show this change in statistical struc- 
ture of turbulent shear flow in a pipe when the Rey- 
nolds Number is altered. 

On the other hand, this same state of affairs must ex- 
ist in constant pressure layers, and in spite of this an 
important property such as mean velocity distribution 
can be correlated on a single plot even when the Rey- 
nolds Number is varying. Since large eddies deter- 
mine the shear distribution, which in turn determines 
the mean velocities, this probably indicates that for 
constant pressure layers the small wave number (large 
eddy size) end of the spectrum is always similar in 
pattern (i.e., similar as a function of the two variables, 
wave number, and distance from the wall). The ap- 
propriate nondimensional variables would be F l’«*6, 
né, and y 6. 

The foregoing leads to the hope that a like state of 
affairs may exist for equilibrium layers with varying 
pressures. If this is true, then it is plausible that 
equilibrium layers should be those layers for which 
the profiles at each cross section are identical when 
(1 — l’)/ux is plotted against y 6. In this report, 
we shall accept this as a tentative criterion that lavers 
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ire indeed equilibrium layers. Eventually, it will be 
necessary to show that similarity also exists for some 
significant set of properties of the fluctuating quantities. 

There also remains the question of obtaining a proper 
of the (6’ ro) (dp dx). We shall 


parameter form 


show presently that the choice of 6’ is not arbitrary as in 


the case for laminar lavers. This point is discussed 


in some detail later 


Experimental Procedures and Results 


It was apparent at the outset of the experimental 
work that it would be necessary to measure wall skin 
friction. Not only is + 
in itself but, as has been shown above, it is also a pri 


a useful piece of information 


mary datum in the correlation of turbulent profiles 
An instrument for measuring 7) 1s being built; however, 
this instrument was not ready for use when the present 
results were being obtained so an alternate method was 
used Ludwieg that 
smooth wall, the mean velocity points fall on the well 


and Tillmann*® found near a 


known universal curve of “ Us VS. Vue v, even when a 
A collection of data has 


It is 


pressure gradient is present. 
been plotted on Fig. 4 to show this agreement. 
sufficiently convincing that the following method of 
measuring shear has been used. Since wu Us 
(u/U yl y) Voc,/2, 
versal curve of Fig. 4 can be replotted as a universal 


5.T 


V2 c, and yux v the uni 


family with c, as a parameter as shown in Fig. 
Experimental points taken near the wall are plotted 
on the figure as shown, and c, is determined by selecting 
the appropriate member of the family which fits the 
points. This method should prove to be of value in 
turbulent boundary-layer research, since it provides a 
method of obtaining the skin-friction coefficient simply 


i Ina private communication, D. Ross has independently pro 


posed to use this same method 
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from the mean velocity profile. This method should 
make available for analysis a great deal of data taken 
in the past for which no skin-friction coefficients were 
obtained but for which accurate velocity profiles were 
taken. 

For the first attempt at establishing an equilibrium 
set of profiles, it was decided to use a relatively mild 
dimensionless gradient. This was done because it was 
known that, with gradients that produce near sepa- 
rating conditions, the velocities close to the wall are small 
and the turbulence levels large, making accurate meas- 
urements difficult because of the lack of existence of 


instruments having strictly linear response. This ini- 
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tial choice was fortunate for other reasons, as will be 
seen presently. The procedure adopted was to set 
up a trial pressure distribution and to measure the 
velocity profiles at a number of stations, determining 
the skin-friction coefficients and making a plot of 
(wu — U)/ux vs. y/6. For the first trial pressure dis- 
tribution, the profiles were not similar, and it was 
necessary to alter the pressure distribution a great num- 
The 


final pressure distribution for this first set of measure- 


ber of times before similar profiles were obtained. 
ments is shown in Fig. 6. The corresponding set of 
similar profiles is shown in Fig. 7. Although the pro- 
files were similar from x = 83 to 887 in., the adverse 
pressure started very near the beginning of the test 
A tripping wire was used at the entrance of 
The 


0 to S83 in. was 


section. 
the test section to induce a uniform transition. 
initial run of the turbulent layer from x 
used to give the boundary layer a chance to settle down 
and attain a reasonably uniform internal dynamical 
The 
down process have not yet been explored. 


state. conditions existing during the settling 
It is ex- 
pected this will be done later in the program. 

It will be noticed that the profiles of Fig. 7 show an 
expected difference from those with constant pressure. 
Compared with this difference the profiles show little 
scatter from a single curve, indicating that they are 
nearly identical members of what must be an equilib- 
rium set. The same profiles have been plotted in the 
conventional Fig. 8. It 
will be seen that there is a consistent variation with in- 
as the Reynolds 


manner of u/U vs. y 6 in 


creasing + caused by the decrease of ¢ 
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TURBULENT BOUNDA 


Number increases. This variation of cy with x is shown 
The change in shape of the profiles on Fig. 8 
Since JZ 


for constant pressure profiles is about 1.3 and for sepa 


in Fig. 9. 
causes a change in // from about 1.4 to 1.5. 


rating profiles is about 2.5, it is clear that this set of 
profiles is far from separation. In spite of this, the 
skin-friction coefficients are much smaller than corre- 
sponding values for constant pressure profiles. 

Next it was decided to attempt using a steeper di 
mensionless gradient. Immediately, serious diflicul- 
ties were encountered. The first of these was the ap 
pearance of strong disagreements with the simple two 


dimensional momentum equation, 


dé C; (T+ 2) 6 dp 


dx Z s 


g dx 


The c, term was insignificant, and the dp dx term was 
larger than the d@/dx term by as much as a factor of 
The 
second difficulty was an inability to obtain a set of con 
stant for the 
Corrections in dp/dx which had previously led to con 


two, indicating a lateral divergence of the flow. 


dimensionless characteristics profiles. 
vergence in obtaining equilibrium profiles now produced 


wild variations at downstream stations. After con 
siderable work these difficulties were understood and 
corrected. A detailed discussion of the reasons will be 
The final profiles obtained for the steeper 
10. The 
pressure distribution itself is shown on Fig. 6, along with 
The skin-friction dis- 
It is felt that these 


two sets of constant history profiles, along with the data 


given later. 
pressure distribution are plotted in Fig. 
the first pressure distribution. 

tribution is also shown on Fig. 9. 


on constant pressure profiles, give a useful cross section 
of such nonseparating turbulent profiles in adverse 


gradients. 
ANALYSIS AND DISCUSSION 


Parameters for Analysis 


A question arises as to what parameters can be used 
to analyze the data that have been obtained. First 
consider the boundary-layer thickness. The distance 
6 to the outer edge of the profile has the difficulty that 
the outer edge of the layer is ill-defined. The displace- 


ment thickness 


(wu, U)| dy 


are not satisfactory for analyzing turbulent profiles 
because they cannot be obtained directly from the uni- 
versal plot that forms the basis for our correlation of 
turbulent profiles. To see this we note that, for each 
member of the family, (4 — UL’), us is a universal func- 
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tion of y,6. Thus if we put (x 


y, 6 = n, then ¢ = /(n) and we have 


. 


9 
6* = —b4/— f(n 
No, J 


Since the integral is simply a number, the ratio of 6* 


dy 


to 6 is not fixed, but is a function of c,—that is, of Rev- 
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nolds Number and roughness. A _ similar statement 
holds for @. 


Number and roughness makes them unsuitable for 


This dependence of 6* and @ on Reynolds 
analyzing a universal family. It is, of course, desir- 
able to have a thickness that does not need these addi 
tional data for its specification. It is clear that the 
simplest universal thickness that is so defined is 


"?U—4 ia 
A= dy = —6 {(n) dn 
J0 Ux J0 


For constant pressure profiles the integral has the value 
5.6. Each set of equilibrium profiles will have a differ 
ent value for the integral, but for each set A is in con- 
stant ratio to 6. This means that we may substitute 


A for 6 in our universal plots and eliminate the uncer- 
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tainty in determining 6. It will be found that A is well 
suited to analyses of turbulent profiles. 


In addition to A, we need a suitable parameter to 


characterize the various members of our universal 

family. Customarily, // = 6* 6 is used; however, this 

suffers from the same shortcoming as 6* and 6. When 8 
* 4 


profiles are plotted using u UL’, 6* is a first moment of 


the profile around the axis ul’ = 1. Similarly, @ is a 
mixed second moment around the axes u Ll’ = O and 
ul = 1, and // is the ratio of these two moments 
As a substitute for //, we can ask for a ratio of mo 
ments on our universal plot. Since the only axis with 
significance 1s (1 l’) us = 0, the proper choice of 
parameter is the ratio of the second to the first moment 


around this axis—-1.e., 





U-—wa 
/ ( yay A 
0 {x 
We can relate our new thickness, A, and integral family 
and // as follows 


parameter, G, to 6, 6*, 


&& = Vc,/2 A 
6= Voe/2(1 —GVc,/2)A | 
TI = 1 (1 - GN €e/2 


In addition to the integral parameter G, we also have 
the dynamic parameter (6’ 79) (dp, dx), which, if 6’ 1s 
properly chosen, will be constant throughout the length 
course, be 


of an equilibrium layer. This should, of 
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true independently of the value of cy However, we 
have seen that the various thicknesses 6, A, 6, and 6* are 
not all in constant ratio independent of cy, Conse- 
quently, the choice of thickness to use for 6’ is not 
arbitrary, and we do not know ahead of time which of 
these thicknesses, if any, is the one that will be ap- 
propriate for equilibrium layers. As yet there appears 
to be no sound theoretical basis for determining which 
thickness is the proper one. Experimentally this can 
be obtained fairly simply by establishing equilibrium 
profiles on both smooth and rough plates, taking care 
that both have the same values of G, and then finding 
what choice of 6’ will make (6’/7)) (dp/dx) constant. 
This has not been done as yet. 


Analysis of Properties of the Family of Equilibrium Profiles 


The data from the two sets of equilibrium profiles, 
plus the accumulated data for constant pressure pro- 
files, give us three points in plotting the characteristics 
of the family of equilibrium profiles. A summary of 
the values found 1s given in the following table. 


Designation G A/é 
Constant Pressure 6.1 3.6 
First Pressure Distribution 10.1 6.4 
Second Pressure Distribution 19.3 12.0 


In Fig. 11 the profiles have been replotted in terms 
of the universal parameters given above. The ratio of 
A to 6 has been plotted against G in Fig. 12. We also 
have sufficient information to show the variation of 
5*, 0, and // with c;. This has been done in Figs. 13, 


14, and 15. For convenience, a scale is also given 


AL 
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showing the variations of Reynolds Numbers on a 
smooth flat plate corresponding to the values of c, 
It is apparent that variations of 6*, 6, and // with c, are 
not negligible, showing that an analysis using them as 
universal parameters would be unsound. 

At this point a word of caution is in order. Thus 
far we have assumed that the entire profile is repre 
sented by the universal function of (#'— UL’)/ux vs 
y/6. This is not true near the wall. Ona smooth wall, 
the profile departs from this universal law in the vi 
cinity of the laminar sublayer. On a rough wall, the 
departure occurs at a distance from the wall compar- 
able to the height of the roughness. The integrals 
for A and G for a given distribution of (1 U) /tx vs. 
y,/6in the outer portion of the layer will have the values 
quoted above only as long as the laminar sublayer or 
the roughness elements are small compared with the 
Fortunately this is usually 
On 


thickness of the layer. 


true. This same effect appears in another way. 
the profiles of Fig. 11, zero velocity occurs at an ordi- 
nate of — V 2/c,. 


greater values of y/6 on the profile. 


As cy increases, this point occurs at 
If this value of 
y 6 is not small, then substantial portions of the uni- 
versal profiles will correspond to negative values of 
u, Whereas actual values of u will be positive correspond- 
ing to the conditions in the laminar sublayer or in the 
roughness elements. This will give false values for 
éand //. On Figs. 14 and 15, the solid curves have 
been terminated where this error is no longer negligible 
The dotted portions show this incorrect prediction. 


Family of Laminar Boundary-Layer Profiles 


Since it is of interest to compare our results for the 
universal family of turbulent boundary layers in ad- 
verse pressure gradients with the corresponding results 
for laminar profiles, we pause at this point to review 
We start by 
asking if a pressure distribution exists such that all 
the profiles are similar—i.e., that wu, U7 = f’(y 
6 is some characteristic thickness of the profile and / 


the properties of the laminar family. 
6), where 
For convenience of integra- 


and 6 are functions of x. 


tion, we start with /’ as the derivative of a function 
(stream function) which can be taken so that /(0) = 0. 
We put 7 = y,6 and insert the expression for u in the 


continuity equation, 
(Ou/Ox) + (Ov/Ov) = 0 


Remembering that U’ and 6 are functions of x only, we 


find 


where prime denotes differentiation. 
Similarity demands that 6U’’/6’U be a constant. We 


can integrate thus 


= 5’ nf’(n) —_ (: os : 
0 
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the result into the 


equation of motion, giving 


substitute boundary-layer 


and 


This equation was studied by Hartree,’ who chose the 


as vet undefined thickness, 6, in such a way that 
(o/U'S? = 1 + (FP U/6U" 


Since 6’'U’/6L” is constant, U6" v must be a constant, 


which Hartree calls @. Falkner and Skan use the 
constant m = 8 (2 — ~g). The equation for f be 
comes 

72-1) = fee 


The two differential equations for L’ and 6 as functions 
of x can be integrated explicitly, and, when fitted to 
*the appropriate initial conditions, the solutions are 


u = Cx” 


2x 2p 


= \ = \ ns 

(m + 1)L (m + 1)Cx 

Hartree has obtained a family of solutions of the 
equation for f. Using them, we can compute 


> . 


s* = J = [1 = (Ss U’)] dy ce sf (] ig dn 
J 0 0 


> . 


6 | u/U[1 — (u/U)] dy = | me f’) dn 
0 0 


We also see that 


6 dp = B J "(1 = f’)dn 
tT) ax f'(0)Jo 
x dU x dp 
= — = Mm 
U dx pU? dx 


We 


have obtained numerical values for 6 6, 6* 6, Z/7, and 


both of which are constant independently of x. 


(6 7) (dp/dx) for the Hartree profiles. These pro 
files are plotted in Fig. 16 with an accompanying leg 


end of useful information. 


Effect of Pressure Gradient on Skin Friction for Laminar 
and Turbulent Profiles 
Most methods of treating turbulent boundary layers 
in adverse pressure gradients assume that the pressure 
gradients have no effect on the skin friction. Lud- 
Klebanoff* 
showed that adverse gradients reduced c,, but neither 


wieg and Tillmann* and Schubauer and 


group was able to make a quantitative statement about 
the relationship of c,and dp/ dx. 
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Early in the course of the present work it was ap 
parent that pressure gradients had an unexpectedly 
large effect on skin friction. Even though our first pres 
sure distribution had relatively gentle gradients, the 
skin-friction coefficients were nearly half of those for 
constant pressure profiles at the same Reynolds Num 
ber. 

By combining the results of our present investigation 
with those of others in the field, we are able to present 
a relatively complete picture of turbulent skin friction. 
As mentioned previously, Fig. 4 shows that for smooth 
plates there is a region outside the laminar sublayer in 
which the velocity distribution can be represented by 
the law 


u/ux = 5.6 logy (usy/v) + 4.9 


even when pressure gradients are present. Prandtl’ has 
shown that, when the wall is rough, the semilogarithmic 
part of the curve is displaced downward parallel to 
itself by an amount Az, ux, which depends on the 


roughness Reynolds Number kuy,v. In Fig. 17 are 
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shown experimental data on the relations between these 
two parameters for the case of sand grain roughness 
and also for transverse square bars and wire screening 
laid on a flat surface. It will be noticed that, when 
uxk/v is greater than approximately 100, the asymp- 


totic form is 
Au ux = 5.6 log (t%k/v) + constant 


Since yu», v at the outer edge of the laminar sublayer 
is approximately 50, we see that, when the roughness 
is greater than about twice the thickness of the laminar 
sublayer, the universal law for the turbulent region near 


the wall becomes 
u ux = 5.6 log (y/k) + constant 


The above results are known to hold for smooth, mod 
erately rough, and extremely rough walls without pres- 
sure gradient and for smooth walls with pressure gra- 
dients. Dimensional reasoning, which shows that the 
effects of pressure gradients on conditions near the wall 
must be small, indicates that the results are also true 
for rough walls with pressure gradients. 

For all these cases, there exists a turbulent region 
next to the wall for which 1, is linearly related to log 
y. Turning now to our universal plot of («# — LU), ux 
vs. y, A, we readily see that, for the inner portion of this 

l’),/u* should be linearly proportional to 
Further, the slope of this semilogarithmic 


mot, (4 — 
log y/A. 
curve should be 5.6 1.e., the same as that on the plot 
of u/usx vs. log yusx/v or log y/k. This is seen to be 
true in Fig. 1S, where the members of the equilibrium 
family have been plotted on semilogarithmic scales. 
The existence of this region of overlapping between the 
two different plots on which data for the inner and 
outer portions of the layer may be correlated is of great 
importance. The overlapping can be used to prove ab 
initio that this portion of the layer must be semiloga- 
rithmic,’ and it serves to establish skin-friction laws for 
turbulent layers. 

On one plot, “ “ is a universal function of vx, v 
independent of the free-stream velocity or boundary- 
layer thickness. The function depends on the surface 
roughness Reynolds Number wk, v. 
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On the other plot, («@ — lL’) us is a universal function 
of y \ with the function depending on the pressure 
gradient but being independent of the viscosity (Rey- 


nolds Number) or of the roughness. We put 


U Ux = 2(VUx Vv 


(4 — U)/usx = f(v/A) 


These two can be rewritten 


u/ux = gl(v/ A): (ted v) 
u/ux = {(y/A) + (U/us) 


Now the experimental data show that there is a com 
mon region near the wall where these two quite different 
relations will simultaneously correlate the velocity 
In the overlapping region, both represent 
How- 


profiles. 
the same functional relation between wu and y. 
ever, one contains the factor u%A,/v as a multiplier of 
y, A inside the function, while the other contains Ll’ ux 
as an additive term outside the function. The equiva 
lence of the two equations is maintained even though 
uxA v and l’ ux vary through a wide range of values. 
This implies both that a relation exists between iA y 
and l’ ux and that the functions f and g must be loga 
rithms, since this is the only function for which a multi 
plicative factor inside the function has the same effect 
as an additive factor outside the function. This, of 
course, is in agreement with the experimental data 
that show that the functions are indeed logarithmic. 


The function g can be written 


u/ux = 5.6 logw vux/v — (Any /te) + 4.9 
and the function f can be written 
(a — U)/ux = 5.6 log (y/A) — (Aue/ux) + 0.6 


for those portions of the curves on Fig. IS that are 
semilogarithmic. The function Am, ux is the shift in 
vertical ordinate of the straight line on Fig. 4 when 
roughness is present. It is a function of the variable 
kux/ vas shown on Fig. 17. The function Ams» is the 
corresponding shift of the straight-line portion of the 
The shift is measured from the line 


It is a function of the 


curves on Fig. 1S. 
of the constant pressure profile. 
pressure gradient, G, or equivalently of the parameter 
(6’/t9) (dp/dx). 

The relationship between l’/ux and ux, vy, men 
tioned above, can be obtained by subtracting the above 
equations, which eliminates both uw and y; thus, 


= (Av; Ux ) + 
(Atto/Ux) + 4.5 


U' ‘ux = 5.6 log (tx%A/v) 


or since ts/U = Vc,/2 and Vc,/2A = 6", 


\ . = 5.6 log R;, — Am (Rial) . S 


Ux 
Air (G) 4 
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sills Phis equation is of remarkable significance. It is the 9!@— a 
aiiniies universal skin-friction law for all turbulent boundary “ou ’ 
eas. lavers having equilibrium pressure gradients. It in- 
; cludes the effects of Reynolds Number, pressure gra- aos! \ 
dient, and surface roughness. It is probably the most — ™ 2 £56 log Rye + 43 
general ‘‘universal’’ law that exists, because other pres- ggg, ~\ f 
sure distributions will have an arbitrariness of defint- ba | 
tion that will make the appelation ‘‘universal’” inap- 004} > 
propriate. ee 
Chis law contains the four parameters c,, R;,, G, we oe 
and R,. Rs, may be replaced by R, and // by G from i wal a | 4 oo oe | I 
i . at P) 4 5S 6 
the following formulas 10 ” Oo Re «2 od - 
— FIG.20 CONVENTIONAL TURBULENT SKIN FRICTION LAW 
erent R,=(1-—G Vec,/2)R; 
re ox FOR SMOOTH PLATES WITH CONSTANT PRESSURE 
wioiaa HW = 1/(1 — GV, 2 
low- Che existence of four parameters makes the skin-fric 
er of | tion law difficult to plot for practical use. In Fig. 19, 
Uy a useful set of curves is given for the case of smooth 
liva plates where k = 0. Lines are shown for both G con- ey Oa - 7 Sa 
pugh stant and // constant. Although // is frequently more Bis 
lues. convenient to use, a given equilibrium boundary layer | 
piv will move along a G curve in its development. Al + 
ES- though it is cumbersome to plot the four-parameter 100} 
ulti- |  skin-friction law, it can fortunately be put in a form 
ffect for ready engineering applications. If we introduce e [ 
3, of two factors /; and fs, the skin-friction law can be ae 
data written | 
1o;-— 
V 2/C, = 5.6 log [R;, Fo, Fi] + 4.3 - 
| 
When F; and fF) are unity, this is the well-known skin- a 
friction law for smooth plates with constant pressure, | 4 l | 
which has been plotted in Fig. 20. The factors <a 10 to, 1000 
a = v 
—E—E——O eee FIG.2| REYNOLDS’ NUMBER FACTOR FOR THE EFFECT 
are have been plotted in Figs. 21 and 22, using the data of OF SAND PAPER ROUGHNESS ON SKIN FRICTION 
t in the present report for pressure gradients and the 
hen Prandtl-Schlichting data for sand grain roughness. 
ible In use, after Ff; and Fy have been obtained from Figs 
the 21 and 22, a fictitious Reynolds Number R, /2, F; is 
the computed. This is used in place of the actual Reynolds — jgg9g——______—— er 
line Number, R;,, to find c,in Fig. 20. 
the ° : eo , s [— 
Referring back to the analysis for laminar profiles, 
iene the formula 
| , or 
en- | 4 
| . ig Bars 7 
x 
; ' \ peste 
‘ —<= 
04 + WY 
: WP leat Sain i5 20 
15 FiG.19 EFFECT OF PRESSURE GRADIENT "a AMINAR ANDO TURBULENT Skit FIG 22 REYNOLDS’ NUMBER FACTOR FOR THE EFFECT 
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FIG.23 EFFECT OF H ON LAMINAR SKIN FRICTION COEFFICIENTS 


2f"(0) f{* . K(/1) 
Cy = (1 — f’) dn = 
R,, 0 R;, 


was obtained. The factor A(//) has been evaluated 
and is plotted in Fig. 23. Laminar skin-friction co 
efficients have been plotted in Fig. 19, along with the 
data for turbulent layers. It will be noticed that the 
effects of pressure gradient on skin friction are large, 
both for laminar and turbulent layers. Gradients 
considerably smaller than those required to produce 
separation are effective in reducing the surface friction. 
It is felt that this reduction is of great technical impor- 


tance. 


Comparison of Present and Past Results 


It is of interest to give a comparison of the results of 
the present paper with those of previous authors. 
Since an understanding of turbulent processes under 
the combined influences of shear and pressure gra- 
dients does not exist, analyses have necessarily used 
the momentum integral equation, 

dé ens (HT + 2) 6dp 
dx 7 2 gdx 


as a Starting point. Since the momentum equation is 
not complete in itself, it is necessary to make assump- 
tions about how the terms in the equation are to be 
obtained. These assumptions may generally be cata- 
loged as follows: 

(1) It is assumed that only a single family of pos- 


sible turbulent profiles exists. This permits the neces- 
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sary relationship between 6* and 6 to be determined. 
Different authors have used different parameters to 
characterize the various members of the family of pro- 
files. 

(2) It is necessary to introduce an auxiliary equation 
to determine the successive values of the family param- 
eter. This auxiliary equation is invariably based on 
empirical correlations of collections of experimental 
data. 

(3) A law for the variation of skin friction must be 
assumed. Most authors have assumed a power-law 
variation with Reynolds Number, neglecting any influ- 
ence of pressure gradient on skin friction. 

In making the comparison we shall follow the above 
three points. We shall deal principally with the method 
of analysis proposed by von Doenhoff and Tetervin,’ 
since this method is one of the better known in this 
country and is generally typical of the rest. 

Von Doenhoff uses a universal family of profiles in 
which w,/ LU’ is plotted against y/@. The family param- 
eter is //. Since all the quantities entering into this 
formulation may be obtained directly from a velocity 
profile, it is an easy matter to check to see if any given 
profile fits the family. This has been done for typical 
profiles from our first and second pressure distributions, 
and the results have been plotted on Fig. 24. The 
great and consistent difference between the two sets 
of profiles, which should agree, was quite unexpected. 
It had been felt that this single-parameter family as 
sumption was one of the relatively reliable points of 
prior analyses. The difference shown is sufficient to 
cast serious doubt upon this assumption for general 
boundary-layer behavior. The difference is probably 
caused by the fact that in our profiles /7 was nearly 
constant whereas, in von Doenhoff’s profiles, 77 was in- 
creasing fairly rapidly. Apparently this has an effect 
on the shape of the profile. 


Von Doenhoff uses, as an auxiliary equation, 


ot 68(H —2.975) 2 6 dp 


To dx 


dil 
iY 


2.035 (T — 1286) 
dx 


In analyzing his data, he assumes that the base profile 
for zero pressure gradient is a one-seventh power pro- 
file for which /7 = 1.286. As explained earlier, this is 
not in agreement with our knowledge of constant pres- 
sure profiles as shown in Figs. 2 and 3. Actually, 
IT varies with c, even for constant pressure, as shown 
in Fig. 15. If the range of Reynolds Number is small 
and only smooth plates are considered, the variation of 
IT with c;is not large, and it can be argued that neglect- 
ing the variation is not serious. Consistently with this, 
one may follow von Doenhoff in viewing any /7 = con 
stant set of profiles as analogous to what we have been 
calling equilibrium profiles. The reasoning here is 
analogous to that given later in connection with Fig. 
26. For // = constant, von Doenhoff’s equation be- 


comes 





lined, 
‘TS to 
f pro- 


ation 
iram- 
‘d on 
ental 


st be 
r-law 
influ- 


bove 
thod 
vin,® 

this 


S in 
‘am- 
this 
city 
iven 
ical 
ons, 
The 
sets 
ted. 
as- 
5 of 
to 
eral 
bly 
irly 
in- 


rect 


5 


file 
TO- 
31S 


es- 


wn 
all 
of 
ct- 
is, 
n- 
en 


ig. 
e- 








TURBULENT 


(6/ ro) (dp/dx) — 1.018 (HT — 1.286) 


This can be compared directly with our data, which 
have been plotted as a band on Fig. 25 along with von 
The results are quite discrepant. 
i.e., for 


Doenhoff’s line. 
The only agreement is for (6/70) (dp/dx) = 0 
constant pressure profiles where they must agree rea- 
sonably well. After that, the trends are quite different. 
Von Doenhoff's equation predicts a small value of 
(9 7) (dp/dx) at separation, whereas in actuality it 
must be infinite because 7) goes to zero. His linear 
relation misses completely the bending over of the 
experimental data. As will be explained in more de- 
tail in connection with Fig. 26, the reason the present 
data are plotted as a band is because the parameters 
H and (6/7) (dp dx) are not consistent with the uni- 
versal parameters we find appropriate for turbulent 
layers. The scatter is caused by the variation of /7 and 
@ with cy. 

For skin friction, von Doenhoff uses a power-law 
variation with Reynolds Number and assumes that 
there is no effect of pressure gradient. As has been 
seen, the present results show that the pressure gra- 
dient effect is large. 

In summary, it is seen that the experimental results 
we have obtained show little or no agreement with the 
predictions of the von Doenhoff method. The same is 
true of other methods with which the author is familiar. 
Since these methods do not agree with the relatively 
sumple and well-defined experimental results we have 
obtained, one can conclude that the field is still wide 
open for the advent of a reliable method of predicting 
the behavior of turbulent layers under the influence of 


pressure gradients. 


Downstream Stability 


It will be recalled that difficulty was experienced in 
obtaining constant values of G for the steeper set of 
pressure gradients. Throughout most of turbulent 
boundary-layer work there exists a fundamental propo- 
sition that is seldom enunciated explicitly but which 
is generally believed to be true. It states that a turbu- 
lent boundary layer has a downstream stability such 
that in a given environment the effect of any small 
perturbation, such as local changes in roughness or 
pressure distribution, will die out downstream, leaving 
no trace other than a possible slight thickening of the 
layer. Thus, the effect of a crack or a wave in the wall 
will have no significant effects on the structure of the 
boundary layer far downstream. Klebanoff and Diehl? 
found, indeed, that constant pressure layers can be 
violently disturbed and yet return to a universal state 
downstream, several times thicker but otherwise none 
the worse for their rough treatment. Similarly, for our 
experiments with the flatter pressure distribution, the 
layer possessed this downstream stability, permitting 
corrections to be made at one point without unduly dis- 


turbing downstream points. 
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FIG.25 COMPARISON OF PRESENT AND PREVIOUS RESULTS FOR 


H versus © $2 FOR EQUILIBRIUM TURBULENT LAYERS 


When work was first started on the steeper gradients, 
it was not realized that the situation might be different. 
Corrections were made simultaneously at several points, 
and it quickly became apparent that something was 
wrong. It will be seen that downstream stability is 
necessary to ensure that such simultaneous changes will 
lead to a converging result. Since we were sufficiently 
far from separation that upstream influences were not 
present, the chaotic response that occurred when this 
procedure was followed indicated that some kind of 
downstream instability might be present. The follow 
ing analysis shows why this was so. 

Intuitively, one believes that (6/79) (dp/dx) should be 
a reasonably well-behaved dynamic parameter. For 
constant pressure profiles it is zero. For separating 
profiles it is ©, because 7» vanishes at separation. In 
between, it should provide a reasonably reliable indi- 
cation of the state of the layer. Particularly for 
equilibrium profiles, it should monotonically increase 
as one proceeds from constant pressure profiles to 
separating profiles. As explained previously, its value 
for turbulent profiles will be somewhat dependent on 
cy, but this will not affect the present analysis. The 
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difficulty in its use is that when one is setting up a pres- 
sure distribution, either for a boundary-layer test or 
for the more general problem of the design of fluid me- 
chanical machinery, one deals not with 6 and 7) but with 
x and p or g. Correspondingly, one would like to use, 
instead of (@ 7) (dp dx), dimensionless parameters 
such as (x,g) (dp dx) or (x Ul’) (dU dx). For specific 
families such as the equilibrium profiles, we can ask for 
the relationship between the two types of parameters. 

To illustrate, consider the equilibrium laminar pro- 
files. We have seen previously that for these profiles 
(x gq) (dp,dx) and (@ 7) (dp dx) are both constant 
throughout the length of the layer. Consequently, for 
equilibrium layers there will be an exact relation be- 
tween these parameters. This is shown in Fig. 26. 
Near zero, (vg) (dp dx) increases proportionately to 
(0. ry) (dp/dx); however, when (@ 7) (dp dx) is greater 
than 1, (vg) (dp, dx) is nearly constant at its asymp- 
totic value of 0.1S1. 
attempt to obtain profiles having a value of (@ 7) 


This means that, if one were to 


(dp dx) greater than approximately unity, the pres- 
sure gradient would have to be set with an almost im- 
possible precision to avoid getting a different (4 7) 
(dp dx). Stated differently, all of the members of the 
equilibrium family of laminar layers with rising gra- 
dients are obtained with practically the same pressure 
distribution, p vs. x, except those for which the gra 
dients are relatively mild. Let us see what the effect 
of this condition is on downstream stability. Suppose 
that an equilibrium laminar layer is disturbed so that 
in one region @, say, is increased, (@ 7) (dp dx) will be 
increased, but (x« g) (dp dx) will be unchanged (we 
suppose the boundary layer has no effect on the outer 
flow). Now if the boundary layer originally corre- 
sponded to a point near the origin of the curve of Fig. 
26, the disturbance will shift the point to the right of 
Any local irregularities in the profile will 
In addition, there will be 


the curve. 
be ironed out by viscosity. 
a progressive change in the gross characteristics of the 
profiles because they no longer lie on the equilibrium 
curve of Fig. 26. Intuition and experience say that 
the profiles will asymptotically return to a point on 
this curve. 

Now suppose the original layer corresponded to a 
point on the flat part of the curve of Fig. 27. Again 
the disturbed point moves to the right, but there will 
be no tendency to return because it is still on the 
equilibrium curve. Such a layer would not possess 
downstream stability. 

Let us see what happens in the case of turbulent pro- 
files. As shown earlier, our family of turbulent pro- 
files is not expressible in terms of 6, which is related to 
the universal length parameter A through the skin- 
friction coefficient c, However, since the Reynolds 
Number effect is not great for turbulent profiles, we 
may establish upper and lower limits for the relation be- 
tween (xg) (dp dx) and (6 7») (dp, dx) if we are willing 


to restrict the range of Reynolds Numbers. We have 
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vdp (6 7) (dp dx 
gq dx 0 NC, 
The momentum equation, 
dé Cr (I + 2)@6dp 6 dp 
= 4 -=¢e11+(H +2 P| 
dx 2 2 gdx ; re dx | 
with appropriately averaged c, (denoted ¢é,) can be 


integrated to give 


ia] ey 
ores 
) 


4 dp 
Ty) AX 


ve, 
Substituting this in the equation above, we have 


2(0 ro) (dp dx 


xdp | 


g dx 1+ (7+ 2 


(0 7) (dp dx 


This equation with upper and lower limits on the vari 
ation of c, (and on //, which also varies with C;) has 
been plotted on Fig. 26, using values obtained from our 
number of 


experimental profiles. This band has a 


interesting features. It shows that turbulent boundary 
layers can withstand about three times the nondimen- 
sional pressure gradient (v g) (dp dx) that laminar 
profiles can. The greatest value for this parameter for 
equilibrium profiles is about 0.58. This means that a 
flow with 
b dp 2vn dU . 
= = 0.0S 
q dx U dx 


or l° = const. v~*-** can be decelerated to zero velocity 
without separation, but this is the maximum rate of 
deceleration which can be accomplished without even- 
tual separation. 

In Fig. 26, it is also seen that the turbulent boundary 
layer behaves similarly to the laminar boundary laver 
dp dy 
for small pressure gradients, but the curve quickly 
levels off for values of (@ 79) (dp dx) between 
Above this point, a surprising thing happens. 


in that (xq) (dp dx) is proportional to (@ 7 


and 5. 
The di- 
actually 
decreases for profiles closer to separation. The peak 
values of (xq) (dp dx) lie between 0.5 and 0.58. As 


mensionless pressure gradient (xq) (dp dx 


corresponding to sepa- 
2). At 


2.6." 


(07%) (dp dx) approaches ~, 
(dp dx) 
) +) 


separation, however, // lies between 2.2 and 


ration, (x q) approaches 2 (// + 
The corresponding values of (xg) (dp dx) at separa- 
tion are 0.43 and 0.48, showing a definite decrease 
from the peak values. At first sight this behavior 
seems to be entirely incorrect, for it implies that the 
pressure distributions producing near-separating pro- 
files are definitely flatter than those producing profiles 
well removed from separation. Further, since the 
curve has a maximum, this says that there are cases 
where two quite different boundary layers will be pro- 
duced by exactly the same distribution. 
Both will be shape-preserving, but one will have a 


shape only slightly different from the familiar constant 


pressure 
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pressure profile while the other will be very near sepa- 
ration. A typical example is (xq) (dp/dx) = O45, 
which can be integrated to give 1° = const x~°-*! 
In such a velocity distribution, Fig. 26 shows that an 
equilibnum boundary layer can exist for which (6 7 


dp dx) = 1.2 with /J = 1.3, a value very near that 


for constant pressure profiles. In exactly this same ve 
locity distribution, a second layer can exist for which 
§ r) (dp/dx) = 15.0 and I] = 2.1, which is near the 
separation limit of 2.2 < // < 2.6. This shows that 
two layers in the same pressure distribution can be at 
opposite extremes of the family for adverse pressure 
gradients. 

This unexpected behavior can perhaps best be ex- 


plained by examining the identity 


x dp (6 To) (dp dx) 


q dx 6 XC, 


When the pressure gradient is small, the rate of growth 
of @ with x is governed primarily by c,, and if x is 
measured from the apparent origin of the boundary 
laver, 6 xc, will be nearly constant. For this condi- 
tion, the formula shows that (xv, q) (dp dx) will be nearly 
proportional to (@ 7») (dp dx), as one would intuitively 


think it should. However, as the pressure gradient 


increases, the pressure gradient quickly begins to 
dominate the small skin-friction term in the momentum 
Under this condition, @ xc, 1s more nearly 


(JZ + 2)/2] (0/70) (dp/dx), and 


becomes nearly equal to 2 (// + 2 


equation. 
proportional to 
(v q) (dp dx 
If JJ were constant, (x g) (dp dx) would be constant. 
However, // increases with (8 79) (dp dx), causing (x 
to decrease for large (@ 1)) (dp dx), explain 
The above dis 


q) (dp dx 
ing the observed shape of the curve. 
cussion shows that, for values of (@ 7) (dp dx) greater 
than | to 1.5, both laminar and turbulent layers will 
have little downstream stability and turbulent layers 
can show instability. 

Our experience in obtaining our second set of equilib 
It will be 
remembered that these profiles were closer to sepa 
ration. The (6 7%) (dp dx) 
Whenever a change in conditions was made, the effect 


rium profiles confirms the above analysis. 
was approximately 5. 


of the change was felt with increasing amplitude down- 
stream. However, the rate of divergence was not un- 
managably large. Once the phenomenon was under- 
stood, the method of obtaining equilibrium profiles 
was changed from a process of simply correcting the 
pressure distribution to a process of ‘‘steering’’ the 
boundary layer along a G = constant course. Experi- 
mentally, the experience greatly resembled the process 
of driving a mildly unstable automobile that always 
wants to wander off the road and must be brought back. 
Once a satisfactory equilibrium layer had been ob- 
tained, it showed little tendency to wander off with 
time. The downstream instability is an instability in 


vand notin /. 
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Difficulties Obtaining Two-Dimensional Flow 


During the course of our investigation, considerable 
difficulty was experienced obtaining two-dimensional 
flow. It is felt that these results have considerable in 
fluence on the question raised recently by several au 
thors about the proper form for the momentum equa 
tion. 

Wieghardt" and Mangler'' conducted experiments 
on turbulent boundary layers in rising pressure gra 
dients. They used the two-dimensional momentum 
equation to deduce the values of ¢c,. Their results indi 
cate that c,, instead of decreasing as separation was 
approached, actually increases. Schubauer’ and Lud 
wieg,* using direct measurements of shear, found that 
c, decreases in rising pressures. However, when the 
two-dimensional momentum equation is used to obtain 
cy from Schubauer’s data, rising values of c, are ob 
Ludwieg has suggested that lack of two-dimen 


On the other 


tained. 
sionality is the cause for the discrepancy 
hand, Newman,'? Goldschmied,'* and Bidwell'' sug 
gest that the discrepancy is caused by the momentum 
equation being in error through the omission of terms 
containing the turbulent fluctuations 

During the course of the present experiments, it was 
found that agreement with the two-dimensional mo 
mentum equation was obtained when the pressure was 
constant, but when the pressure was rising, definite 
disagreement was obtained. For our second rising 
pressure distribution, with near-separating profiles, the 
disagreement was very great. The skin-friction term 
was negligible compared with the d@ dx and |(// + 
2] (6q) (dp dx) terms. (U7 + 2) 2] 
(@ g) (dp dx) term was 50 to 100 per cent greater than 
the d@ dx term. If ¢, 
momentum equation, a large negative value would 


)) 


Further, the 
had been computed from the 


have been obtained. 

An attempt was made to obtain agreement with the 
momentum equation by taking into account the fluc 
tuation terms as suggested by Newman, Goldschmied, 
and Bidwell. However, this procedure led to correc 
tions far too small to explain the observed discrep 
ancies. In the case of equilibrium profiles, with simi 
lar conditions existing at each station, the fluctuation 
terms in the momentum equation are either zero or 
very small. The reason for this is that these terms 
usually occur as an x derivative of an integral of the 
fluctuations. With these 
integrals are substantially constant. 


similarity in conditions, 

An investigation was made to see if the discrepancy 
was caused by nonlinearity in the instrumentation, 
which would give incorrect mean velocity measure 
ments in the presence of turbulent fluctuations. Mean 
flattened hypodermic 
Turbulence measurements 


velocities were measured by 
probes and by hot-wires. 
were also made by hot-wires. The hypodermic probes 
were calibrated against angles of attack and yaw. 


Knowing the levels of the u fluctuations and estimating 
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the levels of the v and w fluctuations, estimates were 
made of the effect of the turbulence on the mean ve- 
locity readings of both the hypodermic probes and the 
hot-wires. Reasonable agreement was obtained ex- 
cept near the wall where the turbulence levels were 
relatively large. It was found that the corrections 
had insufficient effect on the integrated values of @ and 
6* to explain the momentum discrepancy. 

From the above it was concluded that neither the 
additional terms proposed by Newman, Goldschmied, 
and Bidwell nor instrumentation errors could account 
for observed disagreement with the momentum equa- 
tion. The evidence pointed to the possibility that the 
flow was not two-dimensional. Apparently, the flow 
was spreading laterally in such a way as to tend to thin 
out the boundary layer. If w is the local component 
of lateral velocity parallel to the plate (w is assumed 
constant with y), then a term (6*/L’) (Ow,/Ozs) must 
be included in the momentum equation. Calculations 
showed that the lateral flow required to explain the 
discrepancy appearing in the momentum equation was 
At the sides of the tunnel, a divergence angle of 
Attempts were 


small. 
about 1|'/.° would have been sufficient. 
made with yawmeters to measure the lateral flow, but 
these were unsuccessful because the lateral turbulent 
fluctuations were so large (greater than + 10° in flow 
inclination) that they masked any small shift in mean 
flow direction. 

It will be remembered that the 
structed so that the boundary layers on the side walls 
are removed. It was conjectured that this removal 
might be causing the flow to diverge laterally. The 
slots were taped up, and the discrepancy was reduced 
This was surprising because now 


tunnel is con- 


but not eliminated. 
the wall boundary layers should have caused the flow 
on the test floor to converge and reverse the discrep- 
ancy. Next, it was noticed that a crack, about 1/32 
to 1/16 in. wide, existed where the walls joined the 
floor. This was taped up, and the discrepancy was 
again improved but not eliminated or reversed. At 
one section of the tunnel, this last change caused a 
previously well-behaved boundary layer on the floor to 
separate. This indicates that these fine cracks had been 
removing sufficient low-velocity air to prevent sepa- 
ration. An analysis of this situation showed that, when 
a leak such as these cracks exists, there is a strong 
tendency for the swiftly moving central stream to 
travel straight down the tunnel but for the slowly mov- 
ing air next to the floor and walls to move laterally 
toward the leak. 

Next we considered the possibility that the boundary 
layer air on the floor was moving laterally to the walls 
and then up through the wall boundary layers to the 
outlet flaps on the top of the tunnel, which were used 
to control the pressure distribution. Measurements 
of the wall boundary layers indicated they were not 
growing, even though the wall slots were closed. This 
strengthened our belief that Jow-velocity air was mov- 
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ing through the boundary layers to the outlet flaps. 
To prevent this movement, the ends of the outlet. flap 
openings next to the walls were closed. This was 
finally sufficient to reverse the discrepancy. Once the 
discrepancy was under control, it was an easy matter 
to adjust conditions so that the two-dimensional mo- 


mentum equation was satisfied. 


After working to cure this discrepancy as long as we 
did, we came to have a great respect for the ease with 
which air can move laterally in boundary layers sub 
jected to adverse pressures gradients. None of this 
trouble appeared when the pressure was constant. 
After having had this experience, it is easy to believe 
that the results of Mangler, Wieghardt, and Schu- 
bauer should be interpreted as indicating that their 


flows were not two-dimensional. 
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General Instability of Simply Supported 
Rectangular Plates 


J. R. M. RADOK* 


Aeronautical Research Laboratories, Melbourne, Australia 


SUMMARY 


rhe problem of lateral bending of thin plates, stiffened by line 
reinforcements, is reduced to the limiting case of uniform plates 
under triangular loading, distributed over a strip, when the 
width of the strip tends to zero. Using the basic solutions, aris 
ing from this limiting process, the characteristic equations for the 
compressive buckling loads of rectangular plates, reinforced by 
stringers, by ribs, by distributed stringers, and by discrete ribs, 
are deduced. In all these cases closed expressions in terms of 
elementary functions are obtained for arbitrary spacing of the 
reinforcing members, and it is shown that the case of discrete 


stringers and ribs cannot be solved in a similarly simple manner 


NOTATION 


a,b = dimensions of plate 
thickness of plate 

k = 1 + (12/,/h*) (1 — v? 

lateral loading of plate 


iM = arbitrary coefficients for stringer and rib 
loading functions 
= lateral deflection of plate 


ae = coordinates 

vi, = positions of ribs or stringers 

A = area of stringer 

b = flexural stiffness of reinforcing member 

D = flexural rigidity of plate 

/ = second moment of area of distributed 
stringer 

Vr, Mery, Wy = moment resultants 

N; = normal compressive force per unit length of 
plate edge 

p = a/b 

Y = B/bD 

; ji,i= 

6 = kKronecker’s delta—i.e., 6;; = 

(0,7 4) 

6 = A/bh 

v = Poisson's ratio 

A = discontinuity of expression in brackets 


(1) INTRODUCTION 


— PROBLEM OF GENERAL INSTABILITY of rec- 
tangular plates with reinforcing members parallel 
to their edges was first investigated by Timoshenko in 
1913. This work is reproduced in the same author's 
book on Elastic Stability.'. Since then several other 
authors have dealt with this problem, among whom 
mention should be made of Cox and his coauthors?~* 
in England, Seide and his coauthors® '° in America, 
and of Lokshin’ in Russia, whose papers were not ac- 
cessible to the author of this report. 
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In the references quoted above, results have been 
deduced by means of the energy method. This method 
is known to lead always to infinite systems of equations, 
and in the case of reference 1 conclusions with regard 
to the buckling stresses are deduced by considering 
progressive approximations, obtained by neglecting all 
but the first few unknowns of these equations. On 
the other hand, it is shown in references 3 and 4 that 
the variables in these equations may be rearranged 
into infinite series that are easily summed and that in 
this way the infinite system may be reduced to a finite 
one involving new variables. <A similar approach has 
also been adopted in references 9 and 10. 

The method of the present report is equivalent to that 
applied in reference 5 to the problem of general insta 
bility of indirect 


approach of the energy method and leads immediately 


circular cylinders. It avoids the 
to a determinant of order equal to the number of rein 
forcing members parallel to one of the sides of the plate. 
It also shows that the case of a plate with stringers, as 
well as ribs, cannot be reduced to a closed expression 
This last fact was first suggested in reference 4 and also 
in reference 7, which partly deals with the buckling 
problem of plates clamped along all four edges; it is 
easily shown that this latter problem is equivalent to 
the former. 

The present approach deduces the characteristic 
equations directly from the equilibrium equation for 
a plate in compression. It will be seen below that in 
the cases of ribs or stringers it leads directly to the buck 
ling equations of references 3 and 4. Apart from these 
cases, a detailed discussion is given of the case of 
stringers and ribs and of the case of distributed stringers 
and discrete ribs which is of importance for aircraft 
construction when the stringer spacing is often so close 
that each buckling wave comprises several stringers. 


1.1) Discussion and Interpretation of the Method of 

Solution 

Before stating the basic assumptions, a short physical 
interpretation of the mathematical procedure will 
greatly assist the understanding of this report and of 
the significance of the assumptions made (cf. also refer- 
ence 5). 

In order to ensure the compatibility of displacements 
and stresses in a plate reinforced, for example, by 
stringers, certain relations must be satisfied between 
the shear stresses at adjacent edges of neighboring 
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panels, or, in other words, these stress resultants under 
go discontinuous changes at the stringers. As_ the 
equilibrium equation to be used here represents the 
condition of lateral equilibrium of a plate element, the 
stringers, if idealized into lines, must be loci of infinite 
pressures. This was the basic idea of the method of 
solution of reference 5. Accordingly, the homogeneous 
differential equation 


Viw + (N,/D) (07w/Ox7) = 0 (1.1.2) 


which is relevant for the case of a simple plate without 
lateral loads, would have to be replaced by the non- 


homogeneous equation 


Viw + (N,/D) (0?w/Ox?) = g/D (1.12 


where g is infinite along the lines of attachment of the 
stringers while it is zero elsewhere. 

However, the approach of reference 5 suffers. still 
from the disadvantage that it is not analytically rigor- 
ous. In fact, the differential equations of that report, 
as well as Eq. (1.1.1), are based on the assumption that 
the displacements have continuous derivatives up to 
and including the fourth order. For this reason, the 
use of a conditionally convergent series for g which is 
divergent at the reinforcement is not strictly justi- 
fied. In the present report this formal objection has 
been overcome by means of a limiting process that is 
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applied to the solutions corresponding to local lateral 
loading of finite magnitude. 

The above loading g is thus introduced by means of a 
suitably chosen Fourier series, which is zero every 
where except along a strip surrounding the stringer 
where it represents a finite triangular distribution of 
pressure. (Fig. 2.) Eq. (1.1.2) may then be solved by 
substituting for w a Fourier series with unknown co 
efficients, the values of which can be determined in 
terms of the coefficients of the earlier series. In this 
way one obtains the solutions of Section (2.2). 

In Section (2.3) the above-mentioned limiting proc 
ess is applied to the solutions of Section (2.2), which 
leads to the basic solutions to be used in Section (3) for 
the solution of the buckling problems. These basic 
solutions may be interpreted as sources of shear stress. 
So far, solutions of this type have only been used in plate 
theory when dealing with point loads.*® 

Finally, a few remarks will be made with regard to 
the analytical properties of the solutions of Sections 
(2.2) and (2.3). It is easily seen that the solutions of 
Section (2.2) satisfy all the necessary conditions with 
regard to continuity of derivatives which are demanded 
by the differential equation. The basic solutions of 
Section (2.3), on the other hand, do not belong to the 
same class of continuous functions, because certain 
of their third derivatives are discontinuous. However, 
this property of discontinuity is essential for the pur- 
pose of satisfying the “‘internal’’ boundary conditions 
so that the basic solutions are suitable for the solution 


of the problems in hand. 


1.2) Assumptions 


Only those assumptions that refer directly to the 
reinforcing members will be stated explicitly. Some 
of these are not essential and have only been introduced 
in order to simplify the presentation of the theory. For 
example, it is not difficult to relax the conditions that 
the flexural axis of the reinforcing members should lie 
in the middle plane of the plate and that their torsional 
deformation may be neglected. In reference 3, which 
deals with the case of ribs, this last assumption has not 
been made. 

(A.1) The theory of thin plates is applicable. 

(A.2) The plate is simply supported along its 
edges. 

(A.3) The flexural axes of the reinforcing members 
lie in the middle plane of the plate. 

(A.4) The reinforcing members do not possess tor 
sional stiffness. 

(A.5) The stringers or ribs may be idealized into 
lines. 

(A.6) In the case of stringers and ribs the former are 
so closely spaced that each transverse buckling wave 
includes several stringers that may hence be dis- 


tributed. 
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INSTABILITY OF 


2) EQUILIBRIUM EQUATIONS AND BASIC SOLUTIONS 


hroughout this report, whenever possible, use will 
be made of the notation of references | and 6. In 
particular, stringers will be reinforcements in the direc 
tion in which the plate is compressed, while ribs will be 
transverse reinforcements (see Fig. | As far as the 
deduction of the equilibrium equations is concerned, 
the work of the above references will be fully utilized, 
and only the final stages will be reproduced in order to 
illustrate the modifications caused by the presence of 


distributed stringers. 


2.1) Equilibrium Equations 
In terms of the moment resultants, the equilibrium 
equation of a plate subject to lateral loads and com 
pression in the Ox direction is given by° 
OM O° M,, O° .1, Ow 
— 2? +- - "ie A ; 
Ox" Ox OV OV 


For an isotropic plate the following relations hold be 
tween these resultants and the lateral displacement w, 
provided the last is small 
WV, = —D[(0°w Ox?) + v(O°w Oy I] 
M, = —D[(O°w, Oy”) + »(O0°w/ Ox?) | Gig 
M,, = D(1 — v) (O*w/Ox Oy) 
Hence, Eq. (2.1.1) becomes 
ary O'w 
; Ge 
Ox" OY 
‘Wy N, O-w q 
dvi =D dx?) =~D 


N, O-w 
D ox? Ox! 


(2.1.3) 
In the case of distributed stringers, one is effectively 


dealing with an orthotropic plate and Eqs. (2.1.2) 


must be replaced by 


M, = —D[k(0?w/Ox?) + v(O0°w Oy") |) 
M, = —D|(O?w/Oy?) + v(0°w Ox") | (2.1.4) 
M,, = D(1 — v) (0?w/Ox Oy) 
where 
P= 1 + (12.8) (11 — + (2.1.5) 
Instead of Eq. (2.1.3), one then obtains, from Eq. 
2.1.1] 
O'w N, 0*w ma; 
Viw + (Rk — 1 = + 
Ox! D ox* ox" 
> ay 77 N, O-w = q (2.1.6) 
Ox? Oy" Ov! D ox* D 


2.2) Solutions for Lateral Loading over a Finite Width 
Strip 
Following reference 5 the desired loading would have 
to be introduced by means of the following loading 


functions: 


SIMPLY SUPPORTED RECTANGULAR PLATES 1] 


q=D> > sa:sin sin ~~ sin —- 2.2.1a) 
m=1 n=1 a D 0D 


for a stringer at the station y,, 


T 1A 


_ - . umTV . mqmX 17 
gq DS a sin — sin sin 2.2.2a) 


fora rib at the station \ As indicated in the introduce 
tion, these loading functions, while corresponding to 
the basic solutions of the next section, lead to solutions 
of Eqs. (2.1.4), which have certain discontinuous third 
derivatives. For this reason the ensuing analysis 
would not be justifiable, since the solutions of this 
equation must have continuous derivatives up to and 
including the fourth order. 


The loading functions, Eqs. (2.2.la) and (2.2.1b), 


will now be replaced by the functions 


i one . Manx . nwy, . Nerv sin* nd 
g D> YS sn) sin sin sin - 
’ 1 l a D dD °6 
p+) I 
a . ATV . mm Mam sin> m6 
q D> 2 +.,;9 — in sin 
l | ) a ad 42-0 
) ) ) 


which are easily seen to correspond to triangular load 
ings over strips of width 46 surrounding the respective 
reinforcements. The widths of these strips must, of 
course, be restricted in such a way that the loading 
stays within the region of the plate (see Fig. 2 
It is easily seen that the series (2.2.1) and the series 
(2.2.2) are absolutely and uniformly convergent in the 
region of the plate, provided 6 has a finite value. It 
is also readily verified that, in the limit when 6 — 0, 
these series tend formally to those given by Eqs. (2.2.1a) 
and (2.2.1b). 

First consider the case of the isotropic plate, subject 
to the loadings in Eqs. (2.2.1) and (2.2.2). Substi 
tuting for win Eq. (2.1.3) the series 

> C,, sin coin sin : a 2.2.3) 
l l a 0 


and comparing coefficients of sin (mmx a) sin (ny >), 
one finds 
Sm) Sin (nav, b sin* 6 


Com = j (/m\* n\*) ml \ n°6 
ry ( ) +G)i (") " D 


\a 


for the case of triangular loading on a x-wise strip at 


y,, and 


r,, Sim (marx; a sin* 726 


Con = j fm : n\*| m V 176 
ri) +) -(") " D 
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for the case of triangular loading on a y-wise strip at 
x), where in both cases the width of the strip is 46. 

In the case of the orthotropic plate with ribs, one 


finds from Eq. (2.1.6) 


"By J 
rnj Sin (Max ;/a) 
(7) (") bs (h—1 a) my , Nz 
il " , woot} - 7" 
a b/ § a a D 
sin” mé 
oe (2.2.6) 
m*5* 
where obviously 
k—1>0 (2.2.7) 


Substituting for C,,, in Eq. (2.2.3) the above values, 
one obtains the relevant solutions for the plate with 
strip loading of the indicated type. 


(2.3) The Basic Solutions 


The required solutions for the problem of general in- 
stability of plates with line reinforcements can be ob- 


> # 
( mn 


[ti + Gt 


where 


MTX 
Sin 
a ) 


nTry 


W = (2.3.4) 


ps 7. Cus ae 
1 


m=1%n 


These expressions for w are easily seen to have the 
properties required by the present problem. For ex- 
ample, substituting for C,,, the values in Eq. (2.3.1), 
it is seen that the derivative 0*w/Oy* has a discon- 
tinuity at vy = y;, since the corresponding solution con- 


tains series of the form 


——_— NTV . NY; ce 
p> cos sin (2.3.5) 
n=1 1” b b 
The discontinuity of this series at the station y = y; is 
obviously 
w/2 (2.3.6) 


Similar reasoning may be applied to the remaining se- 
ries. 

In the next section the characteristic equations for 
the general instability of stiffened plates will be ob- 
tained by satisfying the “‘internal’’ boundary condi- 
tions that are known to relate the deformations of the 
reinforcements to the above discontinuities. 
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tained from the solutions of Section (2.2) by means of 





the same limiting process that transforms the load- é 
ings in Eqs. (2.2.1) and (2.2.2) into the basic loadings 
in Eqs. (2.2.la) and (2.2.2a). In fact, letting 6 ~ 0, | 
one obtains, finally, for the isotropic plate, 
C. Sng SiN (nay;/b) 
- f(y’ Pe 4p ; =)" ? 
7 —_ T~ 
a b/ f a D 
: Ww 
(2.3.1 
Vi 
p ; m 
for a Stringer at Vj, th 
ta 
Cc! nj SIN (Max) /a) re 
™ j/m\* n\?)\? m\? _N, 
1: + . — Tr 
\\a b/ S a D Uy 
(2.3.2) 
for a rib at x;, while for the orthotropic plate with a 
rib at x;, 
wl 
; tr: 
rnj SIN MTX; A 
: (2.3.3) 
+ (b : (") (“) _N; 
(a — 5) = i 
a a D 
(3) THE BUCKLING EQUATIONS 
(3.1) Case of Discrete Stringers or Ribs 
The basic solutions of Section (2.3) contain the arbi- 
trary constants s,,; and 7,,;, which are required to satisfy 
the internal boundary conditions at the reinforcing 
members. On the basis of the remarks at the end of 
the last section, these conditions are On 
obt 
O*w A,N, 0*w O*w 
2, + : = —DA| — the 
Hi Ox » Gt Jou "lle Pe for 
(3.1.1) (3. 
for a stringer at y,, and 
O*w O*w 
Oy* jr =z; Ox* Jeux; 
I 
for a rib at x,. In these formulas, A indicates the | me 


discontinuity in the respective derivatives that are the 

only discontinuous parts of the relevant shear re- 

sultants. > 
First consider the case of f stringers. The relevant ; 

basic solution is obtained from Eqs. (2.3.1) and (2.3.4) 

by summing with respect to 7 over all stringers. Sub- 

stituting this solution in the f conditions, Eqs. (3.3.1), 

corresponding to the stations y,(k = 1, 2,...,/), one | 


finds 


ns of 
load- 
lings 


—> () 


rbi- 
sfy 
ing 


of 


he 
he 
re- 
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_ N,B?m*b? a sin (mmy;/b) sin (nry,/b) 
6. + 2 (sm ~ % -— arusenones = 0 
Dr? 1 3m? + Bn?}? — (N,B?m*b?/Dx?) 


which is the condition for the existence of nonzero 
values for the constants s,;. Eq. (3.1.3) is a deter- 
minant of order /, and it is easily seen to agree with 
the results on page 374 of reference 1, which were ob- 
tained by means of the energy method. In fact, the 


result given there may be rewritten as 


N,6?m*b* | 


Lmn (Mm? + B*n-)* — 
{ Dr? f 
N,B?m7b? |. may 
2 Xs xm — 6, — 2 sin : > Quy 
] Dr* bh p l 
Try; 
jin” —=0Q@ 

b 


which is an infinite system of equations for the arbi- 


trary coefficients d,,,.. Dividing the above equation by 


7 _ ») * yO, 1 a k (2 ° 
(i, = ee ae 0, = J. * od. 1 0d) 


the coefficient of a,,,, multiplying by sin mmy,/6, and 
summing with respect to m, one obtains a system of | 
equations in the new variables 

; . PRY: 

>> mp sin ; y Cao 2 


pal ‘ 
the determinant of which is Eq. (3.1.3). Thus it is 
seen that the constants s,,; of the basic solutions are 
proportional to the above variables. The same process 
of summation has been used frequently in references 
2-4. 

Next, the case of g ribs will be investigated. The 
same procedure may be applied as in the case of string 
ers, using this time Eqs. (2.3.2), (2.3.4), and (3.1.2 
The characteristic equation is now 


_ MTX; | MTX, 
} sin sin 
P See 4 a | 
by, + 27 B'n'! > a 0, (Vj, 1 = Ss 3 (3.1.4 
Mewt se 4 sy cashes N,B?m?*b? 
5m° + 2B? — - 
Dr’? 
On page 379 of reference 1, a related result has been there if one puts” = 1. It can also be shown that the 


obtained using again an energy method and assuming 
that m = 1. By a process analogous to that indicated 
for the case of stringers, it can be shown that Eq. 


(3.1.4) is equivalent to the special result obtained 


Finally, consider the series, figuring in the determinants (3.1.3) and (3.1.4 


manner used in reference 3, and the corresponding sums are 


coefficients 7,; of the basic solutions are proportional 


to 


. pwx 
> » Ay», Sin f 


p=1 a 


). These series can be summed in the 


nTy; Vi ; v; 

sin —— sin = sinh 7B —- )} sinh ) B 

> b ) . ~ ( b/f ( bf 
= 7 B2m2b? 0 sinh +B 
tn? + pnt}? — N,B?m*b S68 B sinh rf 

Dr?* 
ms 7 
sinh J rA(1 a - ) sinh SaA Yel | ‘ 
| b/I ( bN | (fory,; €y,) (3.1.5) 


A sinh 7A 
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where 
m- N,b? 
C* = , = —, A*=c(c + 20), B* = c(c — 20 3.1.6) 
B° 1{Dr? 
marx marx; , : , Y : ‘ x SEs 
sin sin sin r(F + G) {1 - ) sin 7w(F + G) 
= a a 7 a ad 
— : : N,B?m*b? SIG (F + G) sin (Ff + Gj)r 
jm? + B'n?}? B 
Dr- 
. ; . : ” ee 
sin a(F G ( ) sin rik —G 
a a >. | 4 
(F G) sin (F G)r 
where 
“é a’ N, = " vo 
fe = , G=F n-B?, x, 4X a.18 
tr? D 


Thus, the determinants (5.1.5) and (5.1.4) ean be written in closed form, using Eqs. (3.1.5) and (3.1.7), respec 
tively. The final characteristic will be independent of 7 or m—1.e., of the wave lengths at right angles to the rein- 
forcing members. The equations deduced here are somewhat more general than those of references 3 and 4, which 
assume constant stringer or rib spacing. 


3.2) Case of Piscrete Stringers and Ribs 


This problem is referred to in reference + and also in reference 7, since the analytical difficulties that it presents 
are of equal nature to those encountered in dealing with the buckling problem of a plate clamped along all edges. 
The purpose of the present discussion is to provide a proof that the characteristic equation cannot be presented in 
closed form or, quoting from reference 4, that the ‘‘difficulties in respect of mathematical summations’ are irre- 
solvable. The proof given here is indirect and leads to the conclusion that summation will only be possible if all 
the reinforcing members are parallel. 

Using Eqs. (2.3.1), (2.3.2), and (2.3.4) and assuming / stringers and g ribs, the relevant basic solution is 

WEN... MEY 
sin sin 


a ) = . nay — . Max;) 
>i sin “i+ > 7,,8 


Es | 
m=1n=1 m\* ney m\* N, i=! b s a I 
AC+O)-@% : 
Wa b/ J a D 
Substituting this solution in Eqs. (3.1.1) and (3.1.2), one finds the following set of equations 
> va Sp b B, ("") _ ALN, (=) Pe ra | + |B, ("") _ ALN, (=) . > R.i* gin ME = yy 
i=1 { 2 a h a f a h a j=] a 


> B, (=) sin “"™S," + Yr, )D “+ B,("") Aj =o 
) re 9 


s=1 b b f 
(em bt cad. PEE 2... Be Loss @. BH 1.2...., 9 3.23) 
where 
| nny; , NTY, _ max; . marx 
sin sin sin sin ve 
) ) ; a a ($3.25 
i = + Af= 2) 
é m at ’ 7 
ant wn m= A 
, is NT; “di ‘ MAX | 
RS =? sin => = 2 1 3.24 
n=1 4 b m=1 Am, a 
Am, = m*} (ma)? + (n/b)?}? — (m/a)? r2(N, D o2.0 


In agreement with the deductions at the end of the in order to achieve this state, the new variables R,, and 
preceding section, the first f equations and the remain- S,, had to be introduced so that these equations are no 
ing g equations of Eqs. (3.2.2) (for constant m and”) longer in terms of the same variables. Thus there 
are independent of ” and m, respectively. However, arises the question whether a further summation can be 
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INSTASILITY OF SIMPLY 


applied to these equations so that one obtains a finite 
set of equations in consistent variables. 


Rewrite Eqs. (4.2.2) in the following manner 


Dh 
| + A,,'{ 
} mr \* A,N Nim \* + 
a) = (S) 
a il a 
D WLIW, (29 9") 
R,, sin = 2 
a 
, 
S, sin asad f, Da + A a im 0 
b PB (nm bd)! { 


assuming that there is only one stringer at y, and one 


rib at x,. It is now easily seen to be impossible to re 
duce these equations to two equations with consistent 
variables. In fact, 
R,, in the second set of equations in conformity with 


having introduced the variable 


Eqs. (3.2.4), this variable can be eliminated by means 
of the first set of equations. But this process will not 
lead to a closed expression, since S, already involves 
as, Gt = 1,2... 24 & A similar reasoning applies 


,—) 


SUPPORTED REC 


TANGULAR PLATES 15 


reduce to a closed expression if either all s,, or all r, are 
zero--i.e., if only parallel reinforcing members are 
present. 

The physical interpretation of this result is largely 
contained in the final conclusions of Section (:.1 
It has been seen there that in the cases of stringers or 
ribs the buckling equations involve the number ot 
waves in the directions of the reinforcing members 
In the case of stringers and ribs the wave lengths in 
the x and in the y directions are of equal importance 
If this were not so, it would follow that the case of only 
parallel reinforcements would also be independent ot 


the corresponding wave lengths 


3.3) Case of Distributed Stringers and Discrete Ribs 


It has been seen in the previous section that it 1s 
not possible to obtain a closed expression for the char 
acteristic equation in the case of a plate with stringers 


and ribs. However, in many practical cases, the 


stringer spacing will be sufficiently close as com 


pared with the rib spacing to justify the introduction 


of the assumption (A.6 In that case one obtains, as 


with respect to S,. Hence, these equations will only _ before, using Eqs. (2.3.5), (2.5.4), and (5.1.2 
mmx; | mmx 
sin sin 
‘ , — a a 
6;, t+ 2kRy Bn’ 7. 0 } = 1.2 2 3. 1 
1 N,B?m?b? 
}m? + n27B242 + (Rk m* — 
Dr 
rhe sums of the series in Eqs. (5.5.1) are 
. mex, . maxX 
sin sin 
~ a a vs ‘ \ , \ 
s ; — 2 sinh 7(Q — P { | ) sinh 2(Q P 
n= ; ee V,8°m"b 1kP 1 a 
Vim? + n2B2]2 + (k — 1)m' — : 
Dr 5 QO — P)”* sinh x(Q Y 
: 2 : ‘ ] . | 
sinh 7(Q + F {1 sinh r(Q + I 
ad ad ” : >) 
(0+ P) *sinh r(Q0 + P a 
where 
: Va n'B) 1 / N,a? 
Pp? = ) on3g? cll, », O= (” 3° ) (X%, WA Bide 
ik? | Dr? I k k 2Dr* 


(4) CONCLUSIONS 


The characteristic equations for the general insta- 
bility of rectangular plates with stringers or ribs, with 
stringers and ribs, and with distributed stringers and 
discrete ribs have been deduced directly from the 
equilibrium equations. The method used leads in all 
but one of these cases to a determinantal equation the 
order of which equals the number of discrete reinforc- 


ing members. For the case of stringers and ribs it 


has been shown that the present approach leads to an 
infinite determinant. 

Some of these buckling equations have been pre- 
viously obtained by the use of the energy method.'~‘* 
However, comparison of the present work with that 
of references 3 and 4 shows that the energy approach 
leads to more involved calculations. 

Another advantage of the present approach is that 
the basic solutions of the equilibrium equations. ob- 





116 JOURNAL OF THE AERONA(lU 


TICAL SCIENCES FEBRUARY, 954 
tained in Section (2.3) have a simple physical inter- Cox, H. L., and Riddell, J. R., Buckling of a Longitudine 
pretation in that they represent the mathematical Stuff ned Flat Panel, Aeron. Quarterly, Vol pp. 225-24 
formulation of discrete reinforcements. For this rea- asin cabis ' ; 
son, the analytical process leading to the characteristic et eae en jaadsivasa dagen Lory ite Aine 

: ; drical Shells, College of Aeronautics, Cranfield Rep. 61, June 


equations may be directly related to the physical as- 
pects of the buckling process, and a clearer under- 
standing of the mechanism of general instability thus 
becomes possible. 
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Conference on Supersonic Flow 


An unclassified conference on Supersonic Flow will be held in the Mathematics Build 
ing, University of Maryland, from 9:30 a.m. to 5:00 p.m. on March 19-20, 1954. 

Information concerning the conference may be obtained from Dr. E. L. Resler, Associ 
ate Research Professor, Institute of Fluid Dynamics and Applied Mathematics, University 
of Maryland, College Park, Md. 
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Survev of Friction Coefficients, 


Recovery 


Factors, and Heat-Transfer Coefficients for 
Supersonic Flow’ 


JOSEPH KAYE? 
Massachusetts Institute of Technology 


SUMMARY 


\ brief survey is presented of the progress made in the last 4 


years on theoretical and experimental work on friction coefficients, 


recovery factors, and heat-transfer coeflicients for supersonic 


flow of air 


NOMENCLATURI 


heat-transfer area 

local skin-friction coefficient 

Stanton Number of Van Driest 

specific heat at constant pressure 
D = inside diameter of pipe 

local apparent coefficient of friction, 2g7/ pV? 
( = flow per unit area 

acceleration given to unit mass by unit force 


] = heat-transfer coefficient, (g/A = 
V/ = Mach Number, V/V gkR7 

Vu, = length Nusselt number, AL /k 

Py = Prandtl Number, Cpyug/d 

L = static pressure 


= rate of heat transfer 


= recovery factor, (law — h to — 1 
Re, = diameter Reynolds Number, DG u 
Re, = length Reynolds Number, LG u 
St = Stanton Number, 4/cpG 
l - temperature, °F. abs 
temperature, °F 
| = véE locity 
= density 
r = thermal conductivity 
lu = viscosity 


T = shear stress at wall 


Subscripts: 
0 refers to stagnation conditions 
refers to adiabatic wall conditions 
1 refers to mean-stream conditions 
refers to free stream conditions 
refers to flow of incompressible fluid 


refers to wall conditions 


Presented at the Aerodynamics Session, Twenty-First Annual 
Meeting, IAS, New York, January 26-29, 1953 

* The assistance of Prof. T. Y. Toong and George A. Brown in 
the preparation of Figs. N, O, P, and Qis gratefully acknowledged 
This survey report was written mainly under the sponsorship of 
the Office of Naval Research of the United States Navy under 
Contract Number N5-ori-07805. Permission has been received 
to reproduce the figures taken from the JOURNAL OF THE AERO 
NAUTICAL SCIENCES, the reports of the National Advisory Com 
mittee for Aeronautics, and the reports of the U.S. Naval Ord 
nance Laboratory 

+ Associate Professor of Mechanical Engineering 


INTRODUCTION 


Sf ipn IMPORTANCE OF DESIGNING AIRCRAFT for safe 
operation at higher and higher supersonic speeds is 
apparent to all. Twenty years ago designers of air 
craft were concerned with the effects resulting from an 
aircraft passing through the speed of sound —‘‘the 
sonic barrier.” But having shown that it is feasible 
to design an aircraft that can start from rest and exceed 
sonic speed safely, the designer is at present concerned 
with safe operation of aircraft at supersonic and hyper 
sonic speeds. 

Safe operation of supersonic aircraft implies safety 
for both the inhabitants of the aircraft and its equip 
ment. It is evident that, as the speed increases to 
five or ten times the sonic speed, the function of the 
human beings on board must be reduced considerably, 
since the human being represents a poorly designed 
servomechanism at these high rates of action and reac 
tion. Hence, as the speed increases, more and more 
of the normal human operations will be taken over by 
electronic gadgets. But here the designer is faced with 
an obstacle similar to the ‘‘sonic barrier.”’ 

If the speed of the plane is sufficiently high, its sur 
face temperature may reach a value of the order of 
|,000°F. or greater, resulting from ‘‘aerodynamic heat 
ing.’ Hence, the designer can foresee a ‘‘thermal 
barrier’ requiring a tremendous amount of internal 
cooling in this aircraft in order to keep the inhabi 
tants, the aircraft, and the electronic equipment in safe 
operating condition. The problem of safe operation 
of electronic equipment is made even more difficult 
when it is realized that most of this equipment has 
been designed only for operation at low ambient tem 
peratures and not for continuous operation at 1,000°F 
or at greater temperatures. 

The designers of supersonic aircraft require an ex 
tensive knowledge of friction coefficients, recovery 
factors, and heat-transfer coefficients over a_ wide 
range of Mach Numbers for both laminar and turbu- 
lent boundary layers. It is intended here to present 
a survey of available data and theory for supersonic 
velocities and then to discuss two research programs 
in this field now under way at M.I.T. 

Because of the large number of papers published 
in this field in the last few vears, only some of the latest 
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TABLE 3 


SUMMARY OF ASSUMPTIONS USED IN THEORETICAL WORK ON LAMINAR BOUNDARY LAYER ON A FLAT PLATE 


Specific Thermal Prandt) Pressure 
Author Date Density Heat Conductivity Viscosity No, Gradient Radiation 
Blasius 1908 
Pohihausen 1921 c r < « c o r 
1/2 1/2 
Busemann 1935 {(T) ~ T ~T 1 0 u 
Von Karman 77 77 
1939 {(T) < ~ ~ 
and Tsien : . 
Crocec 1941 (T) ‘ ee ~T” 0.725 ( 0 
(for w = 0.5,.75,1,.0,1.25) 
Brainerd 768 76 
etd 1942, f(T) P ~ 1% i 0.733 ( r 
and Emmons 
1 1 
Crocec 1946 f(T) c a TF ~T 
Cop 4 89 .89 
witboos 1943 f(T) ‘ ~ 0.76 
Hartree 
~ 6 1 
Chapman’ ‘ 1 1 
1949 f(T) c ~ T ~ 7 ( 
and Rubesin ’ wor - 
28 
Lighthill 195 {(T) « {(T) (T) 0.7 (Vv) 
Van Driest’ 1951 {(T) ‘ 1(T) Sutherland 
Rule 
Klunker and ,... nT nt 1 P a 
Matioest © 195 (T) (tT) {(T) (T) (T) 
Young and 
1952 (T f f 
Seeall 5 ) f(T) (T) f(T) f(T) 
Moore? 1952 "(T) f(T) (T) f(T) (T) 
= constant 





Mach Frictior Recovery 
Dissociation No. Coefficients Factor Remarks 
¢ = 0.664Rej! /? 
664Rej!/? rs pr!“ v 


Arbitrary surface temperature 
Arbitrary velocity of main 


stream and arbitrary surface 


temperature 


t(T,P) 


<> variable examined 


—= variable not examined 





ones are reviewed; no attempt is made to cover the 
field completely. In addition, although the problem 
of transition is closely associated with the above prob- 
lems for laminar and turbulent boundary layers, it 
was decided to exclude the large number of papers 
treating this difficult problem. 


OBJECTIVES 


The objectives of the present discussion are three- 
fold: 

(1) To present a brief survey of the state of the art 
on theoretical and experimental work on friction co- 
efficients, recovery factors, and heat-transfer coeffi- 
cients for supersonic flow. 

(2) To present a brief description of the program in 
supersonic flow of air in a tube at M.I.T. under the 
sponsorship of the Office of Naval Research. 

(3) To present a brief description of the program on 
supersonic flow of air over a flat plate under way at 
M.1.T. under the sponsorship of the Office of Air Re- 


search. 
REVIEW 


(A) Theoretical Work on Laminar Boundary Layer 


The analytical investigation of the laminar boundary 
layer began in 1908 with the thesis of Blasius, who de- 
termined the velocity profiles for incompressible flow 
on a flat plate with zero pressure gradient. In 1921 
Pohlhausen used these velocity profiles to calculate 
temperature profiles in the laminar boundary layer on 
a flat plate. Since then numerous papers have been 
published which extended and modified, first, the basic 
assumptions used in the analysis; second, the mathe- 


matical techniques of reduction of the partial differen- 
tial equations; and third, the methods of obtaining 
numerical results. The variety of assumptions used 
in some of these papers is seen from the summary in 
Table 1. 
| are given in detail in an excellent review by Kuerti,! 


Most of the references before 1950 in Table 


which also contains the references to the earlier reviews 
of Lewis and of Rubesin and Johnson. 

Since the papers listed in Table 1 represent only a 
fraction of the available papers, it is obviously possible 
to give here only a few samples of theoretical work. 
These will be selected from the most recent papers. 

(1) Skin-Friction Coefficient for Laminar Boundary 
Layer.—The calculated values of local skin-friction 
coefficient, c;, for flow of air over a flat plate are com- 
pared in Fig. A (Moore*— Fig. 10). 
for the insulated plate with all air properties evaluated 


These curves are 


at the free-stream temperature. 

The effect of evaluating the air properties at the wall 
temperature is shown in Fig. B (Van Driest*— Fig. 3) 
by the curve marked C,’V Re. The effect of heat 
transfer on the mean skin-friction coefficient for a flat 
plate is shown in Fig. B, in terms of the ratio of the 
wall temperature to free-stream temperature, 7,, 7°... 

(2) Recovery Factors for Laminar Boundary Layer. 
The various theoretical analyses for recovery {factors 
for both subsonic and supersonic flow over a flat in 
sulated plate lead to the simple result, 


gm V Pr (1) 


This approximate rule for laminar boundary layers 
agrees within | per cent of the more exact and com- 
plicated function if the Prandtl Number lies between 
1.2 and 0.7 and for Mach Numbers less than about 8. 
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FRICTION AND 


For large values of the Mach Number, however, it 1s 
not clear what temperature should be used to evaluate 
the Prandtl Number. Note that this rule predicts that 
the recovery factor is independent of Mach Number 
and Revnolds Number. 

If dissociation of the fluid at high Mach Numbers is 
taken into account, as Moore’? has shown for air, the 
Prandtl Number varies strongly with the degree of dis 
sociation, and the above rule does not hold. For this 
case, the temperature rise of the insulated plate over 
the free-stream temperature, rather than the recovery 
factor, is shown in Fig. C (Moore*—-Fig. 15 For 
sea-level conditions, it 1s seen that dissociation begins 
to play an important role at a Mach Number greater 
than Sif radiation effects are ignored. It is also evi 
dent that the effects of dissociation at these high Mach 
Numbers will be diminished if the cooling effect of 
radiation is considered at these high temperatures. 

3) Heat-Transfer Coefficient for Laminar Boundary 


Layer. For incompressible flow, the local coefficient of 
heat transfer for a flat plate is given by 
Nu, = hL/k 0.33Re, °Pi (2a 
or by 
ot = 2/66 = O.33he, rr (2b 


where the coefficient, /7, is defined by 
h (q/A)/(to — t (3 


For compressible flow in the boundary layer of a flat 
plate, the same form of theoretical result as Eqs. 
2a) and (2b) is obtained for the case of isothermal plate, 
heat-transfer coefficient is 


provided the ‘‘effective’’ 


defined by 
h, (qg/A)/(te — tes (4) 


and provided the values of the fluid properties are 
evaluated at the temperature just outside the boundary 
layer. The ‘‘constant’’ of Eqs. (2a) and (2b) changes 
with Mach Number for compressible flow. 
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Moore? has shown that the calculated value of the 
heat-transfer coefficient based on dissociation increases 
with increasing Mach Numbers and simultaneously 
decreasing absolute pressure to such an extent that these 
values can become double or triple the values of the 
corresponding heat-transfer coefficient without dis 
sociation. 

The variation of the product of Stanton Number 
and square root of length Reynolds Number versus 
Mach Number is shown in Fig. D (Van Driest® 
Fig. +) for several values of the ratio of wall to free 


stream temperature. The effect of evaluating the air 
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properties at the wall temperature is also shown by the 


curve labeled Cy’V Re. 


(B) Experimental Work on Laminar Boundary Layer 


(1) Data on Friction Coefficients for Laminar Bound- 
ary Layer.—The amount of experimental data avail- 
able on skin-friction coefficient is limited in range of 
Mach Number and of length Reynolds Number. 

The data of Blue‘ obtained by two methods (one an 
interferometer study and the other by a total-pressure 
probe), for a Mach Number of 2.0, showed that the 
measured average friction coefficients were from 7 to 
39 per cent larger than the theoretical value of Crocco 


AERONAUTICAL SCIENCES 


FEBRUARY, 1954 





over the range of Reynolds Number from 0.3 to 1.1 X 
10°. 

The data of Higgins and Pappas,* obtained by bound 
ary-layer velocity profiles, agreed well with the inde- 
pendent work of Blue, although their Mach Number 
was 2.4. The average friction coefficients of Higgins 
and Pappas were from 32 to 48 per cent larger than 
the theoretical values. 

Further work by Maydew and Pappas,® obtained by 
impact-pressure surveys of the boundary layer, showed 
average friction coefficients for the flat plate which 
were from 37 to 94 per cent larger than the theoretical 
values at a Mach Number of 2.4. 

Potter,’ at the Naval Ordnance Laboratory, found 
that average skin-friction coefficients for laminar bound 
ary layers on cones and cone-cylinder combinations 
agreed with flat-plate values with a maximum deviation 
of about 20 per cent when appropriate conversion was 
made to a flat-plate geometry. These measurements 
covered the range of Mach Number from 1.86 to 4.24. 
A typical result is shown in Fig. E (Potter’——Fig. 7 
Potter's friction coefficients are based on force meas 
urements and not on interpretations of data on ve 
locity profiles. 
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Liepmann and Dhawan,* Dhawan,’ Coles,"’ and 
Coles and Goddard'' presented data for local skin- 
friction coefficients on flat plates based on a sensitive 
force measuring element floating in the flat plate. 
The measurements covered the range up to a Mach 
Number of 4.5 and up to a Reynolds Number of 10. 
These data are probably the best local friction coeffi- 
cients available at present and agree within a few per 
cent with the theoretical values predicted for a laminar 
boundary layer over a flat plate with zero pressure gra 
dient. A typical result is shown in Fig. F (Coles” 
Figs. | and 3). 

From this group of experimental data it is evident 
that the available measurements of skin-friction co- 
efficients extend to a Mach Number of 4 and that these 
data agree with the latest theories over this range of 
Mach Number probably within the experimental error 
inherent in the particular type of measurement. 

2?) Data on Recovery Factors for Laminar Boundary 


The experimental data on recovery factors 
summarized in 


Layer. 
for laminar boundary 
Table 2. 

The local recovery factors of Wimbrow”? for one cone 
at a Mach Number of 2.0 and for one paraboloid at 
Mach Numbers of 1.5 and 2.0 agreed within | per cent 
of the theoretical value based on the square root of the 
Prandtl Number evaluated at the adiabatic wall tem- 
perature. These recovery factors were independent 
of Mach Number, Reynolds Number, and body shape 
but appeared to increase slightly with surface rough- 


layers are 


ness. 

The recovery factors of Stalder, Rubesin, and Tende- 
land'* for a flat plate were independent of Reynolds 
Number and equal to 0.881 + 0.007. This value is 4 
per cent larger than the theoretical value evaluated at 
adiabatic wall temperature. The experimental data 
are shown in Fig. G (Stalder, Rubesin, and Tende- 
land'8—Fig. 5). 

Some typical results of the recent experimental values 
of Eber’ for cones and cone-cylinders are shown in 
For the 


Fig. H (Eber!*—Figs. 3, 4, 5, 7, 8, and 9). 
laminar boundary layer, Eber’s recovery factors agreed 
with the square root of the Prandtl Number evaluated 
at the temperature within 1 per 


cent. 


adiabatic wall 


TABLE 2 
Data on Recovery Factor for Laminar Boundary Layer 





Author Date Model 

Wimbrow!” 1949 Cone 
Paraboloid 

Stalder, Rubesin, and Tendeland! 1950 Flat plate 

Eber! 1952 Cones (10°-80° 
Cone-cylinders 

des Clers and Sternberg"® 1952 Cone 

Slack!é 1952 Flat plate 


Stine and Scherrer’ 1952 Cone 
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Nominal Mach Number = 2.4 
Stagnation Pressure 
5 psia | 
10 psia | 
15 psia 
6 20 psia ] 
4 25 psio 
he 30 psia 
5 © 35 psio 
ee 40 psia | 
u. 0 45 psia 
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REYNOLDS NUMBER 
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LOCAL TEMPERATURE - RECOVERY FACTOR ON FLA 
PLATE 
FIG. G - STALDER, RUBESIN, AND TENDELANC 


The recovery factors of des Clers and Sternberg,'® of 
Slack,?® and of Stine and Scherrer” are also shown in 
Table 2. 

The M.I.T. recovery factors on a flat plate agreed 
with the cone values of 0.85 and not with the flat-plate 
values of 0.SS of references 13 and 16. 

In summary, up to a Mach Number of 3, the experi 
mental recovery factors agree within | per cent of the 
theoretical value given by the square root of the Prandtl 
Number evaluated at the adiabatic wall temperature. 
An unexplained discrepancy of about 4 per cent exists 
in some of the experiments made with flat plates. All 
evidence indicates the laminar recovery factors are 
independent of Mach Numbers up to 4 and of Reynolds 
Numbers up to the beginning of transition of the lami 
nar boundary layer. 

(3) Data on Heat-Transfer Coefficients for Laminar 
Boundary Layer.— Experimental data on heat-transfer 
coefficients for laminar boundary layers in supersonic 
flow are scarce. The following papers contain some 
results for supersonic flow. 

The local heat-transfer coefficients of Slack,'® for a 
flat plate at a Mach Number of 2.4, are shown in Fig. I 
(Slack'*—Fig. 9). The uncorrected heat-transfer data 
(not shown in Fig. 1) scatter by a factor of 5, whereas 


» 


Reynolds Mach Recovery 
Number Number Factor 
2.7 X 10 2.0 ) $55 + 0.008 
BS x iC 1.5 ) 845 + 0.008 
2 0 ) $55 + 0.008 
oe 4 x 108 94 Q S81 + 0.007 
6 x 103 0. 88—4.65 0.845 + 0 O08 
- 5 «x 11% ; - 
0.1 — 1.3 X 10 2.18 0.851 0 007 
0.15 — 3 xX 108 2.4 ) 884 + 0.006 
0.2 — 1.3 X 10! 2.0 845 
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RECOVERY FACTOR ALONG A 40°CONE - 
CYLINDER. 


FIG. H-EBER'® 
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COMPARISON OF SEVERAL ANALYTICAL METHODS 
WITH EXPERIMENTAL SKIN-FRICTION COEFFICIENTS 
ALONG A FLAT PLATE, M=2.5 


FIG. J-RUBESIN, MAYDEW, AND VARGA®°-195) 


the data shown in Fig. I (Slack'®— Fig. 9), corrected 
for variable surface temperature, scatter from —15 to 
+100 per cent relative to the theoretical line of Chap- 
man and Rubesin'* for constant wall temperature. 

The heat-transfer coefficients of Eber,’ some of 
which are shown in Fig. H, for cones with Mach Num- 
bers ranging from 0.SS to 4.2, scatter from his correla- 
tion line by about +20 per cent. His correlation for 
cones agreed well with the theoretical value for flat 
plates when the appropriate conversion was made from 
cone-type flow to flat-plate flow. 
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The local heat-transfer coefficients of Scherrer and 
Gowen" for a cone at a Mach Number of 2.0 indicated 
agreement within 10 per cent of the theoretical value 
at the base of the cone but showed a large difference of 


50 per cent at the nose. 


C) Theoretical Work on Turbulent Boundary Layer 


The problems of turbulence in general, of the turbu 
lent boundary layer in particular, and of transition 
from a laminar to a turLulent boundary layer are prob 
ably the most important problems in fluid mechanics 
whose detailed nature and mechanism are not under 
stood. Hence, all theoretical investigations on turbu 
lent boundary layers are based on arbitrary and simpli- 
fied models or sets of assumptions, which in turn per 
mit the calculation of desired quantities to a first ap 
proximation. The extraordinary feature is that the 
results based on these arbitrary models agree fairly 
well with the available experimental data in many 
cases. 

(1) Skin-Friction Coefficient for Turbulent Boundary 
Layer. Several reviews have been presented of the 
many and various models used for the turbulent 
boundary layer in supersonic flow. In order to keep 
the size of this survey within bounds, the excellent 
review given by Rubesin, Maydew, and Varga” will be 
the only one discussed here. The basis of these models 
is to use one set of properties, velocity profiles, etc., for 
the laminar sublayer and another set of properties, 
eddy viscosity, velocity profiles, etc., for the turbulent 
portion of the boundary layer. Most of these models 


depend on slightly different methods of integration of 
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the von Karman momentum integral for the turbulent unity. He developed a general formula for skin fric 
boundary layer. The resulting formulas for local and _ tion including heat transfer. One of his charts for skin D 
average skin-friction coefficients for a flat plate are not friction is shown in Fig. L (Van Driest??—Fig. 18). ( 
simple; for this reason the theoretical results are Donaldson** rederived the form of the incompressible Bo 
shown here only in chart form. turbulent skin-friction law for a plate so as to extend it ski 
The friction coefficients of Rubesin, Maydew, and to compressible flow with an arbitrary set of velocity in 
Varga,” obtained by extending the original analysis of _ profiles. His theoretical values compare very well up not 
Frankl and Voishel, are shown in Fig. J (Rubesin, toa Mach Number of 5 with those of Van Driest?? and tur 
Maydew, and Varga*’—Fig. 9) fora Mach Number of | Rubesin, Maydew, and Varga” based on more com- nec 
2.5. Wilson*! extended von Karman’s analysis for in- _ plicated analyses. of 
compressible flow to include the effects of compressibil- (2) Recovery Factors for Turbulent Boundary Layer. the 
ity. His results are shown in Fig. K (Wilson*'—Fig. The theoretical work on recovery factors is limited to per 

6). Van Driest*? derived the continuity, momentum, the assumption of constant properties in the turbulent 

and energy equations for turbulent flow of a com- portion of the boundary layer—i.e., to a constant value 

pressible fluid and used an eddy Prandtl Number of of eddy viscosity and eddy conductivity. The early 

TABLE 3 
Data on Skin Friction for Turbulent Boundary Layer 
Reynolds Mach Wit 
Author Date Model Number Number Method 

Wilson?! 1950 Flat plate —19 X 10° 2.2 Velocity profiles Sta 
Rubesin, Maydew, and Varga?® 1951 Flat plate —6 X 10! 2.5 Velocity profiles Hil 
Potter’ 1952 Cones, cone-cylinders —S xX 16 1. 86-4 .24 Drag forces Eh 
Bloom? 1952 Flat plate —10 X 16 5.5 Velocity profiles : 
Brinich and Diaconis*® 1952 Hollow cylinder —-H x i¢ 3.05 Velocity profiles, Schlieren ob cs 
servations Sla 
Coles" 1952 Flat plate —10 X 108 1.5 Floating element Stir 


Coles and Goddard!!! 1952 Flat plate —10 X 10 4.5 Floating element 










































work of Ackermann, using a kinetic-theory model, 107 * or TT 3 
showed the recovery factor to Le equal to the cube root E e 
of the laminar Prandtl Number for Prandtl Numbers - Equations 66 and 45 4 
greater than 0.5 and less than 2. Other analyses gave a 
similar results except that a decrease of recovery factor 
with increasing Mach Number was indicated. = 
The work of Tucker and Maslen** extended the in- - 
compressible analysis of Squire for a flat insulated a Ma = 0 
plate to compressible flow. Their results can be rep ¥ 
resented by the following approximate rule within | 9 toe 
per cent = os F Equations 66 and 45 
~ + 1+ 0.5281, ‘| — = 
mr = - lan Pr z S 
3N + 1+ M.,’ ws 
rT) 
where Pr is the laminar Prandtl Number evaluated at " 1073 
the adiabatic wall temperature, .1/. is the free-stream ) 
Mach Number, and _N is the reciprocal of the exponent o 
of the turbulent boundary-layer velocity profile ap- « 
proximated by the power law. This relation holds for -" oi 
Prandtl Numbers greater than 0.65 and less than 0.75. - |e J 
This relation reduces to the cube-root rule of Acker- = a Equations 66 and 45 = 
mann for zero Mach Number and large value of .V or ry — Tw 1.0 Tw. 3.0 = 
9 large Reynolds Numbers. < a “Teo” Too \ 4 
(3) Heat-Transfer Coefficients for Turbulent Boundary = s 
Layer.—The theoretical work on heat-transfer coeffi- ad a 
cients for a turbulent boundary layer is limited to the = 
analysis by Van Driest*? and to the use of the Reynolds = 
analogy for the work on skin-friction coefficients Ly 3 a q 
Rubesin, Maydew, and Varga.” The results of Van = 
Driest’s analysis are shown in a series of charts in Fig. a ' 
M (Van Driest?*—Figs. 1, 2, 3, and 5) covering a range a z Equations 66 and 45 : 
of Mach Numbers from 0 to 6 and several values of the E iw =2.0 Tw. 6.0 ; 
ratio of wall temperature to free-stream temperature. — f Too Too : 
ric. y [ 4:40 M:82 ] 
skin D) Experimental Work on Turbulent Boundary Layer 10-3 | /Te Too 
). (1) Data on Skin-Friction Coefficient for Turbulent 
ible Boundary Layer.—Most of the experimental work on - 
d it skin -friction coefficients for a turbulent boundary layer W a cities 
“ity in supersonic flow is summarized in Table 5. Since ee 
up none of the experimental models had a completely oailite =e : mana ae A 
ind turbulent boundary layer from the leading edge, it was 10 10 10 10 
mii- necessary to make a correction for the laminar portion REYNOLDS NUMBER, Roy 
of the total boundary layer. This correction altered TURBULENT HEAT-TRANFER COEFFI- 
the measured value of the local coefficient by 10 to 500 CIENT FOR AIR 
to per cent. FIG. M-VAN ORIEST?*-195! 
ent 
lue 
ty TABLE 4 
Data on Recovery Factors for Turbulent Boundary Layer 
Reynolds Mach Recovery 
Author Date Model Number Number Factor 
Wimbrow!” 1949 Cone —2.7 X 10 2.0 0.885 + 0.008 
Paraboloid —4.8 X 105 1.5 0.902 + 0.005 
2 + 
Stalder, Rubesin, and Tendeland 1950 Flat plate —7 X 10 2 y ° an - a 997 + ) (WT 
Hilton?’ 1951 Flat plate —10 X 10 20 0. 880 + 0 004 
Eber?* 1952 Cone and cone-cylinder —1 X 108 2.87 0.92 
—0.25 X 108 25 0.97 
»b des Clers and Sternberg’ 1952 Cone and cone-cylinder —7 X 108 9-3 .4 0. 882 + 0.007 
Slack '!* 1952 Flat plate —3 X 10% 2.4 0.906 
Stine and Scherrer’? 1952 10° Cone 0.4 —4 X 10 1 .97-3.77 0.882 + 0.008 
40° Cone-cylinder 0.3 —1 X 108 3.10-3.77 0.885 + 0.011 
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With the exception of the work of Potter’ on cones, 
the experiments of the other investigators in Table 3, 
showed consistent values of turbulent skin-friction 
coefficients which agreed within about 10 per cent of 
the theories of Wilson,?! of Rubesin, Maydew, and 
Varga,” and of Van Driest.*? In order to test these 
theories to better than an error of 10 per cent, it ap- 
pears necessary to measure turbulent skin-friction co- 
efficients to a much smaller error. If, however, the 
variable effects of transition are considered, this re- 
quirement appears to be a difficult condition to achieve 
in tests. 


(2) Data on Recovery Factors for Turbulent Boundary 
Layer.—It is evident from the data in Table 4 that, with 
the exception of Eber’s data on cones and cone-cylin- 
ders and possibly Slack’s data on a flat plate, the re- 
covery factors for a fully turbulent boundary layer 
agree to within 1 per cent with the theoretical approxi- 
mate rule—i.e., equal to the cube root of the Prandtl 
Number. Furthermore, the data indicate that these 
recovery factors are independent of Reynolds Number 
and Mach Number over the ranges tested. The theo- 
retical work of Tucker and Maslen*‘ does not agree with 
these experimental results. 
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(3) Data on Heat-Transfer Coefficients for Turbulent 
Boundary Layer.—Experimental data on heat-transfer 
coefficients for a turbulent boundary layer in super- 
sonic flow are almost nonexistent. Eber'* has meas- 
ured a few heat-transfer coefficients on a cone at a 
Mach Number of 2.87 and over a limited range of 
Reynolds Numbers from 0.4 to 0.9 X 10°; he found 
his values of the Nusselt Number to scatter from + 10 
to —40 per cent relative to the well-known Colburn 
equation for subsonic flow in a turbulent boundary 
layer—.e., 

Nu = 0.029Re® Pr** (5 

Slack'® has also measured some heat-transfer coeffi- 
cients for a flat plate at a Mach Number of 2.4 and over 
a Reynolds Number from 1.8 to 3 X 10° The 20 odd 
points are shown in Fig. I (Slack'® Fig. 9). 


M.I.T.-ONR PROGRAM 


A research program, sponsored by the Office of 
Naval Research, has been under way for several years 
in the Department of Mechanical Engineering, with 
the objective of obtaining reliable data on the rate of 
heat transfer to air moving at supersonic speeds in a 
round tube. This program consists of two separate 
parts. In the first part, a well-insulated apparatus 
was used to measure the values of the local adiabatic 
wall temperature and local static pressure of a super- 
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sonic stream of air in a tube with an entrance Mach 
Number of 2.8. In the second part, an apparatus 
using steam condensing outside a round brass tube was 
used to measure local coefficients of heat transfer to a 
similar supersonic stream of air from the same nozzle. 

The results of the first part of this program are avail- 
able in references 31, 32, and 33; hence, these data will 
not be discussed in detail. The results of some pre- 


measurements of heat-transfer coefficients 


In addition, preliminary results will be 


liminary 
will be given. 
given of some theoretical work for the compressible 
laminar boundary layer in the entrance region of a 
round tube. 

(1) Adiabatic Flow in Tube. 
sent some typical results for supersonic flow of air in a 


Figs. N, O, and P pre- 


round tube under adiabatic conditions for the case 
where a laminar boundary layer appears to exist in the 
tube. The calculated quantities are determined from 
two flow models: 


one-dimensional flow model based on constant prop- 


the first corresponding to the usual 


erties at any cross section, and the second based on a 
two-dimensional flow model where the properties vary 
over the cross section of the tube. In addition, the 
preliminary calculations obtained from the M.I.T. 
Differential Analyzer in the solution of the theoretical 
problem, from the Doctoral Thesis of Toong,** are also 
shown. 

The local coefficients of friction based on the two- 
dimensional flow model are in fair agreement with the 
values for a flat plate and the theoretical values for a 
tube, where it appears that a laminar boundary layer 
exists in the tube. In Fig. P, this agreement is excel- 
lent. 

The recovery factors based on the two-dimensional 
flow model vary from 0.87 to 0.89 for that portion of 
the tube flow where it appears that a laminar boundary 
layer exists. These values are 2-3 per cent larger than 
the theoretical values for a laminar boundary layer on a 
flat plate with zero pressure gradient. 

The calculated thickness of the laminar boundary 
layer in the tube, based on the two-dimensional flow 
model, is also shown in Figs. N, O, and P. The agree- 
ment between this quantity and the thickness deter- 
mined from Howarth’s flat-plate formula and from the 
Differential Analyzer solution is excellent for this type 
of calculation. 

(2) Heat-Transfer Measurements in the Tube. 
Preliminary heat-transfer data are shown in Fig. Q for 
the case where a laminar boundary layer appears to 
exist in the supersonic flow in the tube. The scatter 
present in these data is typical of most measurements 
of local heat-transfer coefficients. The data agree 
well with the theoretical value for a flat plate with zero 
pressure gradient and with the Differential Analyzer 


solution for tube flow. 


M.I.T.-OAR PROGRAM 


A research program, sponsored by the Office of Air 
Research, was started about | year ago as a joint effort 
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FIG. P 


of the Departments of Aeronautical Engineering and 
Mechanical Engineering with the objective of obtaining 
reliable data on heat-transfer coefficients to air moving 
at supersonic speeds over a flat plate. The experi 
mental program is being carried on at the M.I.T. 
Naval Supersonic Wind Tunnel over a range of Mach 
Numbers from 1.9 to 3.1 and up to a length Reynolds 
Number of 17 X 10°. 

This program has been divided in two parts. The 
first part consists of an attempt to secure a flat-plate 
model that corresponds as nearly as possible to an in 
sulated flat plate. On this model the local recovery 
factors and velocity profiles in the boundary layer could 
Thus the type of boundary layer pres- 
The second part 


be measured. 
ent on the plate could be identified. 
of the program will be concerned with construction of 
a separate flat-plate model made of brass to secure 
local heat-transfer coefficients by use of condensing 
steam. Velocity profiles will be measured to deter- 
mine the type of boundary layer present on the heated 
plate. The effects of turbulence promoters and of 
large temperature differences between the free stream 
and wall will be investigated. 

To date the insulated plate tests have indicated that 
the recovery factor for a laminar boundary layer, 
properly corrected for the effect of the heated nose 
resulting from the shock at the nose is 0.850 + 0.005 
for Mach Numbers from 1.9 to 3.1. This value agrees 
within 1 per cent with the theoretical value based on 
the square root of the Prandtl Number evaluated at 
the adiabatic wall temperature and also agrees with 


similar data on cones. 
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STANTON NUMBER VS. LENGTH REYNOLDS NUMBER 


The recovery factor for a turbulent boundary layer 
is 0O.SSI_ + 0.008 for Mach Numbers of 1.9 to 3.1 and 
up to a Reynolds Number of 17 X 10®. Within ' 2 per 
cent, this recovery factor is independent of Mach Num- 
ber and of Reynolds Number over the ranges measured. 


CONCLUSIONS 


In 1948, the author examined the problem of the 
transient temperature distribution in a wing moving at 
supersonic speeds; he found it necessary to review the 
state of the art relative to heat transfer at supersonic 
speeds. One conclusion reached at that time was that 
the state of the art was relatively undeveloped and 
that a considerable amount of theoretical and experi- 
mental work was required to improve it significantly. 

In the past 4 years significant contributions have been 
made to both the theoretical and experimental work 
on laminar boundary layers up to about a Mach Num- 
ber of 38. Small discrepancies between theory and ex- 
periment still exist in this range for friction coefficients 
Additional experimental data 
The 
greatest need is to extend the experimental data beyond 
a Mach Number of 3 and to investigate problems such 
as the effect of a pressure gradient, the effect of radi- 
ation, the effect of transition, the properties of air at 


and recovery factors. 
are greatly needed for heat-transfer coefficients. 


temperatures of the order 3,000°F, ete. 


F 
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In the past 4+ years several significant contributions 
have been made to the theory of the turbulent boundary 
layer for supersonic flow; these have been based on 
arbitrary models or assumptions because of the lack 
of knowledge of the details of turbulence. The experi 
mental data on skin-friction coefficients and recovery 
factors for turbulent boundary layers have been ex- 
tended to a Mach Number of about 3 and a length 
Reynolds Number of about 20 X 10°. It is curious 
that for this case where the basis of the theory of the 
turbulent boundary layer is fairly uncertain, the agree- 
nent between the predicted and measured values of 
the skin-friction coefficient and recovery factor is the 
best of all available comparisons at the moment. There 
is a great need for reliable data on heat-transfer co- 
efficients for turbulent boundary layers, since, prac- 
tically speaking, none is available. The experimental 
work for the turbulent boundary layer should be ex- 
tended to higher Mach Numbers and to larger Reynolds 
Numbers; the effects of pressure gradient, radiation, 


etc., should be investigated. 
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Note on an Interaction Between the Boundary 
Layer and the I nviscid Flow 


Antonio Ferri and Paul A. Libby 
Department of Aeronautical Engineering and Applied Mechanics 
Polytechnic Institute of Brooklyn, Brooklyn, N.Y 


November 2, 1953 


"acm TO THE CLASSICAL boundary-layer theory the flow 
about bodies at Reynolds Numbers of aeronautical interest 
can be considered as composed of two régimes: an outside in 
viscid flow and a thin boundary-layer region adjacent to the 
body. This point of view leads to the approximation that, on a 
slightly curved surface, 0p/Oy throughout the layer is negligibly 
small. The additional assumption that the inviscid flow is irro 
tational leads to the requirement that Ou/Oy is zero at the outer 
edge of the boundary layer. In this theory any interaction be 
tween the two régimes is accountable by a simple correction to 
the body shape based on the boundary-layer displacement thick 


ness 


Recently, in connection with hypersonic laminar boundary 
layers, this classical point of view has been modified; an inter- 
action between the two flow régimes leading to a self-induced 
axial pressure gradient has been considered. It is the purpose 
of the present note to point out another type of interaction which 
may be of practical importance and of fundamental interest even 
at Mach Numbers below those considered in the hypersonic 
boundary-layer theory and which may have to be considered 
in that theory 


This interaction may be discussed from the following point of 
view: It is evident that the division of the flow field into the 
aforementioned two régimes is somewhat arbitrary but well 
defined for low-speed flow. However, in what would correspond 
to the inviscid régime for supersonic flow, it is possible to gen 
erate, by means of upstream highly curved shock waves, a flow 
region adjacent to the surface with large entropy gradients and 
large vorticity. In such regions the pressure and velocity are 
not uniquely related, and thus for slightly curved surfaces 
Op/oy can still be small, although Ou/Oy can be on the same 
order of magnitude as exists in classical laminar boundary layers. 
This possibility seems to require a modification of the classical 
approach when highly curved shock waves exist at the leading 
edge 

In order to indicate the practical importance of this inter- 
action, consider two slender bodies of revolution. The first is 
sharp-nosed and has the usual almost conical shock wave. The 
second has the identical shape except for a rounded nose that 
fairs rapidly into the same contour as the first body. It is pointed 
out that the ogive of minimum pressure drag is just such a 
rounded nose body. At the nose of this second body will be a 


highly curved shock that generates a layer of large entropy gra 


dients and high vorticity close to the surface of the body. Soi 
what downstream of the nose of this second body the pressur« 
as determined by inviscid flow theory, will be essentially uniform 
through this vorticity layer and will be equal to that at the cor 
responding point of the first body. However, the velocity 
the wall will be less in the second body 

If, now, boundary-layer calculations are made for the two 
bodies, it is clear that this vorticity layer has two effects Phe 
“outer edge’’ of the boundary layer in the classical 
Both alterations 


velocity at the 
sense is reduced, and Ow/Oy there is not zero. 
are functions of the local displacement thickness and of the vor- 
ticity distribution, are interdependent, and tend to reduce the 
skin friction and heat transfer if the boundary layer remains 
laminar. It is not evident what the influence of the outside 
vorticity on the transition would be. If the boundary layer is 
turbulent on both bodies, the outside vorticity would probably 
have a marked influence on the turbulent skin friction and heat 
transfer 

The fundamental significance of this type of interaction be 
tween the essentially inviscid and the boundary-layer flows ap 
pears to be that the division of the flow into separate régimes is 
not well defined for this case. The boundary between the re 
gions in which shear stresses are and are not important is not 
well defined. 

To develop some feeling for this interaction, the following 
analysis was carried out: An incompressible, constant pres- 
sure, laminar boundary layer was superposed on an inviscid vor 
ticity layer of a given thickness. The velocity in the vortex 
layer was assumed to vary from a value at the wall of 90 per cent 
of free stream to the free-stream value within the thickness of the 
vorticity layer. The boundary layer was required to be thin com 
pared with the vorticity layer so that displacement effects were 
small. Consequently, the values of Ou/Oy at the outer edge 
were small although included in the analysis. As in the usual 
integral approach, the boundary-layer equation was integrated 
from the wall to a boundary-layer thickness 6 = 6(x). Within 
this layer, the viscous stresses were thus taken to be significant 
At y = 6, the velocity and velocity gradient 0u/Oy were matched 
to the corresponding values in the inviscid flow at y = 6. At 
The skin 


friction was thus found as a function of Reynolds Numbers based 


the wall, the usual no-slip conditions were imposed 


on vorticity-layer thickness and on distance from the leading 
edge and of the per cent deficiency in the inviscid wall velocity 
The results of this simple analysis indicate a reduction in skin 
friction of about 15 per cent. It is pointed out that this defi 
ciency is on the order of what can be expected in high supersoni 
flow. These results appear to be insensitive to the thickness 
of the vorticity layer probably because the boundary-layer thick 
ness was required to be small compared to the same 

These considerations seem to indicate that the influence of up 
stream highly curved shock waves on boundary-layer character 
istics warrants further investigation 
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Aerodynamic Flutter Coefficients for an Airfoil 
with Leading- and Trailing-Edge Flaps in 
Two-Dimensional Supersonic Flow 


Hermann J. Hassig 
Dynamics Engineer, 
Calif 


November 9, 1953 


Lockheed Aircraft Corporation, Burbank, 


S' PERSONIC TWO-DIMENSIONAL AERODYNAMIC FLUTTER CUEFFI 
CIENTS for a wing with leading- and trailing-edge flaps ar 
presented in terms of the coefficients for a wing without flaps 
With the help of the results in reference 1, the coefficients, de 
scribing the aerodynamic forces and moments on a harmonically 
oscillating wing with leading- and trailing-edge flap in a two 
dimensional supersonic flow, can be expressed in terms of quan 
tities tabulated in references 2 and 8 
used for a strip-type flutter analysis of a wing with arbitrary ratios 
of leading- and trailing-edge flap chord to wing chord, numerous 
Because inter 


If these coefficients are 


rather inaccurate interpolations are required 
polations for hinge axes position and for the value of & are neces 
sary, they are not amenable to high-speed digital computers 

Che purpose of the present note is to express all the required 
coefficients in terms of the coefficients for the harmonic oscillation 
of a wing without flaps, thus avoiding the coefficients that depend 
on the location of the hinge axes as well as their interpolation 
An indication that this should be possible for a wing-aileron com 
bination is found in reference 2, where it is stated that fo(.V, &) is 
the only quantity that needs evaluation. It can be shown from 
reference 2 that Lie, Msc, Nie, Ns, and Nse,* all depending on the 
location of the hinge axis, can be expressed in terms of Ly, Lys’, 
V..’, and Mas’, which determine the force and moment due 
With the help of refer 


” 


to h and @ motion for x» = O (Fig. 1 
ence 1, this can be extended to include the leading-edge flap 

However, the required coefficients can be obtained in a much 
simpler way by taking advantage of the characteristic of super 
sonic flow that a downstream element cannot have any effect on 
in upstream element and by using a simple technique of super 
position of displacements and loads 

When P, Ma, Mg, and M , are, respectively, the complex am- 
plitude of the downward force, the tail-down moment about a 
reference axis, the tail-down hinge moment of the trailing-edge 
flap and the nose-up hinge moment of the leading-edge flap due to 


* Lee Ls + ils, etc.; Ls, Le, etc., are defined in reference 2. 
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the complex amplitudes 4, a, 8, and y of the displacements (Fig 


1), the following matrix equations can be written 


tig An, Ane Ang Ar | i ) 
M.. secne | ten Awe Aap Ae a 
Ve ws Ag, 1 ser 1 3g ly (° 
r | ce wee ee we ; | 


where p is the air density, v is the forward speed, k& is the reduced 
frequency wh/v, and 6 is the semichord 
In order to demonstrate the method two examples ure given 


EXAMPLE (1) DERIVATION OF Aygg 


Ihe force on the airfoil due to trailing-edge flap deflection is 
equal to the force on the flap itself, since a downstream element 
cannot have any effect on an upstream element. From refer 


ence 2, the expression for the force due to a is 
P)a = 4 pbv*k* Lg (R ja 2 


Replacing 6 by (1 — x,)b, the semichord of the flap, and a 
by 8, the flap deflection, leads to the following expression for the 


force on the flap due to 8 
P)g = —4p(1 — x, )bv*(1 v; )*R?L 3 (1 x kip 3 


It will be noted that L.,’ is now evaluated at (1 — x))k rather 


thanatk 
Comparison of Eq. (3) with its equivalent term in Eq. (1) leads 


to 


tnal(k) = (1 — m)8La’ [1 — xk ' 


EXAMPLE (2) DERIVATION OF Ay, 


Leading-edge flap deflection y can be obtained from super 
position of the three displacements shown in Fig. 2 
The force on the entire section is determined by the following 


coefficients 


Due to I Lsa'(R } 
Due to II —(1 — x2)®L,,'(01 xe )k 5 
Due to III: —2x:Ay,(k \ 


It follows that 


k) = Lu’(k) — 


1), (1 — x)*Ly’(Q1 x2)R| — 2x2Ap,(k) (6 


Using the technique shown in the above examples, the re 


maining coefficients in Eq. (1) can be determined. The results 


ure as follows 
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ILLUSTRATION OF THE DERIVATION OF An yl) 





Ann(R) = Lirlk) 

Analk) = Las'(R) — 2xoAny(k) 

Angl(k) = (1 — x1 )8Lag’ (1 — xR] 

Apny(R) = Las'(R) — (1 — x2)8 Lay’ [ — xk] — 
2x2A pn(R) 

Lan(R) = My'(R) — 2xoAnn(k) 

Aaalk) = Ma,'(R) — 2xot Las'(R) + Aan(k)} 

1 

1 





aplk) = Agelk) + 2x1 — xo )Anglk) 

ay(Rk) = Mas'(k) — 2xolss'(kR) — (1 — xe) (1 — 
x )R| + 2(x9 — x2) (1 — xe)3Ls,’(1 — 
x )R] — 2x2A ay(R) 

12,(k) = My'(k) — 2eiLwl(k) — xFtMie'(aik) — 


2L12(x1k)} 

Agalk) = Mss'(k) — 2xiLas'(k) — x14} May'(mik) — 
2134'(x1k)} — 2xoA Balk) 

Aga(k) = (1 - x1 )*Ms4’[C1 — x1 )R] 

1g4(k) = Agalk) + 2A py(k) — (L — x2)*Mss[ — 
x )k| — 2x2A g,(R) + 2x, — x2) (1 - 
x2)3L34" [(1 — x2)R] + (x1 — x2)41 Ms," [(rr — 
x )k| — 2Lss'[(x1 — x2 )RI} 

Ay,(R) = x21 Myo'(x2k) — 2Lo(x2k)} 

AyalRk) = X24t Mgy'(xok) — 214’ (x2k)} — WoA yp(k) 

A,g(k) = 0 

Ay+(k) = 2Ax — X»)A +,(R) + A valk) 








The coefficients Ly», Lsy’, My’ and V3,’ in the above equations 
are tabulated in reference 2 for several Mach Numbers and for 
a large range of the argument k. The lowest k value listed is of 
the order of k = 0.01. If a 10 per cent flap chord is used, the 
lowest k value, based on the semichord b for which the A;; co- 
efficients can be obtained from these tables, is of the order k = 
0.1. Thus if lower values of k are needed, the tables have to be 
extended to lower values of the argument & if Eqs. (7) are used 


Although the presented method involves interpolation of Ly, 
Ly’, Mi’ and M3,’ for nontabulated values of the arguments k, 
it seems, however, that the amount of work involved is much less 
than when interpolations for hinge axes location, as well as re- 
duced frequency k, are necessary, as is the case when the co- 
efficients in reference 3 are used. This is especially true if Eqs. 
(7) are used in combination with high-speed digital computers 
that can do the interpolation. In addition, more significant fig- 


ures can be obtained. In view of these advantages it appears 
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that the tables for the plain wing given in reference 2 should be 
extended to very low values of reduced frequency & in prefer 
ence to extending the tables in reference 3 for additional hinge 
axis locations and reduced frequencies 

It is also worth noting that the same principles can also be 
applied for determining the associated indicial functions for a 
wing with leading- and trailing-edge flaps in terms of the indicia] 
functions for a plain wing. The required indicial functions for a 
plain wing have been evaluated in closed form in reference 4 
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On a Differential Equation Occurring in the 
Theory of Heat Flow in Boundary Layers with 
Hartree’s Velocity Profiles 


Henry E. Fettis 


Aeronautical Research Laboratory 
Base, Dayton, Ohio 


October 30, 1953 


Wright-Patterson Air Force 


I A RECENT NOTE IN THIS JOURNAL, Schuh! derives the equa 
tion 
d*6 £2 do 
a = 66 = 0 ] 
dg? 2 dé 
for the temperature variation in the boundary-layer flow about 
a body. In this equation, @ represents the variation of temper 
ature normal to the body when the variation along the body con- 
tour is assumed proportional to a power of the distance. The 


variable £ is defined as 


£ = g!/3p1/8(2 — B) (U;/vx)!/4y 
with 
o0 = Prandtl Number 
U = speed of laminar flow = x” 
8B = 2m/2m+1 
vy = kinematic viscosity 
e = m2 — B) 
x = distance parallel to wall 
y = distance normal to wall 
b = parameter defined by Schuh;! (proportional to velocity 


gradient at wall) 
The purpose of this note is to present the solution of Eq. (1 
in closed form, subject to the conditions 
ao) = 1, Oo) = 0 (2 
It is first convenient to change the independent variables from 


Etot, where f = £/6. The new equation becomes 


d*0 ( 2 ) doe 
+ Es + —-—¢=0 3) 
dt? 3t] dt 3l 


Introducing a new dependent variable defined by 


g = 1l/Bet/% (4 
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Eq. (3) is transformed to 
dy 1 l+2e 2 : 
+ © 4 a ) 
dt? rt» @ Qt? 


which is the standard form for the Whittaker confluent Hyper 


geometric Equation (reference 2, Chapter VI, Paragraph 2), 


d’y l k (1/4 m*| , 
a” © ae as f? J — " 
Comparison of (5) and (6) indicates that 
m = 1/6, k= (1 + 2e)/3 


A solution to Eq. (6) which satisfies the conditions at infinity is 


the Whittaker function W’;, »(t 


¢ AWs, m(t 7 


Thus 


where A is a multiplicative constant. From Eq. (4) it then fol 


lows that 


12 W_ (1 426)/3,1/6 (€3/6 (8 


6 = 6)/34¢ -e1 -€ 


To determine the constant -1, it is convenient to resort to the 


integral representation of Hy, ..(2) (loc. cit 
We, m(z) = 
1/2 ” F 
=} | 4 ” - ds 
lin ke+(1/2 Py 
9 
so that Eq. (8) may be written 
6S Ae~ e/6 ' é /3 
6 =e 9 — 2/3)(l+e d 
r((3 + 26€)/3 a 6 
10) 
For & = 0, this reduces to 
6! 44 se 
0) = sary d 1] 
I 3 + Ye 3 JV 
“13 ’ > 9 > 
6°" ATC1/3)/P[(3 Ze/3 
Hence, 
= 1(3 + 2e)/3]/6'/4r(1/3 12 
ri(3 + 2€/3 
4= t g/12 yy 26) /3. 1/6 (£3/6 13 
M1/3 
or 
rf s 6 : 
ee 2e 3 (23/6) + 5 2/3)(1+e) , de 118 
(1/3) Jo 
Eq. (14) is suitable for numerical calculation when € is small 
For large values of &, the change of variable 
= (£3/6)u 
gives 
6} 3 te £46 : 
0 = ‘ 2 ni 1+ u 2/3)(I+e é s/o ‘du 
(1/3 0 
15 
lhe temperature gradient at £ = 0 is best evaluated by expand- 
ing HW,, »{z) in the vicinity of z = 0. We obtain (see reference 
2, Chapter VI—2 
r( —2m 
HW s) = Vi.m(z) 4 
rj{(i/2 k—n 
I(2n : : 
Vi, (s) = ¢—3/22!/2 » 
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and with k (1 + 2€)/3, m = 1/6, and £3/6, this gives 


s 


We Gee 


so that 
0 Pil + (2€/3)jI 1/3 
. 18 
dtjg-9 6)/3r[(2 +2) /3)P(1/3 
The resulting value of the heat transfer coefficient is 
9 
" = > 
r(i +: ji 1/3 
} abso 4 3 
19 


The comparison between this value and that given by Schuh de 


pends on the closeness of the approximation 


which may be shown to be asymptotically correct for large values 


of «€ 
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On a Result of Nikolskii and Taganov 
Concerning Transonic Flow 


Marvin Stern 
Republic Aviation Corporation, Farmingdale, L.!., N.Y 
November 24, 1953 


B' ELEMENTARY CONSIDERATIONS We can strengthen the basi 
Nikolskii and 
their important result It 


theorem derived in the paper of Taganoy 
and thereby immediately conclude 
is impossible to have potential transonic flow if a straight portion 


of the physical contour lies in the local supersonic region.”’ 


We strengthen their basic theorem to read In potential 
transonic flow, moving along the sonic line in a clockwise direc 
tion, the angle that the velocity vector makes with the horizontal 
steadily decreases.’’ This differs from the original by omitting 


the equality sign in their statement of the theorem—-viz., 6; < 0 


We consider potential transonic flow over a convex contour 
It is assumed that both the contour, A, and the sonic line, ¢ 


c*, are continuous curves that intersect at two points. The 
sonic line then separates the flow field into two regions, g < c*« and 
q > cx (Fig. 1 

Consider a point P within the supersonic region. Through 
this point there pass two characteristics making an angle with 
A and B, 


the velocity magnitudes are equal, while the velocity angles must 


each other. At their intersections with the sonic line, 


satisfy one of the three possibilities 


O04 < Op l 
64 = On 2 
64 > On 3 
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A _B q« C. 
‘ q=C. 


I-III q>C\ 
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jp 
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Fie. 2 
I 
p' 
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‘i Pi. = I' 
Fie. 2. 


Ihe three possible image systems on the hodograph plane are 
shown in Fig. 2. The directions of J’ and JI’ are set by their 
orthogonality with // and J, respectively.” 

Without any assumption on the functional determinant, j7 = 
0(q, 8)/O(x, y), we will be able to exclude two of the above three 
possibilities. We can state that characteristics J and J] must 
map in some manner onto at least portions of the epicycloid 
characteristics I’ and JI’ on the hodograph plane. Further, 
merely the assumption of potential transonic flow guarantees 
that the point P of the physical plane maps onto one point P’ 
on the hodograph, this point necessarily being on the intersection 
of the I’ and JI’ characteristics. 

Case 1.—We assume only that q and @ are continuous. For 
point P on the sonic line, the angle between the characteristics 
is zero. It follows that, as we bring point P close to the sonic line, 
the angle between the characteristics becomes arbitrarily small 

As the angle between the characteristics becomes small, the 
intersection points A and B come arbitrarily close together 

From the continuity of 6, as the distance between A and B be 
comes small, the difference between 64 and 6g becomes arbitrarily 
small, and consequently the distance between the images A’ 
and B’ becomes sufficiently small. 

But as A’ and B’ approach each other, the intersection point 
P’ must approach the limit circle q This is a contradiction and 
Case 1 must be excluded. 

Case 2.—From the exclusion of Case 1, it follows that Case 2 
cannot allow for 6 to take on the same value at the end points of 
the finite interval AB without being constant throughout this 
interval. The validity of this statement is evident, since, 
if 6 were not constant throughout an interval, it would be de 


1954 


creasing strictly monotonically and could not therefore take on 
the same value at isolated points 

We will now show that the assumption of 6 being constant aloiug 
i finite segment of the sonic line dlso leads to a contradiction 

It is evident that the image point P’, which falls on the inter- 
section of epicycloids J’ and JI’, would in this case map coinci 
dentally with A’ and B’ on the sonic circle. Point P then has 
sonic velocity 

Similarly, the two characteristics drawn from any other 
point within the characteristic triangle APB will intersect the 
sonic line between A and BA (the segment along which @ is as 
sumed constant). It follows that all points within the character 
istic triangle APB have sonic velocity 

This extension of a sonic region into the assumed supersonic 
domain represents a contradiction, and Case 2 must therefore be 
excluded. 

It follows that Case 3 is the only possibility, and therefore the 
theorem has been proved 

It now remains to show that: ‘‘It is impossible to have po 
tential transonic flow if a straight portion of the physical contour 
lies in the local supersonic region.”’ 

Consider a point P on a straight segment of contour in a locally 
supersonic region. The two characteristics emanating from P 


terminate on the sonic curve Their images on the hodograph 


plane are as shown in Fig. 3 
On either epicycloid, the magnitude of the velocity at point ? 


can be represented as an increasing function of its argument 


qe = f(0; 6; 
dp = 1( A, . Ap 


* 


If ? were to move along this flat portion of the physical con 
tour, then @p would remain constant But if 6; < O, then 6p,, 
and 6p,, would both decrease. The first equation above would 
show qp to be increasing, while the second shows gp decreasing 
Here is the contradiction, and therefore the theorem proved on 
strict monotonicity of angle variation cannot hold. It follows 
that potential flow for this case is impossible 
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Leading-Edge Separation from Delta Wings 


R. H. Edwards 

Hughes Aircraft Company, Research and Development 
Laboratories, Culver City, Calif 

November 10, 1953 


R ECENTLY, LEGENDRE! AND ADAMS? have presented approxi- 
mate analyses of the problem of leading-edge separation 
from delta wings at supersonic speeds. In both cases the slender 
body or low aspect ratio approximation is made, the flow is con- 
sidered to be conical, a Kutta condition is applied at the leading 
edge, and the vortex sheet that separates from the leading edge 
is replaced by a concentrated vortex 

Since the flow is conical, the strength of this vortex is propor- 
tional to the downstream distance from the vertex of the wing, 
and a multivaluedness of the pressure in the fluid results. This 
multivaluedness can be removed only by the introduction of a 
branch cut that connects the separated vortex to the leading edge 
of the wing. Along this branch cut, however, a discontinuity in 
pressure appears. In effect, then, the branch cut represents a 
sheet of bound vorticity which feeds the concentrated vortex, and 
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FIGURE | 


the true vortex sheet has been replaced by this fictitious arrange 
ment in which the pressure continuity condition has been relaxed. 

Some selection principle is then necessary in order that the 
strength and position of the concentrated vortices may be deter 
mined. The Kutta condition yields one of the determining 
equations. A further assumption used by both Legendre and 
Adams is that the concentrated vortex should move with the 
fluid. Since the configuration is conical, this assumption pre 
scribes the local mean velocity at the vortex, and the problem is 
then completely defined. This procedure leads to an ambiguity 
in the determination of the lift.? 

A different assumption is suggested here 
net force on the true vortex sheet, the first mean value condition 


Since there is no 


on the simplified representation would seem to be that there 
should be no net force on the vortex and branch cut. In contrast, 
the condition of Legendre and Adams is equivalent to the stipu 
lation that there should be no net force on the vortex alone. One 
idvantage of the method suggested here is that the ambiguity 
in the determination of the lift is removed 

In order to formulate the foregoing assumptions analytically, 
let the xyf coordinate system be introduced as in Fig. 1, and 
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let us suppose that Fig. 2 represents the plane / ] Let 
x+y 

The condition of Adams and Legendre is expressible in the 
form: V(zo) = Uso. Here V is the mean velocity vector in 


the cross-flow planet = latz == 
The condition proposed here may be formulated readily. The 


force on the branch cut is given by plT(z, — ajt. The force 
on the vortex has the form —p|Viz Uz Ti The re 
quired condition is, then, 

V(zy) = U( 2s a 


A different way of looking at this assumption is the following 

At each cross section, the cross flow due to a vortex sheet has 
been replaced by the cross flow due to a single vortex. Asymp 
totically, such a replacement is achieved with a single vortex of 
strength equal to the total vorticity, located at the center of 
gravity of the vorticity of the sheet Now, for the problem at 
hand, the center of gravity of vorticity of the true sheet moves 
conically, but vorticity is being fed continuously into the sheet 
from the leading edge of the wing. Consider then the velocity 
of the center of gravity of the vorticity that has been shed for 


ward of f = 1, evaluated at the plane / = 1 

Let (Mo, 2*), (Ti, 21), and (T represent the strengths and 
centers of gravity of the vorticity at ¢ = ¢, of the vorticity that was 
shed prior to ¢ = 1, after ¢ = 1 but prior to/ = 4, and prior to 


{ = ft, respectively 
From conical considerations, the following relations are tound 


to be true 


m= tS l tly, | h — L)Peo 
Phen 
Voo* + (4h 1 Pot, 
lite: Pot, 
dz* dz, 
= 2isy — t ] 
dt, dt, 


Now as f; @ 1, 2) @ a, since the vorticity is shed from the leading 
cdge Hence, 
ie 


V(s* ,= 
" dt 


This form is, of course, identical to that given previously for 
V(sy) and would indicate that the assumption given here is equiva 
lent to making the vortex motion similar to that of the c.g. of the 
vorticity. 

For his assumption, Adams gives two values for the lift co 


eflicient for small angles of attack They are 
* 
C. = A.R.) a+ 1 + 4 A.R "=? 
: 2V2 
T i : 
C; = \ A R at 7 4A R a? 


For the hypothesis indicated here, the ambiguity is removed, 
and the single value obtained for the lift coefficient has the form 
. us i ‘ 

Cy = 5 (AR a + 7(A.R a 
The general form of the lift coefficient is similar to those given 


by Adams, but the coefficient of the a “* term 1s smaller than 


either of his values 
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Transport of the Boundary Layer in Secondary 
Flow 


F. F. Ehrich and R. W. Detra 


Aviation Gas Turbine Division, Westinghouse Electric Corporation 
November 6, 1953 


NOMENCLATURE 


1 = area = SS dodn 


A.R. = passage aspect ratio 


( = blade chord 

Ge = secondary flow energy loss coefficient 

n = length measured normal to the secondary flow stream 
lines 

q = perturbation velocity normal to the main flow direc 
tion 

R = passage radius of curvature 

= pitch 

t = time 

l = main flow velocity 

l = maximum value of the main flow velocity 

¥ = coordinate in pitch direction 

y = coordinate along the blade height 

€ = turning angle 

& = total turning angle 

a = length measured along the secondary flow stream line 


l = perturbation stream function 


QUIRE AND WINTER! have effectively described the perturba 
S tion velocity pattern in secondary flow through a bend in a 
rectangular channel or cascade blade passage. In terms of a 
perturbation stream function in the plane perpendicular to the 
direction of flow (in normalized form), 

Oy/Uos) 9 y/Uos) —«- HU/U) 


_ = oG (1) 


O(x/s)? O(y/s)? O(y/s) 


where «& is the total turning angle. A first-order estimate of the 
transport of the boundary layer may be had by considering the 
movement of the boundary layer due to this velocity pattern and 
neglecting the cross effect that the distortion of the boundary 


layer has on the perturbation velocities 


Actually, Eq. (1) is pertinent only to the stream function in 
the exit plane. But, if the assumption is made that the secondary 
velocities (and, hence, the stream function) build up linearly 
through the cascade, then, 

07(y/ Ups) : o(U/ Uo) 


O7(y Uos) 
4. = Ve (2) 
O(x/s)? O(y/s)? O(y/s) 


where ¢ is the turning in the cascade. Then the quantity 
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is invariant in the direction of flow 
Considering a particle of fluid in the cascade, its movement in 
the plane perpendicular to the direction of flow as it progresses 
through the passage is given by 
do =qdt } 
where g is the perturbation velocity normal to the main direc 
tion of flow. Then 


q = dy dn = elys|diy elys)/dn 7) 


where » is the direction normal to secondary stream line in the 


plane normal to the main flow. The time element can be ex 


pressed in terms of the radius of curvature of the passage, 
dt = (R/U )de 6 

Thus, 

dn \e de ] 


do = sR(U,/I d(y/els 


Considering a short fluid filament in the plane perpendicular to 


the main flow: 


SS Uo/s)d(n/s) = A/s? = 
SS (Uo CU) (R/s) 


where A is the area swept by the stream filament in the normal 


d(y eL/os) je de Ss 


plane 
Assuming that the filament maintains its mean velocity through 


the cascade, we find 


A/s? = (R/s) (Uo/U) (Al p/eUys)} (2/2) 9 


In the exit plane (e = «) we get, finally, 
A/s? = (c/s) (Uo/l A(y/eoUos)| (0/2) 10 
A very straight-forward graphical procedure is implied by Eq 
(10). Taking a short stream filameut between two stream lines 
of the secondary flow calculated from Eq. (1), the position of the 
filament in the exit plane may be located by positioning it such 
that the swept area is equal to the amount given in Eq. (10) as in 
Fig. 1. Fig. 2 gives the stream line picture in a semi-infinite 
(i.e., pitch/height ratio <1) passage for an inlet boundary 
layer whose thickness is equal to 0.2 times the pitch with a linear 
variation in velocity. 
The stream function was calculated by using an electrical 
analog for a coarse solution and refining the data by relaxation 
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OUTLET 


Fic. 5a (left Contour map of 6/6. Fic. 5b (right) 


methods. Taking 


(c/s) (6/2) = 0.375 (11) 
and plotting stream function at increments of 
A(y/eoUos) = 0.02 (12) 


it was found that Bernoulli surfaces would be calculated to cross 
one another, because the Bernoulli surfaces traveling at low ve 
locities are presumed to remain in the cascade for excessive 
amounts of time and, hence, are subject to excessive distortions 
In order to obtain a reasonable picture of the outlet flow then, 
we take 

U/Uo l (13) 
for purposes of Eq. (10) 

Fig. 3a gives the results of tracing the position of a large num- 
ber of points associated with various Bernoulli surfaces of the in 
let plane by the graphical technique. An interesting but some 
what deceptive comparison is given by the experimental results 
of Carter and Cohen? given in Fig. 3b. While Fig. 3a incorpo 
rates no loss phenomena and is merely the inlet boundary layer 
transported, Fig. 3b is an actual loss survey with the losses at 
each point being reported as the difference between inlet and out 
let total pressures (normalized on mean dynamic head). Hence 
the comparison in Fig. 3 is not quantitative but merely verifies 
the equivalent regions of secondary flow phenomena. 

Eq. (9) allows the graphical procedure to be applied to regions 
within the cascade and, hence, to trace the gradual development 
of the distortion of the boundary layer as it starts to turn as in 
Fig. 4b 

A possible origin of secondary losses might be the high-velocity 
gradients associated with sweeping high energy fluid onto the 
hub wall near the pressure side. 

If a measure of the shear stress is taken to be the vertical height 
of the boundary layer at any station, Fig. 4a gives a contour map 
on the floor of the cascade of the shear stress. The effect is ap- 
parently negligible on this basis of computation since the mean 
value of this shear figure on the floor of the cascade is practically 
equal to that of the undistorted boundary layer. Fig. 5 describes 
another possible mechanism of the secondary flow. 
a map of the kinetic energy of the secondary flow and Fig. 5b 
shows the map of the kinetic energy mass weighted with the axial 


Fig. 5a shows 


flow 





HUB WALL 


Distortion of the boundary-layer edge through the blade passage 


The energy loss coefficient resulting from the integrated mass 
weighted kinetic energy of the secondary flow (Fig. 5b) and non 
dimensionalized on the kinetic energy of the inlet flow is 


C. = 0.11786?/A.R.[1 — (0.2/A.R.)]* (14 


where A.R., the passage aspect ratio, = blade height/pitch 
This secondary loss coefficient is small for typical values of as 


pect ratio and turning angle 
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Remarks on ‘‘Interactions Between Wholly 
Laminar or Wholly Turbulent Boundary 
Layers and Shock Waves Strong Enough to 
Cause Separation’’® ! 


Seymour M. Bogdonoff 

Associate Professor, Gas Dynamics Laboratory, James Forrestal 
Research Center, Princeton University 

November 17, 1953 


N REFERENCE 1, Gadd presents experimental results and a 
I theoretical analysis of interactions between laminar or turbu 
lent boundary layers and shock waves strong enough to cause 
separation. For the turbulent boundary interaction, the prime 
justification for the theory is its approximate verification by his 
experimental studies. The writer would like to point out several 
relatively new sets of data obtained at the Gas Dynamics Lab- 
oratory which deviate considerably from Gadd’'s experimental 
results. References 2 and 3 present studies of the separation of 
a turbulent layer as caused by shock waves of varying strengths 


and steps of varying height. These investigations were carried 


* This research was sponsored jointly by Office of Naval Research, Me 
chanics Branch and Office of Scientific Research, under Contract No. N6 
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out in some detail at f ~ 3 and both show separation occurring 
at a wall static pressure ratio between 2.0 and 2.1, as compared 
with 2.5 predicted by Gadd. Of even greater importance are 
comparisons of the effect of pressure ratio on the “influence re 
gions,”’ Fig. 1 The vertical scale, the influence region (using the 
nomenclature of Gadd), is defined as the distance from the theo 
retical point of impingement of the incident shock to the begin 
ning of the interaction divided by the displacement thickness of 
The horizontal 
scale is the overall pressure ratio across the interaction 

A likely cause of the discrepancy may be the use by Gadd of 
shock generators that did not completely span the tunnel. This 
has been shown to have large effects on the phenomena under 
investigation.? Other results in reference 2 show that the as 
sumption of the basic velocity profile at the reattachment point 
is also questionable. The throw considerable 


doubt on the results of the analysis as applied to the prediction 
a shock 


the boundary layer ahead of the interaction. 


above results 


of the interaction of a turbulent boundary layer with 


wave. 
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Aerodynamic Determination of Flame Speed 


Robert A. Gross 
Harvard University, Cambridge, Mass 
November 6, 1953 


_ FUNDAMENTAL PROBLEM iti combustion studies is the de- 
termination of flame speed. Intensive efforts in the theoreti 
cal analysis of this problem and the experimental determination 
of flame speeds have been carried out. Important fundamental 
and practical difficulties still plague the theorist and the ex- 
perimenter. 

Rather than review here the extensive and rapidly growing 
literature on this subject, I refer the interested reader to two re- 
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cent reviews of the subject Here you can read of some of 
the problems involved and various attempts to overcome them 
It was from reading these reviews that I decided to make more 
widely known an entirely different method of measuring flame 
This method, based upon aerodynamic pressure meas 
urements, has been employed at the Harvard University Com 
Research Laboratory mid-1951 
difficulties in this technique, it 


speed 


bustion Aerodynamics since 
Although there 
has several distinct 


ment permitting the determination of local flame speed. The 


are, of course, 
advantages. It is a point-type measure 
physical quantities measured are different from those of other 
existing methods, permitting further checks on existing data 

In addition, information has been obtained on the total pres 
sure variation through the entire flame front 

It is first necessary to define what is meant by flame speed 
Here we restrict ourselves to a laminar flow of a homogeneous 
gaseous mixture of combustible material The flame speed S is 
defined as the dot product, 


where g is the velocity vector immediately prior (upstream) to 
the flame front and is the unit normal vector to the flame-front 
surface Now, what surface do we mean? We shall define this 
surface as that area where, upon approaching the flame front, 
we are first able to physically detect a change from the fre« 
stream temperature, composition See Fig. | 

In the temperature data that can be found in the literature and 


pressure, OF 


in Our pressure measurements, this point is well defined. Let 
us designate the properties at this point by the subscript | 
The element of gas thus proceeds from the free stream to the 
point 1, whereupon its properties start to change because of 
significant heat transfer, mass transfer, and chemical reaction 
Finally, thermodynamic equilibrium is established. Let us call 
this point of thermodynamic equilibrium point 2. The gas then 
flows away from the reaction zone. At atmospheric pressure, 
for common combustible mixtures, the distance from point 1 to 
point 2 is of the order of 1 mm, The volume enclosed betweer 


the surfaces 1 and 2 constitute the flame front 
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Since the flame front is so thin, it is reasonable to assume that 
it behaves as a discontinuity. We can therefore write mass and 
momentum conservation equations for a stream tube that passes 


through this flame front. Thus, mass conservation gives us 


where gq», is the normal (to the surface 2) component of the ve 


locity vector leaving point 2, and ) is the density ratio 
; I 


Momentum conservation normal to the flame front is given by 


wo 


Pi + piS? = po + pogen? (; 


where p is the static pressure. The tangential velocity remains 


unchanged; 
Gia ™ Ga * Wa (4 


We now introduce the total pressure p;, defined as 


Ph = prt pq” 2) (5) 


4 4), together with the definition of total 


Equations (2), (3), and 
pressure, permit the determination of the total pressure change 


across a flame front 
(A — 1 )p.S? 1 qa 2 2 
Pu —- Py = — 5 1+ " (“ (6) 


The flame speed may be solved for in terms of measurable quan- 


tities using Eq. (6). Thus, 


1 2 
S = i Date, as 1% Pew.) (7) 
“ee V,. iAop; pi 


where ;A.f; is the total pressure change between points 1 and 2 
and ,A,p; is the difference between total and static pressure (dy- 
namic pressure) at point 1. 

Eq. (7) is the basis of the flame-speed measurements to be pre- 
sented. It is based upon momentum and mass conservation 


across the flame front, and these appear to be reasonable assump 
tions 

Measurements have been carried out on propane-air flames in 
the Combustion Tunnel at Harvard University. A flat two 
dimensional V flame was employed. A small pitot-static probe 
was used to measure ;A,p;}. Radiation to the surrounding room 
was the chief agent that kept the probes from melting. The 
probes were rotated to obtain maximum readings at both points 
1 and 2 to eliminate errors because of yaw. A typical traverse 
through the flame front is shown in Fig. 1. A sensitive micro 
manometer, accurate to +0.0002-in. manometer fluid, was used 
to obtain pressure magnitudes. Values of the density ratio \ 
were computed from the best thermodynamic data available 
Using Eq. (7), values for the flame speed of a propane-air flame 
at atmospheric pressure are presented in Fig. 2 

Flame speed is a function of the initial thermodynamic state 
of the combustible mixture and the radius of curvature of the 
flame front.’ The initial state of our combustible is defined by 
the air-fuel ratio, the temperature, and the pressure. A flat \ 
flame was used to eliminate the problem of flame front curvature 

The experimental difficulties are primarily those encountered 
in accurate pressure measurements. Eq. (7) is extremely sensi 
tive to small errors in either pressure measurement and herein 
lies the most serious problem of using this technique. Great 
care must be exercised in making all measurements. Primarily 
because of the nature of Eq. (7) and the limit of the resolution of 
our measuring instruments, this technique is insufficient to dif 
ferentiate among flame-speed values given by other investiga 
It does, however, give additional flame-speed data ob 
Further information of 


tors.” 
tained by an entirely new technique. 
pressure measurements within the flame front itself can be read 
in a report by Quinn.* This work has been supported by the 


Office of Ordnance Research 
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Errata to ‘‘The Aerodynamic Forces on Low 
Aspect Ratio Wings Oscillating in an 
incompressible Flow’’! 


H.R. Lawrence and E. H. Gerber 
Cornell Aeronautical Laboratory, Inc., Buffalo, N.Y 


December 1, 1953 


bt HEADING: ‘‘Imaginary Part’’ in Table VI should read 
“Negative Imaginary Part.” 
The heading of Table VII should read 
‘Values of the Integral P?*(p, ¢ 


ind the second equation should be 


= bP*(p,q It ¢ L 
, k Jp \ 
The writers wish to thank H. M. Voss for pointing out the cor 


rections presented herein 
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Note on ‘‘Calculation Method for Three- 
Dimensional Rotationally Symmetrical 
Laminar Boundary Layers with Arbitrary 
Free-Stream Velocity... ”’ 


J. C. Cooke 
Professor of Applied Mathematics, University of Malaya, 
Singapore 


November 17, 1953 


T° THE SUBJECT PAPER, Drake! gives a method for rapid pre 
diction of the important coefficients in this case, based on 
Mangler’s transformation? and the two-dimensional Hartrec 
profiles. Seban* has deduced from these profiles that it is ap 


proximately true that 
(U,/v) (de?/ds) = B — A(U,'6?/y 


in Drake's and indeed the usual notation, and Drake gives B = 
144and A = 5.23 for 0 > m —0.09 and A = 5.32 for 0.068 
m>QO. Here we have written m = — U,'6?/v following Thwaites, 


to whose paper we now wish to refer By analyzing a large 
number of two-dimensional flows, Thwaites arrived at the con 
clusion that the coefficients A and B could be written 4 = 6, 
B = 0.45. These formed the best linear approximation to cover 


ill the flows whose solutions were known. Thwaites also stated 


that the work of Young and Winterbottom led to values A = 
6.28, B = 0.47 but did not consider these to fit the known flows 
so well as his own. It should be noted that the Hartree flows 


were, of course, included in the investigations of Thwaites. It 
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would probably therefore be better to take the values of Thwaite 
in preference to those given by Drake in his approximate method 
rhe value of m at separation is 0.068 for the Hartree profil 


be taken arbitrarily in his 


Thwaites, pointing out that this car 
ipproximate method, chose 0.082 for various reasons. Although 


of course, Seban’s approximation fits the Hartree profiles better 


solutions better 


than that of Thwaites, the latter fits all known 
ind should be used in the approximate method under consider 
it10n 

Drake then proceeds to make a test by comparison with Ho 


mann’s exact solution for flow at a stagnation point It is sur 
prisingly evident that it still does not seem to be generally known 
that Homann’s solution is in error, though this was pointed out 
by Mangler Again, Schlichting® pointed out that the Homann 


differential equation 


Qo tT <OO oO 


primes denoting differentiation with respect to the independent 
variable y’! could, by putting 
n= 2y, AG) = V 26 y 
be transformed into the Hartree differential equatio1 
: = Bf? — 1 
primes denoting differentiation with respect to 7) with 8 1/2. 


Schlichting also refers to some work on the subject by Fréssling 
As a matter of fact Miss Hannah’ solved Homann’s equation di 
rectly and obtained results different from those of Homann 
Cooke!’ pointed out the transformation and gave results sub 
stantially the same as Miss Hannah’s by interpolation from 


Hartree’s results The possibility of the transformation its in 


deed implied in Drake’s work, and he actually gives the value 
8 = 1/2in the course of his work Asa result of all this, we know 
that the coefficient in Drake’s Eq. (25) should be 1.312 


1.3172 


However the closeness of this agreement cannot be considered 


ind not 
, comparing closely with his approximate result 1.311 
as a test of Drake’s approximation because Homann's problem 
can be exactly reduced to Hartree’s problem Since Seban’s ap 
proximation was taken from Hartree’s solutions, we should ex 
pect the agreement to be close It would indeed be exact but 
for the fact that Seban’s approximation aimed to fit a// the Har 
tree profiles as closely as possible and not merely that for 8 
1/2. Still the approximation is bound to be close if Seban’s 
approximation is close, and the test is in fact a test of Seban’s 
work and not Drake's 

It may well be pointed out here that two errors have been de 
tected in Hartree’s table for y’ in his notation hey were found 
by differencing. They are as follows 
Bp = 1/2,x = 1.6, for 0.8860 read 0.8760 


B = 2,x = OS, for 0.7858 read 0.7958 


These corrections are made on the assumption that one digit 


only is in error in the values given 
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Errata—‘‘Quadruple Velocity Correlations and 
Pressure Fluctuations in Isotropic 
Turbulence”’ 


M. S. Uberoi 

Department of Aeronautical Engineering, University of Michigan 
Ann Arbor, Mich. 

November 23, 1953 


_ MISTAKES were overlooked in the original manuscript of 


my paper The correct formulas are 


In the ‘‘Introduction,’’ 


l Op 3 ol 1 , 
= l;? 0.18 K 
p* \On OX 


Selow Eq. (20), 


Iq. (22) should be 
(A\2/V"*) (dV'/dt)? = 25/R) 


Iq. (23) should be 


p? I (> 2 83h" 
\2(Op/dy)? Ay? RA? Ry 
Changes in Eq. (23) have negligible effect on the final calcula 


tions. 
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Further Remarks Concerning Integral 
Transforms of the Wave Equation 


Max. A. Heaslet and Harvard Lomax 
Aeronautical Research Scientists, Ames Aeronautical Laboratory 


NACA, Moffett Field, Calif 
November 19, 1953 


cpap ' HAS INVESTIGATED auxiliary terms that appear in 
an operational form of the potential equation of supersonic 
flow. The purpose of the present note is to supply a physical 
meaning to these terms and to indicate a similar case that has 
occurred previously. 
Consider the normal form of the linearized equation of steady 


state supersonic flow 
(07/0x?) — (07¢/dy?) — (0?¢/0z?) = 0 (1 


where g(x, y, z) is the perturbation velocity potential and the x 


axis is aligned with the direction of the free stream. The La- 
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place transform of any function F(x, y, 2) is 
F(s; { ail F(x, y, 2) dx 2 


e 


rhe transform of Eq. (1) is, after straightforward calculation 


OF dO¢ ide ; O¢ Of oe Of O¢ 
> ey eed A + A } A 
oy" Oz? = OX oy Oy O Oz 


where discontinuity surfaces y, 0, 1 , nN) exist 


in the flow field and the delta notation indicates the jump in tl 


value of the function at a discontinuity surface, the value Ag 


being assumed equal to zero. Considering, in particular, sur 


faces of discontinuity which are envelopes of Mach waves or 
characteristics of Eq 1), it follows that the bracketed terms i; 


the right member of Eq. (3) are directional derivatives of Ag(x 


y, along these surfaces. This direction, usually termed the 


conormal, lies in a plane determined by the free-stream velocity 


vector and the normal to the surface The directional deriva 


increments in tangential velocity and vanist 


tives represent 
from physical considerations rhe same condition follows for 
mally in linear theory from the assumed continuity of ¢, and th 
right-hand side of Eq. (3) can thus be suppressed 
We originally encountered this effect in the study? of unsteady) 
motion based upon the equation 
Org org ory Oy" 


Ox? ov Oz? of? 


Lagerstrom had suggested in his lectures at the California Insti 
tute of Technology that certain simplifications could be achiever 
in the study of wing plan forms with supersonic edges through 


the introduction of a transformed variable 3, where 


A(x, 2,t) = | Q(x, y, 2, t) ay 5) 


Let y = gi(x, 2, ¢) and y = g.(x, 2, t) be envelopes of the foremost 
discontinuity surfaces induced to the left and right sides of the 


Since ¢ is continuous across these surfaces, Eq. (4) be 


wing 

comes 

es daw 2p . O¢ Og 

eer ~ = S° (-1)'] 4 

Ox- Os* ot =f OX OX 

o¢ O¢ Og ¢ 0%; ; 
4 A 6 

oy Oz Oz Of Of 


The bracketed terms again are directional derivatives along the 
conormals to characteristic surfaces in four dimensions, and the 
continuity of g on these surfaces leads to the fact that the right 


side of the equation is zero 
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Approximate Solutions for Heat Transfer with 
Convection Flows 


M. Herbeck 
Max Planck-Institut fur StrOmungsforschung, Gottingen, Germany 


November 19, 1953 


NOTATION 


~_ 
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= temperature 
A = thermal conductivity 
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specific heat 
density 
kinematic viscosity 


absolute viscosity 


u 
a Prandtl Number 
coordinate parallel to the wall 
coordinate normal to the wall 
yy = stream function 
velocity of flow in x-direction: y, = p-u 
velocity of flow in y-direction Yr =p 


arbitrary constants 


1,4", B, D,a, a, p,? 


1) INTRODUCTION 
F: CERTAIN BOUNDARY CONDITIONS, Lighthill,! Schuh,” and 


Punnis*® have given solutions of the differential equation 


ple) Tes © ba Tx — b0:Te 11 
with 
y \ \ 1.2 
especially 
Y = const: y*:x% 1.3 
Eq 1.1) holds for the two-dimensional temperature field over 
a heated wall if the simplifying assumptions exist a) constant 


fluid properties, (b) negligible heat of friction and compression, 


ind (c) heat conduction in x-direction negligible against heat con 


9 


vection rhe application of the stream function, Eq. (1.2) or 
1.3), respectively, presupposes that the temperature decreases 
essentially in that part of the flow field where the velocity pro 
file may be approximated by the tangent at y = 0 

Now we will briefly demonstrate that the methods of solutions 


1.2) proposed in references 1 and 2 may be ex 


for Eqs. (1.1) and 
tended without any difficulties for the case that instead of Eq 


1” 


1.2) the more general function 


B-y "aa 


The special solutions for ¥ ‘ 
will be given in brackets below the more general solution 


with r > 2 is used 


Thus the heat transfer for separation flows, where r = 4, is in 


cluded 


2) SOLUTIONS OF THE DIFFERENTIAL EQUATION (1.1) FOR 


y =y r(x 


(2.1) Proposition of Solution with Exponential Function 


Introducing into Eq. (1.1) the new variables 


. 
mv 
t= Ux = r—1yf t(x ax 2.11 
o . 0 
Mh 
t ? a -B' we , 
o 
- pm Pg 2.12 
= (rp — 1)*/'.B y-x'* 
leads to 
I, = | 2 
with the particular solution 
», 
T = const-e~ *"- ‘Z) ( *) 2.14 
Z = Bessel function of the order +1 
lhe application of integral transformations to Eq. (2.14 


yields additional solutions of Eq. (2.13) which are not to be dis 


cussed here 


FORUM 14:5 


(2.2) Solutions with Confluent Hypergeometric Functions 


Solutions in the first quadrant which fulfill the following con 


ditions are of practical importance 


y= 0,% ) T(x, 0 l i 2.2 
1 x 1-a for B-x@ 
y => \ () T(x >i) Bi on 2 
X const 1 vr, ¥ 2.4 
decreases monotonously with y Using the proposition 
1 ee i 2.24 
1 A -x7-9(4 
with 
o a T r I 
(= B y-x® y Be 2 
u , 
yields 
T = A*-f ¢ ( ( 2.26 
where 
’ l | 
= ,Fi| p+ r 2.27 
F (» T © ( ) > IQ 
l 
p for 7 \ 
a r 
The boundary condition, Eq. (2.21), requires that ( | In 
order to satisfy Eq. (2.22), we apply the asymptotic expansions 
of Eq. (2.27) and (2.28), see reference 4, page 113, which yield 
that 
l/r)jtp! 
ese A a - 2.29 


is required if rp > 0 
iccording to the 
| 


Eq. (2.23), too, is satisfied 


Furthermore, 
theory of confluent hypergeometric functions, see integral rep 
resentation in reference 4, page 115 

=e > p 


and (2.23) are satisfied according to reference 4, pages 
§ pag 


1 — r, the conditions of Eqs. (2.21), (2.22), 


113-115 


by 
C<4<e@ 2.210 


However, C2 = Cy» yields the solutions that are physically plausi 
ble If0>rp > —1, only C C,; 


boundary layer 


gives a finite thermal 


] 
5 , I -d 
/V 
Bo-iors ge i- r,t Cy, yields the temperature field with 
the minimum 67 among the one-parametric family of solutions 
with Eq. (2.210), which tulfills Eqs. (2.21), (2.22), and (2.23 
together with a finite 67 If r-p l r, the region of Eq 
2.210) is reduced to C. (). For this interpretation, see refer 


ence 7, page 243 

If r-p <1 — r, the theorems for the zeros of confluent hyper 
geometric functions (reference 5, page 179) and the theorems for 
the zeros of solutions of differential equations of second order 
reference 2, page 211) yield that no solutions of the form of Eq 
2.24) exist which satisfy Eqs. (2.21 )-(2.28 

When 7-p > 1 — r, we get, with Eq. (2.29), for the temperature 


gradient at the wall 


2.211 
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- a + 
(Ty)y A-(p)-( B. - 
r 


We get a constant (7%,),~9, independent of x, if 


p—1 


U ”.(dt/dx) = const 


(3) APPLICATION 
In addition to references 1 and 2, which investigate solutions 
with ry = 3 for the conditions (2.21 )-(2.23), we calculate as ex 
ample with r = 4, the temperature field of a heated wall which is 
cooled by a flow with the stream function 
0.0375 p 
y = ° -D?2-x 1,18 
2-8 pv 
Eq. (3.1) is the first approximation for the stream function of the 
D-x 


2.29) 


boundary-layer profile with separation where U. = 7 
(reference 6, page 60). Applying Eq. (3.1) to Eqs. (2.26)-( 
we get 

.* (0.78: ft 


-{,Fi(1.65y + 0.75); 


C(p) 6-1 F\ (1.657 


with 


$)!'[1.65y — (1/4)]! 


$)'(1.65y)! 


Furthermore 


( D ) 
= A-C(p)-0.311-0 (3- 0-545 
' ‘ ; v(2 — B) , 


= ().545 


(Ty)y 


We compare the above approximation with computations of 


, is independent of x for 


Levy,’ who has given solutions for Eq. (1.1) by means of nu- 


merical integration where, instead of Eq. (3.1), the stream func 


tion of the exact solution of the boundary layer equation for 


separation is applied. Levy calculates, instead of Eq. (3.5), 
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the ‘‘reduced”’ temperature gradient (@,), which is defined by 


T(x, y) = A-x?-O(n 


t 


) -x~0-5.y with B = —0.198 
: : 


la: 
v(2 — 


Applying Eq. (3.1) yields 


(Oy)n-o = C(p)-0.311- 3.8) 


Values according to Eq. (3.8) and “improved” solutions of 


reference 7 are represented in Table 1. They indicate only few 
per cents of difference. Thus the application of the approxi- 


mation Eq. (3.1) into Eq. (1.1), which allows a simple computa- 
tion of the heat transfer, seems to be as useful as corresponding 
2 is chosen 


calculations for r = 3 in references | and 2, if ¢ > 1 


according to the simplifying assumptions 


TABLE | 
Values of (@y)n—0 for o = 1 
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0 0 
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—() 401 —0 
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